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OVERVIEW 


In October of 1976 Columbia University undertook an investi- 
gation of a class of space vehicle control problems relating 
to vehicle or component flexibility. The sole Principal 
Investigator initially was Peter Likins, Professor and Dean 
of the School of Engineering and Applied Science. In August of 
1978 Professor Richard Longman joined the project as Co- 
Principal Investigator. Under the direction of the two co- 
principal investigators, the substantial tasks of this study 
were largely the responsibility of Mr. Chittur Viswanathan, 
and the report that follows is the principal content of his 
Ph.D. dissertation. Appendix B of Part I of this report (see 
Section 8) is a technical paper written by Professor Longman 
for journal publication, and as such it may be read independently 
of the body of the report. 

The study progressed in two phases, represented by Parts 
I and II of the Final Report. Part II is submitted here in 
more abbreviated form than is planned for the Viswanathan 
dissertation, both because of the size of the total docu- 
mentation package and because of delays in preparation of the 
report. The expanded presentation of Part II from the dissertation 
includes an appendix of derivations, available upon request. 

In technical content the study addresses two quite distinct 
engineering questions: 

Part I . Where should actuators be placed on spacecraft 
with distributed flexibility? 
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Part II. How should we specify the dynamic characteristics 
of flexible instruments to be pointed by a prescribed control 
system? 

The class of problem considered in Part I is of wide-ranging 
concern, since it is raised by many proposed spacecraft; the 
Solar Power Satellite is perhaps the most dramatic example. 

The class of problem treated in Part II is of specific 
applicability to the Instrument Pointing System on the Space 
Shuttle. In this application an automatically controlled gimbal 
mechanism is designed to point a variety of interchangeable 
instruments toward their targets, with the gimbal system mounted 
on an actively controlled spacecraft. Future instruments must 
be designed after the control system characteristics are estab- 
lished. 

Parts I and II of this report describe work performed in 
reverse chronological order. The task initially undertaken 
(Part II) was an attempt to apply the methods of parameter plane 
stability analysis to the Shuttle-based Instrument Pointing 
Mount (IPM). By establishing stability boundaries in the plane 
defined by two design parameters of the pointed instrument, 
one might hope to provide the instrument designer with a useful 
tool for the selection of these two design parameters within 
the range for which the control system is stable. This concept 
(originally suggested by Dr. Sherman Seltzer, then of the NASA 
Marshall Space Flight Center) is easily implemented for suf- 
ficiently simple mathematical models of the spacecraft system, 
and the resulting parameter plane plots (see Part II) offer the 
prospect of genuine utility for preliminary design. However, 



attempts to extend the useful range of the method to realistically 
complex mathematical models are less successful. The parameter 
plane approach can still be applied to system equations of con- 
siderable complexity, but the parameters indicated for the 
portrayal of stability boundaries are less useful for practical 
design. 

After exploring the parameter plane concept for the IPM to 
what seem to be its practical limi cs, the Columbia group turned 
its energies to the actuator placement problem. The results of 
this effort (Part I following) are very encouraging, although 
not yet definitive. 

In the design of the spacecraft of the past and present, 
the question of actuator (and sensor) location has received 
remarkably little attention. It has generally been assumed 
that both sensors and actuators should be placed on the 
essentially rigid, central body that comprises the core of most 
spacecraft of this era. In those exceptional cases in which 
sensors and actuators have instead been attached to flexible 
appendages, the objective has been to address specific 
physically defined requirements. For example, gas jets were 
located at the tips of the solar panels for the Mariner space- 
craft series in order to maximize the "lever arm" producing a 
moment about the mass center of the spacecraft; the f Iv^^xibil ity 
of the solar panels was ignored in this sele-.tion. 
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Future spacecraft may differ significantly from their ante- 
cedents in both scale and the distribution of flexible components; 
in some cases it may be necessary to treat the entire spacecraft 
as a flexible body, with no rigid, central core. In such cases, 
it becomes quite unclear where one should locate sensors or 
actuators, or even how many should be employed. If one compares 
two alternative designs, it is not even clear what criteria one 
should employ in comparisons. 

The chief contribution of Part I is the development of a 
criterion for the comparative evaluation of alternative actuator 
locations. This criterion is a measure of a quantity introduced 
here as the "degree of controllability." 

In what follows a definition is generated for the "degree 
of controllability" concept (after alternative definitions are 
considered and discarded) ; techniques are established for 
calculating this measure; and application is made to flexible 
spacecraft, both genetically and specifically. 

Although it would be unrealistic to argue that the criterion 
advanced here provides a unique measure of the quality of an 
actuator distribution design choice, it does seem to be as good 
a basis for evaluation as any measure yet devised. 
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PART I 


THE DEGREE OF CONTROLLABILITY CONCEPT 
AND ITS APPLICATION IN THE LOCATION OF 


ACTUATORS ON VERY LARGE FLEXIBLE SPACECRAFT 




r 


1. INTROOXTION 



This rasearch addresses problems broadly connected with the 
dynamics and control of flexible spacecraft. From the dynamics point of 
view, the modeling of a spacecraft consists of treating some components 
as rigid bodies and some as flexible. Thus, solar panels, booms, etc. 
must usually be modeled as flexible, but they may have a rigid core as a 
base. The sizes of the spacecraft can vary tremendously according to 
their intended purpose. For instance, one design for a solar power 
station (SFS) satellite calls far a 4.9 km x 14 km solar array. In such 
applications, it will be necessary to control not only Lhe attitude, 
i.e., orientation, but also the shape of the various components of the 
spacecraft. The need for shape control arises due to the flexibility of 
the spacecraft or its components. Stringent requirements may be imposed 
on the accuracy of attitude and shape control. For instance, accuracy 
of the order of thousandths of an arcsecond are sometimes needed in the 


attitude control of certain scientific instruments. On some future 
spacecraft which function as large antennas, accurate shape control will 
be necessary to obtain distortion free signals. The accuracy 
requir&nent is a function of the wavelength of the signal transmitted 
from the antenna. In order to achieve this goal for large spacecraft 
with distributed flexibility, sensors and actuators must be distributed 




Proa the control point of view, spacecraft can often be classified 
into tvo loain j roups: 

1. those for vyhich ^e spacecraft is fully defined before the 
control system is designed; and 

2. Those £ar which the control systam must oe specified before 
certain interchangeable parts of a multi-purpose spacecraft are 
selected for future missions. 

In '^hat follows consideration is given to both classes of problems. 

Spacecraft in the first class includes tiie usual case, and the 
design of stable control systems for conventional spacecraft is a mature 
discipline. However, the extreme size and distributed flexibility of 
many future spacecraft pose various challenging new problems in their 
attitude and shape control. One such problem is the development of a 
rationale for the distribution of sensors and actuators throughout the 
entire body of the spacecraft. How to compare different actuator 
distributions and vliat criterion to apply to be able to decide on a 
particular choice of distribution of actuators are the specific 
objectives in Part I of this work. 

^ typical mcample in tha second class is the Space Shuttle's 
Instrument Pointing System (IPS) which must be designed to accurately 
control the pointing of any of a family of scientific Instrunents, some 
of which have yet to be designed. These instrvments are typically 
lightweight and their flexibility effects are important due to pointing 
requirements. The actual instruments involved and Lheir purpose will 
vary from one shuttle flight to another. A spacecraft compatible with 
'ny one of such a set of instriments becomes a multi-purpose spacecraft. 
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9ut it is not possible or practicable to design a control systam for 
I c every instrument the spacecraft night carry, since some of these 

instruments may have to be planned for future needs. What is done, 
instead, is to predesign a coraprehensive control system for attitude 
control of a rigid spacecraft. This control system is specified for the 
corresponding multi-purpose spacecraft. And tiie objective is to seek 
the characteristics in the most general terns of those flexible 
. ^ instrixnents the given control system can control with appropriate 

stability margins. 

This research vork is organized in two parts. ?art I deals with 
spacecraft belonging to the first group, in which the system is fully 
known and the distribution of actuators in the design of control systems 

is the subject. Part II deals with spacecraft belonging to the second 

• r 

group, in which the system is not fully known and tdie characteristics of 
the class of flexible instrunents compatible with a given control system 
are investigated from the stability point of view. 
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PART I 


THE DEGREE OF CONTROLIABILITY CCNCEPT ITS 
APPLICATION IN THE LOCATION OF ACTUATORS 


ON VERX LARGE FLEXIBLE SPACECRAFT 



2. THEORY 


2. 1 Introduction 




In Part I attention is focused on an aspect of design of the 
control systems vihose purpose is to control effectively the attitude and 
shape of very large flexible spacecraft. In the past, the approach to 
the control system design wes usually based on the philosophy that a 
spacecraft was essentially a rigid body with Lhe flexibility of 
antennas, booms^ solar panels, etc. treated as extraneous disturbances. 
The purpose of the control was only attitude control of the spacecraft 
and the design vas based on the rigid spacecraft with adjustments made 
for the influence of the attached flexible components on the behavior of 
the entire spacecraft. The term Dynamic Interaction was then used to 
describe this influence due to the component vibrations and localized 
energy dissipation, which were treated as second order phenomena even 
though occasionally they proved to be destabilizing for the entire 
spacecraft. A good historical review regarding the prevailing attitudes 
in this field for the last couple of decades can be fbisid in tH* 
Tbday, the planners in the spacecraft industry are contemplating 
structures of huge dimensions (of the order of kilometers) for future 
spacecraft of the 1990 *s or 2000 and beytnd. For these vehicles, 
flexibility can no longer be considered a second order phenomenon and 
must be treated as an intrinsic property cf the entire spacecraft. In 
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designing control systems for these spacecraft the objective is no 
longer just the attitude control but simultaneously the shape control of 
these huge structures. By shape control we .mean the restoration of 
these huge structures to their nominal shape vi^enever they are disturbed 
due to any causa { thermal stresses, gravity gradient torque, aerodynamic 
torque, solar wind, etc. etc.). This is done by controlling the 
vibrations of a finite number of modes of a truly infinite set of modes 
'.lihich describe the physical shapes of these distributed parameter 
systems. Because of the strong coupling of the attitude iynamics and 
vibration .modes, their control must be treated simultaneously. 
Classical control theory is unsuitable for such large dimensio.nal 
multiple input multiple output systems, and it becomes essential to 
adopt modern multivariable control theory techniques, 

^ In order to achieve attitude and shape control it will be necessary 

to distribute actuators throughout the entire body of the spacecraft. 
How should the number and locations of actuators be chosen in order to 
best control the flexible spacecraft? This problem has been recognized 
for some time, but to date little has appeared in the literature that 
would help guide the control system designer in placing the actuators. 
>tost of the known results id(3itify the minimun number of actuators 
needed far a given set of modes to be controlled, and identify certain 
specific actuator locations which cannot be used because they result in 
an uncontrollable system. 

Cnee the designer chooses a set of actuators it is necessary to 
make sure that whatever distribution pattern .he chooses t.he system is 
controllable. The concept of controllability in modem control theory 

« 

sA , 
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is a binary concept, either a systsa is controllable or it is 
uncontrollable. Starting from a set of actuator locations vrtich produce 
an uncontrollable system, but for vihich the nunber of actuators is 
sufficient to produce controllability, it will usually be the case that 
moving one of the actuators by a distance S > 0 can produce a 
controllable syst&\, no matter how small the £. One expects that for a 
small € , even though technically the system is controllable, in some 
sense it will not be very controllable. It then seems natural to ask 
"how controllable is the system with a particular actuator 

distribution?" It is, therefore, reasonable to seek extensions of the 
established definition of controllability so as to permit a precise 
definition of the "degree of controllability" siiich would prove useful 
for actuator placement. 

It is the purpose of this work to generate, starting from basic 
physical considerations, a rational definition of the degree of 
controllability. The definition obtained is certainly not the only 
possible definition, but * . does have the advantage over a definition 
based on singular value de jmposition that the physical reality of 
actuator saturation limitations is included in a fundamental way. 

The definition is then applied to the actuator placement problem 
for flexible spacecraft. With this tool the control system designer can 
rank the desirability of various candidate actuator distributions, and 
thus he would have a systematic way of pickitig which distribution to 




use 



2.2 Definition of the Degree of Controllability 


Let us consider any general linear time invariant systen in state 
variable form 


x*(t) «c/4x*(t) +^u^(t) (2.2-1) 

vhere xcR and ueR***. It should be noted that although we focus our 
attention on this system, the degree of controllability definition vhich 
we adopt is also applicable to more general systems of the form x*(t) » 
f(x*(t) ,u*(t) ,t) having a solution x*(t)»0 (f{0,0,t)»0) . 

Several candidate definitions were scrutinized in the course of the 
search for a suitable definition of the degree of controllability. It 
is instructive to discuss some of tliese candidate definitions which were 
considered and discarded — the process of starting \7ith a blind attempt 
at a definition and progressing to a wall formulated concept highlights 
the characteristics that a workable definition must have. We discuss 
these candidate definitions in the following: 

1. Eigenvalues of ^ 

It is tempting to try to connect the degree of controllability 
to properties of the standard controllability matrix 

Q - «4iB .4V ), and define degree of controllability as the 

square root of the minimun or maximun eigenvalue of QQ^’ . Five apparent 


^ The superscript T is used to denote transpose of a matrix 
or vector throug)»ut this text. 
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difficulties with this definition .*nust som^ow be handled before the 
definition becanes viable. These are: 1) The degree of controllability 

is affected by a transformation of coordinates (since the eigenvalues of 
33 '*’ are not invariant under changes in state variable representation) . 
2) Although this candidate definition satisfies the basic requirement 
that the degree of controllability is zero when the system is 
uncontrollable, it is not immediately clear what other physical meaning 
can be attached to 02^ and hence to the size of its eigenvalues. 3) It 
is not clear how the stability of the system is reflected in this 
definition. It is much easier to control a stable system with the 
objective of returning the system toward the origin (x**0 solution) than 
an unstable system with the same objective. (The Reverse is true when 
the objective is to reach out from the origin.) -i) The candidate 
definition does not involve a dependesice on the amount of time T alloted 
to acconplish the control tasic. It can be much easier to control the 
system state in some directions in the state space at one time than at 
another time, so the degree of controllability should depend on T. 5) 
It is not clear that the amount of control effort needed to accomplish 
^ the control tasl< is reflected in this definition. In the satellite 

described in section 2.1 where one actuator has been moved by a snail 
amount l to produce controllability, one expects the ’’weak 
f controllability" of the system to be manifested in the need for very 

large control actions to accomplish certain small changes in the state, 
and hence the control effort required is of fundanental importance in 
^ making a definition. 



2. Impulse response 

It is clear that some type of limitation or standardization of 
the control effort must be included in the definition. Consider a 
standardization vhich restricts the control to a tsiit impulse / and 
consider systems with «4 in diagonal form and with u* a scalar. For 
distinct eigenvalues the system is controllable if none of the elements 
of the column matrix 3 are zero. Furthermore, these components 
indicate how far a unit impulse control will move each state component 
instantaneously (impulse response) , so one might suggest the m^ { icl as 
a degree of controllability. Here we are trying to generalize a second 
standard test for controllability to obtain a degree of controllability 
definition. The difficulties with this candidate definition are: 1) 

The control actions are so restricted that the components of the state 
cannot be affected independently. The control of all states by a single 
control u'*’ relies on the differences in the dynamic behaviors of the 
states. 2) The candidate definition does net involve a dependence on 
the eigenvalues of the system, wliich means the information regarding the 
stability of the system is absent. The rest of the argument is the same 
as in (3) of candidate definition 1. 3) The candidate definition does 

not involve a dependence on T and exactly the same argiment as in (4) of 
candidate definition 1 holds. 

3. Qiergy due to an impulse 

When an impulse is applied to a physical system there occurs a 
change in the total energy of the system, and one might vender if this 
change in energy due to an impulse could not be used to define the 
degree of controllability. Since this is dependent on the impulse input 
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it ivlll suffer from all th« dravfeacks of tha impulse response discussed 
under the candidate definition 2. 

Itiese candidate definitions do not seem to include the effects of 
all pertinent variables. Hence, it will be necessary to build the 
definition from more fundamental considerations. It is interesting to 
note that, in certain special cases, the resulting degree of 
controllability definition will be a modified version of the candidate 
definition 2 (and by employing a different approach involving singular 
value decomposition of matrices something of the general form of the 
first candidate definition can result) . 

It is now evident that the definition of the degree of 
controllability, besides being in some sense a measure of hov^ easy it is 
for the controller to control the system, must in some way handle five 
things: 


1) It must have the property that the degree of controllabliity is 
zero when the system is uncontrollable. 

2) It must reflect in some way the stability information of the 

% 


system. 

3) It must som^ow consider dependence on total time T. 




4) It must standardize the control effort in some my, 

5) The control objective must be restricted. 

Concerning the last point, certainly different control objectives should 
influence the choice of the control system design, and hence the degree 
of controllability of a candidate design should be keyed to the 
objective involved. In a large class of problems (regulator problems ) , 
the equilibriun solution x*»0 to equation (2.2-1) is of primary 
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importance, and the control objective is to return x^to zero after a 
disturbance. Since this is the most common attitude and shape control 
problem for flexible spacecraft, we will restrict ourselves to this 
objective. Concerning the standardization of the control effort we will 
require that the control components satisfy |ut,)4 I for i-l,2,...,m, 
v^ich represents realistic physical limitations of the actuator 
capabilities. Note tliat the use of one as the boind for all control 
components implies normalizing each component of u* to produce a new 
•control vector u, and adjusting the matrix ^ to produce a new matrix B. 

Controllability requires the existence of a control function vhich 
can transfer any initial state to any final state in finite time. With 
our more limited control objective, the degree of controllability should 
be related to the volume of initial system states (or states resulting 
from disturbances) vghich can be returned to the desired state x*=0 in 
time T using the bounded controls. Consider the nature of this volone 
in more detail. In a controllable system, more initial states in any 
direction can be returned to t!ie origin if the system is stable than if 
it were unstable. Bence, the volume of initial system states is greater 
in the case of a stable controllable system than in the case of an 
unstable controllable system. Now, in an uncontrollable system there 
will be at least one direction in the state space for which initial 
conditions in this direction cannot be returned to die origin, and the 
volune will lose one or more dimensions, fbr a controllable system 
wiiose paraneters are such that it is nearly uncontrollable, only initial 
conditions very close to x*»0 along the above mentioned direction could 





be returned to the origin in time 7 using the bounded controls. Hence, 
we will generate a definition of the degree of con':rollability based on 
the minL'nua distance from the origin to a normalized state that caanot 
be brought to the origin in time T. *4ore loosely it is the minimun 
disturbance from which the system cannot recover in time T. 

The coordinates of a state space will very rarely all have the same 
physical units, and hence it is clear that comparing distances in the 
state space will require that each coordinate must be made uiitless by 
normalization. How should one choose the normalization to use? 
Recognize that whan comparing two controller designs for controlling the 
same dynamic system, the needed minimun distance for each design will 
usually correspond to a different direction in state space. Hence, 
ranking of the degree of controllability of the two systems will depend 
on comparison of distances in different directions, and this implies 
that we must be equally interested in controlling deviations of the 
state from x’’»0 in all directions in the state space. in order to 
accomplish this the control system designer must specify n-l nunbers 
which represent his degree of interest in controlling each component of 
the state. This could be done, for example, by datermining the 
daviations of x* , ,..., x*_, which would be considered of equal 

importance to a deviation of x* The reciprocals of t.hese nuabers 
would then be used to produce normalization factors for each of t.he 
coordinates of the state space giving a new state vector x. The system 
equations expressed in terns of the normalized state x and normalized 
control u are then written as 



X(t} ^ Px(t) + 3u(t) 


( 2 . 2 - 2 ) 



lu.J < 1 

Just as in optiaial control viihere the control systen designer 

must be specific about his goal by specifying a cost functional, in 
order to define bhe degree of controllability, it is necessary to be 
fully specific not only about the objective of keeping 2 C» 0 , but also 
about the relative importance of keeping earu component of x near zero. 

Relative to the normalized system (2.2-2) we are now ready to make 
the following definitions: 

Definition 2.1; The recovery region for time T for normalized 
systan (2.2-2) is the set 

^ » ( x(0) I 3 u(t) , t« C0,T], luL(t)I < 1 for 

i*l,2,...,m 9 x(T) » 0 } 

Definition 2.2: The degree of controllability in time T of the x*0 

solution of normalized system (2.2-2) is defined 
as 

P » inf II x(0) I! -9^ x(0) ^ 

where fl*l) represents the Euclidean norm. 

Ihus, the recovery region identifies all of the initial conditions 
(or disturbed states) vhich can be returned to the origin in time T 
using the bounded controls. And the degree of controllability is a 



scalar measure of the size of tihe region, ’.^cre the scalar is chosen as 
the shortest distance from the origin to an initial state which cannot 
be returned to the origin in time T. 

The degree of controllability, as defined, is keyed to the state 
vector X employed. tto transformations of coordinates can be allowed 
once the normalization has been ^cified (unless the norm used in the 
definition is adjusted to compensate for the resulting distortion of the 
state space) . 

It should be pointed out that althoiqh Definition 2.2 incorporates 
all the properties viiich were identified as necessary in the definition 
of the degree of controllability, it is not necessarily unique in doing 
so. For example, a standardization of bhe control effort in terms of 
energy can also be employed, but the inec^jality saturation constraints 
on the controls used here represents the more realistic situation. 



2.3 Concepts for Approx iaatin-i the -tecovery Region 

In order to ,'nake the definition of the degree of controllability 
useful, it is necessary to develop a simple algorithm to generate at 
least an approximation to the distance (* . This necessitates 

approximating the recovery region 

Note that the solution of (2.2-2) is given by (see [2], [3D 

t 

x(t) -^(t,t.)x(t.) ^-i(t,t.) J §(t,,s)Bu(s) ds (2,3-1) 

vhere f(t,t«) is the state transition matrix for (2.2-2). A complete 
explanation of f can be found in [2]. Without any loss of generality v« 
can assixne the initial time t,«0 and the final time t>^. Also *we will 
denote x(0) by x, and x(T) by x^. The solution (2.3-1) then is simply 

x^-f(T,0 )X, +i(T,0) I $(0,t)Bu(t) dt (2.3-2) 

0 

r.ie displacement in the state space in time T is x^- x^, and let 

- x^ - X. (2.3-3) 

and the distance moved during time T is the Euclidean norm || Prom 

(2.3-2) the displacement 6 can be expressed in either of the following 
form: 

T 

c E - §"\t,0) ]x^+ J^(0,t)Bu(t) dt 

0 


(2.3-4) 
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« tf(T,0) - S ]x, +t(T,0) «(0,t)Bu(t) dt 


vdhere E is the unit matrix and the subscripts on o are used to identify 
the two forms. In (2.3-4) the displacement is expressed in terms of the 
fir»?l state x^ and in (2.3-5) it is expressed in terras of the initial 
state X,. Our interest lies in in which ;<^»0 and ws termed tiie set 

of states X, as a recovery region (see Def.(2.D) . A companion region 

which is the set of states x^ for wtiich x^«0 can be defined in a 
similar way to the recovery region The region a then represents the 

f* 

reachable states from the origin. Putting x^»0 for a, , and :<^«0 for 
in (2.3-4) and (2.3-5), respectively, we can obtain 


#(0,t)8u(t) dt 


(2.3-S) 


I 

-^(T,0) j $(0,t)Bu(t) dt 


(2.3-7) 


where the (*) on (i«l,2) has been remaved to indicate that the 
displacement is about the origin (x^ or x^ is zero) . 

(tote that for linear time invariant systems the fundamental matrix 


can be written as 


8 e 


A(t-r) 


(2.3-3) 


from -.liiioh 


If 


AT / \ -At 

‘§(To)s e e (2.3-9) 

Two remarks are worth making regarding tdie comparisons laetween the 
regions A and corresponding to the displacements and 
identical control u(t) is applied during the same time T for both cases* 

1. The region .A. of the reachable states from the origin in time T 
is the same as the recovery region A if the control process 
were taking place in a time domain where time runs backward 
(interchange t^(^) and T in and obtain the same magnitude 
as of ) , 

2. The regional is largt than the region A for tLme T if the 

system is unstable about the origin (i.e., the real part of at 
least one eigenvalue of A is greater than zero) . This is 
because there is always at least one component of ^ whose 
magnitude is greater than the corresponding component of , 
(due to the influence of the exponential terms, , in 

§(T,0)) * If the system is stable the converse is true. 

In the following analysis we will be concerned only with the 
recovery region A and the corresponding equation for the state space 
displacement is (2.3-6)t 


‘^The superscript 1 will be omitted henceforth and d' , 
unless otherwise mentioned, corresponds to che recovery 
rogion A , i.e., in (2.3-5). 
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3y ths Caley-Hamilton theorem the state transition matrix can be 
written as 

n- I 

f(0,t) - (2.3-10) 

•<«0 

where the are scalar functions of time. Partition the B matrix into 
coltran matrices t»j, and define the following matrices 

B » [ b, b, b^ ] (2.3-11) 

• [%W V. J (2.3-12) 

* [b As (2.3-13) 

A bp] (2.3-H) 

Than ^ can be represented in the following alternative forms 

i -111 ,,. 3 - 15 , 

w\ ^ 

m ^ ifYU tMl'. 

* ^ ) L^o bp T -r, p *J»-I P (2.3-15) 

P*1 0 

• Z J (2.3-17) 

P'l 0 ^ 


Ebr Lhe purposes of illustrating certain concepts, let us restrict 



ourselves to the case of a scalar control so that the surcroations over f 
as Mill as the subscripts in the above can be dropped, and 3 is a 
colunn matrix b. Also let n>2 for simplicity. Suppose the recovery 
region is as shown by !tegion I in Fig.l. The maximum x, component of 
any state in the recovery teg ion is obtained by using xhe control u 
equal to minus the signun function of the first Cixnponent of the vector 
[qy] in (2.3>17), since this maximizes the x, component of the integrand 
at each time t.lhe right hand side (and left hand side) of the rectangle 
enclosing this recovery region in Fig. 1 can thus be found by 
integrating the f:rst component in (2.3-17) using this control. If 
desired the point at which th'^ recovery region touches this side is 
obtained by integrating the second component of (2.3-17) using tiiis 
control. Ihe top and bottom of the rectangle are found similarly. 

The rectangle obtained in this manner might be considered an 
approximation to the recovery region, and then the shortest distance 
from the origin to one of the sides might be considered an 

A 

approximation, , to the degree of controllability, ^^3t 

this necessarily produces a ^ which is an upper bound for the degree of 
controllability. In some cases this approximation is a tight one, but 
often it is not. Suppose the recovery region was Region II of Pig. 1. 
This corresponds to a system which has a much poorer degree of 

A 

controllability, yet the approximation f remains the same. In 

fact, suppose that ^ Region II degenerates to a 

line forming a diagonal of the rectangle in Fig. 1. Then the system is 

A 

an uncontrollable system, but the approximation (* still predicts a 
good degree of controllability. Hence, this approximation must be 



rajectsd. 

For the case of a scalar control being considered# this shortcoming 

can be eliminat{ 2 d by using ^ as expressed in (2.3-15) and maximizing 

.•r, 

components along AW. The control 

aCv) s (2.3-18) 

eactremizes the coefficient of the vector A W in (2.3-2). It will 
simultaneously produce some components along &e other vectors A^W 
for / . This is a maximization of a component of the vector S but 

it is a component as seen in a nonorthogonal set of coordinates. Bence# 
the upper bounds obtained in the various directions define a 

parallelogram (more generally an n dimensional parallelopiped) which can 
be considered as an approximation to the recovery region, as shown in 
Pig. 2. As before there is some point on each side of the 

parallelogram which is in the recovery region, but no point outside the 
parallelogram is in the region. 

Die minimuR distance to a side of the parallelogram, l.e., the 
minimum perpendicular distance to a side, is an approximate n f to the 
degree of controllability, f . When the system becomes uncontrollable, 
the columns of Q become linearly dependent, and hence the perpendicular 
distance to one of the sides becomes zero. Ihis means that this 
has the essential property that (***0 whenever (^«0. 

We conclude that for the scalar control case we have a viable 
method of approximating the degree of controllability. A simple method 
will be presented in a later section to detarmine the needed minimum 




L*^ 


perpend j-jular distance. 

“Riis approx nation is still an upper bound, and it can be improved, 
in fact .made arbitrarily good, by considering more directions in the 
state space. Let e be any desired unit vector expressed as a column 
matrix of components. By examining (2.3-17) the state in the recovery 
region having a maximixn component along the direction e is obtained 
using tte control 

M. - - s^»v ^^4^3 (2.3-19) 

and hence no points in the recovery region lie beyond the line 

perpendicular to e and a distance 

T 

/ I ® (2.3-20) 

from the origin (but at least one point in the recovery region lies on 
the line). Figure 3 illustrates how use of three e's (e,, e^ , and 
e^ ) identifies three tangents to Lhe recovery region, and vhen taken 
together they begin to approximate the region boundary. Let ^ be 

the minimun value of (2.3-20) for any set of directions e considered. 
Then an Improved estimate of the degree of controllability is 

■®**»rain (^* ^ , and (*** can be made arbitrarily close to the true 
degree of controllability (^ by picking a sufficient nonber of 
directions e. This method of improving the approximation to the degree 
of controllability will also be generalized to the multiple control 


case 



2.4 fundamental Squation for Recovery Region Approximation 
in the multiple Control Case 

Ihe previous section presented a procedure for generating an 
approximation (®* to the degree of controllability C in the case of a 
scalar control u. The procedure required the use of n carefully chosen 
directions in the state space, b,Ab,,.,A b, in the approximation to 
the recovery region in order to insure that (®* had the property that 
f*«o if and only if the system is uncontrollable. If the control 
vector is m dimensional with »m>i it is no longer obvious how to obtain 
this property, since the colunns of the Q matrix necessarily contain 
linearly dependent vectors. Instead, we will consider the eigenvectors 
and generalized eigenvectors of the A matrix of (2.2-2) as the n 
linearly independent directions in the state space. Certainly some 

modifications must apply when these vectors are complex. In the single 
control case the value of became zero when the system became 

uncontrollable because linear dependence of the vectors 
b,Ab,....,A b implies the collapse of at least one dimension of the 

parallclopiped . The vectors chosen here for the multi-dimensional 

control case do not exhibit this reduction. Nevertheless, it will be 
shown in the next section that the desired property of the resulting 
C can be demonstrated under fairly general assvinptions. This 
section is devoted to generating bie appropriate expression for S 
equivalent to equations (2.3-6, 2.3-15-2.3-17), expressed in terms of 
components in these eigenvector directions. 

Let J be the Jordan canonical form of the matrix A, and let P be 



( 2 . 4 - 1 ) 


the matrix of eigenvectors and generalized eigenvectors so Lhat 
P"'.\P»J; J » diag r J, , J^] 


vyhere the 7^ are the square Jordan blocks of dimension , 
k>l,2f ....,r. A diagonal Jordan block is of order one. Associated with 
each Jordan block is an eigenvalue so that Y§n and r is 

greater than or equal to the number of distinct eigenvalues. Also, let 

A T 

Pj be the n colimns of P, and be the n rove of P (the left 

eigenvectors and generalized left eigenvectors) . Every Jordan block 
TJj corresponds to one independent eigaivector of A and 
generalized eigenvectors, and the same is true for Lhe left eigenvectors 
and generalized left eigenvectors. 

It is appropriate to list here a few results on matrix algebra (see 
Appendix A for derivation) vhich we vail use in our further analysis: 


p-'e-"* 


e 


-Tt 


P 


i 


e 


-rt 


diag [ e”^'* 


-T.-k 



] 


e 


-T,t 



( 2 . 4 - 2 ) 


vhere 




\ 

' ^ With these we are ready to derive the desired fundcmental equation 

for i , the state space displacement. The form far this fundamental 
equation is defined in the following: 

^ Definition 2.3: h spatial form for the state space displacement 

is defined as the form 

T 

;• ^ » C I ®< (^3 

in which C is a set of X ( X s n ) tLme invariant nxl colunn vectors 
c, , in the real n dimensional space, and «< (t) is a set of X 

S 

time dependent real scalar elements «^,(t) ,«(j(t) ,...,«<j^(t) , so that one 
element of *^(t) is associated with every colunn vector of C. In matrix 
notation, if 

C = [c, c^3 ; -<“^te)» 




tiien 


■ 

c. . cj / 


C-h) 




L 


oLt 


' f 

- c.,/»<,C't)<it '♦• I o<'^C't-)a.t- -V ^ 


2.4.1 Reduction of /-Equation to the Spatial Ebcm 


The / -equation from (2.3-6) using (2.3-9) is 

rf./' e"'% uwAt 

0 

Premultiplying the integrand on the right side by a unit matrix E » 
PP“’ , where P is the transformation matrix of A for the Jordan canc.iical 
form (see (2.4-1)) , 

T 

I P P' ‘ B k C-t) A-t 
0 

Substituting for P”*o"^* from (2.4-2) 

T 

s I P P*'b U.W <1* (2.4-5) 

0 

It is convenient to partition the nxn matrices P, P'**' r column wise 
arid rowwise, respectively. Let , k»l,2,3....,r denote a typical 



c 


partition of P consisting of 0^. columns of P, and let denote a 


typical partition of P 


consisting of rows of P' 


Also, 


it is easier to identify the column vectors p. (or the row vectors p. ) 

in relation to the partition to which they belong as well as to die 

global matrix P (or P"‘ ) if we adopt a modified subscript notation. So 

1st p ^ ^ p ^ 2^ 3/ • • • • # r f / 2f 3 f • • • • ^ ^ d6not6 ths pj ^ 

L L 

A T 

0 . vectors Wiere 

' j 


L s f 1 , X » I, 

< 


(2.4-'^) 


The index k refers to the partition P^ , P^ and i refers to the 
location in a partition (e.g., k=2,i*4 refers to the 4th coliinn in 
and (L'+4)th column in P. The same is true for row vectors of P"*'). 

k-> , K K .1 

AJ.90, note that L +*/^«L . The L th row vector of P , k=l,2,...,r, 
is of great significance in our analysis and we will refer to it often. 


So it is expedient to define the following: 


K ± 


Definition 2.4: The colunn vectors {*„ , L = ^ / k»l,2,...,r, 

* T 

which are the transoose of the row vectors p of 
P“* are defined as foundation vectors. 

We give this name to these vectors because these are the (left) 
eigenvectors on which the chain of generalized eigenvectors is built 


corresponding to each Jbrdan block, 
with these notations 

^ ~ Ip k_i » ^ K-» > ^ ic-t 

l\, L L -f 









(2.4-7) 


A 

Pw = 


A » 

K-i 

A T 
L -+2. 


aT 


P. [p, p. Pj P,] 


.1 

P = 


r 


A 

P 


L'-J 


(2.4-8) 


Therefore, in partitioned form the integrand of (2.4-5) using (2.4-2) is 


Pe'^^p-'siL -[^P- 


-Pj 




_Tb 

e * 


O 


-Tt 

e ^ 


r ^ 


A 

P. 


A i 


B U. 


(2.4-9) 


-Tt -I T'- - ‘X.l 

Pc P 8»A s[ p, P^. - - P^J 




P. Bi. 


_rt 


-Tt 


p eu 


(2.4-lC) 


In the above tne mathematical identities , P^^ , could in general 
be complex due to complex eigenvalues of the system. 3ut our system 
(x,u,A,B) is real and so is the displacement S . La order to be able to 
attach real meaning to our results the recovery region must be 
constructed in real space. 

2.4.2 Reduction of the ^ -equation to Real Domain 


Because our system is real ve know that the eigenvalues are either 
real or occur in complex conjugate pairs. For real eigenvalues we will 
choose real eigenvectors and for pairs of complex conjugate eigenvalues 
vje will choose eigaivectors in conjugate pairs; i.e.^ if is an 
eigenvector of A corresponding to an eigenvalue , then is an 

eigaivector corresponding to the eigenvalue because 

k y K (2.4-U) 


Similarly, for left eigenvectors the conjugate property is true, i.e.. 



AT 





o-r 




(2.4-12) 


Ebr real eigenvalues the left eigenvectors are also real. The conjugate 
property by extension is true for all generalized eigenvectors (left or 


"^Throughout this text a bar will be used to denote the 
complex conjugate of a mathematical object. 
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right) . 

Ln (2.4-10) the right side is a predict of two matrices each in 
partitioned form according to the order of . The location iiet cs of 
any partition in a matrix does not matter so long as the partitions In 

the tvo matrices correspond to each other i ^ and 

-T t * 

So we can assime that in (2.4-10) the partitions 

-Xfe * 

e, Pk®**’ if' the second matrix are arranged such that those 
corresponding to real eigenvalues (real ) are on top, followed by 
pairs of partitions, each pair corresponding to a pair of complex 
conjugate eigenvalues (complex ) , ftnd the partitions \ in the 
first matrix correspond to tlwse in the second matrix . 

we will let C denote the class of all real quantities, and C ^ 
denote the class of all complex quantities. For class C the partitions 
^ feal, and for class C they are complex. 

It is now convenient to write (2.4-10) as follows: 

where 

!*• ^ ScV Pa.® ic^ ^ 

is B. se^ of 

Jtow, we can formulate the reduction of S to the spatial form in real 


'I 

« 


L‘1 


(2.4-13) 


domain by means of the follo^in-g thsoran. 

Iheorem 2.1: For any linear time invariant system (2.2-2) the 

displacement / in the state space can always be reduced to the spatial 
form defined in Def.{2.3) in which C, cx{t) are real valued. 

Proof: Consider the equation (2,4-13). But for the sets 
this spatial form of S would have been in real space. Mow, consider 
the expansion of (2.4-13) in terms of real and complex classes, i.e., 

P''e« . (Tj + C* cT 

(2.4-14) 

A ^ 

It is therefore necessary to express ^ ^ in real terms to obtain the 
real spatial &rm for i . Eicpand partiton wise and group the 

resulting partition nnatrix products in pairs according to the complex 
conjugate eigenvalues; i.e., consider pairs like 

-T.t A — 

p. S '* P. 6u -h P. S P. Bw. (2.4-15) 

J J J^l Vi.** ij/ 


(;»hich correspond to the pairs of eigenvalues ■'bte that the 

order - of T. and of X are equal ( %). - V. ) because of the 
conjugate property of the eigenvalue \ and eigenvector js>^. , etc.. 
Furthermore, from this conjugate property (see (2.4-11), (2.4-12)) we 
have the following; 



(2.4-15) 


and hence 





pj., 



•t 


A 


P. S^K 






PjS'^) 





(2.4-17) 


Thus, the second tem in (2.4-15) is the complex conjugate of the first 
term; therefore, (2.4-15) becomes 




Pj S w. + 


^ -T. t A 

P. . , e p. 

^ jfi 



(2.4-13) 

viiere Re[.] denotes the real part of a quantity. We now proceed by 
separating tlie real and imaginary parts as shovn below. Let 




- 1. 


+ u 


Pi 




P" 


Pj 


A 0 

P: + 



(2.4-19) 


>here the superscripts R,I stand for tie real and imaginary parts, 
respe'“tively, of the quantity. Substituting for Aj from (2.4-19) in 
(2.4-2) 


->;t _N-t - C- ^ j tO; )t „ Nj-b 

e *e Q. s e, e” 


iit -ciO.t .N;t 

= e e ^ 

Mjt -N-t \ ] 


(2.4-20) 


Using (2.4-19) and (2.4-20) on the right side of (2.4-19), 






i R*[ P; e ^ p. [ ( 

( p/*i p/)s'*} 

= ^ I Pj^Su s«’»vujt^ 

+ ^ e” ^'* ^ ( B SiV»vjj t — Fj^ 8 n. <=^ )J 

"*■ '• ^ i m a ^ i n«c r j **■ t ^ ^ 

+ p/sv.S,nu:;t- P,^B.c^O,t) 


S«_lt -N;t ^ p/Bt,».P,*J;t) 


(2.4-21) 


This can be expressed in matrix form in terms of partitions, and hence 
(2.4-15) can be written as 





A 


P. Bu. 

J 


1 



PjBu. 


= rj 


5 ^■’ ^ ( ^* *8 uCrJ u)jt 

•t- R^Bw. S.niOjt^ 
:^e'V'"^\p/Bus;nui.t 
— p.^eu c*fcjjt3 


(2.4-22) 


Thus, we see that the complex partitions 


on Lhe left side are simply 



replaced by similar real partitions on the right side. 

Proceeding in a similar manner (2.4-15) can ba expressed in real 
form, as on the right side of (2.4-21) , for all pairs of complex 
conjugate eigenvalues. Now, we can rebuild the matrices In partition 
form to replace and ^ in (2.4-13). Denoting these replacements by 

mm 

for respectively, we can write (2.4-13) as 

P P 6 U. • ^ ^ 

A A 

(2,4-23) 




A 


viiere ^ ^ are sets of old partitions which were real to start with 

A A 

and d'x , ij- are sets of modified real partitions replacing the sets of 

^ A 

complex partitions ^ . Now, the right side in (2.4-22) is 
completely real, and therefore using this in (2.4-5) the state space 
displacement, expressed in real domain, becomes 



(2.4-24) 


where 



•[p, K 


- 7,t * ^ 

€ ' P, BU. 


- ^ 
e ^ 


7^ ^ 

e ‘ Pi®'^ 


, •»* 


4 =[p.' p." p/ r- • 


P. P. 




^ P|*Ce4U),t + P|^3'.ncj,‘t ) Bu^ 

•^ £. ' e ' (P,^si»i«0|t — p|^c«u),'t J bi*. 

^ ( p/c-*o^V + 

( p/s.„^^t - P,^omv,H: ) 8.^ 

I ' 




I 

4 s c y «*< ct) d't 


in which we de£in« 


c s [t, s • • *'0 * C^i 


o< = 


C’t') 

•<^C%) 





(2.4-27) 


Thus, the displacemait / can alvays be expressed in spatial form 
(Def. (2.3)) in rral ^ce R" for the system (2.2-2) as shown in 
(2.4-23). At this point it is helpful to summarize the important steps 
involved in obtaining idiis expression. 


2.4.3 Suiimary 

The procedure to obtain the spatial form (Def, 2.3) in real space 
for the state space displacement <T for a linear time invariant 
system (2.2-2) is outlined in the following steps: 

1. The equation to work with is 

<T s / B uCt-) «it- 

0 

wher" A, B, u € x » Ax + 8u. 

2. Obtain the eigenvalues of A and choose real eigenvectors (and 
generalized eigenvectors) for real eigenvalues and complex 
conjuoate eigenvectors (an1 generalized eigenvectors) for 
complex conjugate eigenvalues, and construct the transformation 
matrix P so that P”*AP • J. Also, obtain P~* . The equation in 


step I can now be written as 



I 



» 


t 


I 


< 




P p”’q 


3. Partition P column wise and row wise so that each partition 

has columns and rows, respectively, and corresponds to the 

-Tt 

diagonal block in e .In partitioned form; 







P,et(.W 



<U 


» 

% 
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If all eigenvalues are real then this is the spatial form for 
S in bT given in (2.4-25), (2,4-25). On the other hand, 
if there are comolex eisohvalues present then proceed 
to step 4. 

4. Consider the pair of partitions in either matrix in the 
equation which corresponds to a pair of complex conjugate 
eigenvalues. Replace each one of this pair by a partition from 

A A 

the same location in or depending on the matrix 

considered. Repeat for the other matrix . htow, £ is reducad 
to the sioatial form in r” given in (2.4-25), (2.4-25). 



2.5 Test for Complete Controllability of the System 


In article 2.4 we were concerned with obtaining the spatial form of 
a in real domain R . In this article we will examine the spatial 
form of ^ anri rejate It to the control lahl 1 1 ty of the system 
( 2 . 2 - 2 ). 

2.5.1 Characteristics of the Spatial Perm from the Point of 
View of Controllability 

Let us consider the spatial form of ^ as defined in Def,(2.3). 
The time invariant columns Cj of C, j=l,2,...,£, span a space vhose 
dimension is less than or equal to t vhere JLSn. . It JL<yi then it 
i= clear that the c. vectors cannot span an n dimensional space. So 
any choice of the control vector u(t) for any arbitrary but finite time 
T can influence only those components of the initial state (or 
disturbed state) vector which lie in the subspace spaaned by the Cj 
columns. As such wa have no control over those components of the state 
vector vhich lie outside this subspace of the n dimensional space. 
Hence, an initial state (or disturbed state) which is such that it has 
one or more of these uncontrollable components cannot be brought to the 
origin in finite time. Therefore, the system is not completely 
controllable. This leads us to conclude that to avoid the certainty of 
uncontrollability X not be less than n, or X • n and 

Cj , j»l,2,...,n are all linearly independent (abbreviated L. I.). **5w, 

by demanding that C be a full set of basis vectors for an n dimensional 
^ace we have created access to all n directions in the space. Is the 



system now completely controllable? Since an initial state can have any 
arbitrary component in each of the n directions it is apparent that the 
components in the n directions must be affected independently over a 
finite time T if the svstem were to be completely controllable. The 

fT 

integral j influences the component in the Cj direction, 

j*l,2, ...,n. Hence, the scalar functions must be linearly 

independent over the period T for complete controllability of the 
system. Otherwise the system is uncontrollable. Thus we can conclude: 

1. It is necessary that C be nonsingular so that the system may be 
controllable. 

2. It is necessary that be linearly independent over the 

period T so that the system may be completely controllable. 

3. The t\o conditions (1) and (2) together are sufficient for 
complete controllability of the system. 

itote that t.here is no single sufficient condition here. That is 
because we are discussing a hypothetical case in \«hich C may be 
singular. In the next section we will show that the matrix C in our 
spatial form is always nonsingular thereby making the second condition 
sufficient for complete controllability. 

will illustrate the second condition for complete 
controllability by a couple 'of examples. 

a. Consider a system with scalar control u(t) with all eigenvalues 
real and n independent eigenvectors (no generalized vectors) , so that 
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'<^.ere |»; ,'^J , 3tc. have the usual sieaning 
and b*8. itow, if there are two identical eigenvalues, \ / and 

the rest distinct, then 
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Thus, two dir ctiwns c. , c. collapse into one new direction d and 

j j-<i ^ 

the spatial form is otherwise unaffected. This indicates loss of a 
basis vector and therefore leads to mcontrollability. Here, the fact 
that made ®<j (4') and linearly dependent (see 

Gramian test later) . 

b. In example (a) consider an mxl multiple control u(t). Tnen 

-r T 

C; J o(. CWAt + j At 

i 0 



The bracket consists of m vectors which are linear combinations of Cj 
and c. . If there are less than two L. I. vectors in this m-set then 



it indicates loss of a basis vector, i.e., tbe systsra is uncontrollable. 
Othen/ise it is controllable. The independence of »<j Ct; is rela*-^»d to 
the number of L.I, vectors in the rt>-sat. 

Another observation here is, if there are /. repeated eigenvalues 
then we could combine X vectors of C in the same manner as in 
examples (a) and (b) . In this case, in example (o; the result is always 
a single vector c (indicating collapse of i.— I vectors) with the 
conclusion that it is an uncontrollable system. In example (b) as long 
as .4 * tiiere is a possibility that the m-set will not degenerate 

below jl L.I. vectors. If ^ > m tdien jf-m vectors have already 
collapsed with the certainty of uncontrollability. 

In this discussion have highlighted the characteristics of the 
spatial form / from the controllability point of view. In the next 
couple of sections we will derive rigorously the necessary and 
suff)cijnt conditions for complete controllability for linear time 
invariant systems {2.2-2}. 

2.5.2 Nonsingularity of the C Matrix 

In section 2.S.1 '.ve examined the spatial form of ^ in its most 
general form (Oef.2.3) for characteristics of controllability. We 
concluded that tvo necessary conditions have to be satisfied 
simultaneously to guarantee complete controllability of the system. In 
the following theorem it is shown that for the spatial form of ve 

derived in (2.4-23) (or (2.4-25)) the C matrix is always nonsingular. 
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Th-3orem 2.2: In Lhs spatial form of t.‘ie displacament in re'-l 

space, derived in theorem 2,1 the C matiix is alvays nonsingular. 

.-'roof: If all the eigenvalues of the system are real then the C 

matrix is the transformation matrix P and hence is always nonsinrjular. 
But, vmen toere are complex eigenvalues present the C matrix is sligtttly 
modified, i.e., a pair of complex conjugate eigenvectors (or generalized 
eigenvectors) are replaced by the real part and the imaginary part of 
the vector in the pair which corresponds to = -T; + i ^1 we 
need here is to prove that all the n vectors of C are linearly 
independent (their positions in C do not matter) . 

Suppose that there are q pairs of complex conjugate eigenvalues 
(n-2q real eigenvalues) . Let us rearrange the colunns of P so that 




(2.5-1) 


Others P| (nx(n-2q)) is a set of all real eigenvectors (and generalized 
eigenvectors) and P^ (nxq) , P, (rxq) are sets such that for every vector 
in Pj^ its conjugate is in P^ ^conjiqate property of eigenvectors) . 
Thus, Pj* P, • "Sie inverse of can be obtained by corresponding 

rearrangement of the 'ows of P”‘. The matrix (rearranged form 

of C ) is obtained simply by replacing P and in by 

P and P , respectively, where P«P+:(^ . We now have the 

lu X X A ^ 

following matrices: 


P..- “ [h P. Pxl 

• [p, a' p."] , p- 


(2.5-2) 



e 


3ut P_, * 2 (nxn) , the mit niatrix , and hen-re substituting for 

P from (2.5-2) 
new 

= [e, (2.5-3) 


v^ere E, , E^, are unit matrices of order (n-2q) , q, q, 

respectively. Equating t±ie left side and the right side partition wise 
in the above equation, we have the following result: 


p..’’ p, * 

wew ! 


. P"‘ P 

I r«w X 


P'* P 


Now, substituting for C from (2.5-2), 


p„:J c... = p,:’ [p, p/ p.^] = p. p".- P-* 


(2.5-4) 


(2.5-5) 


-in -IX 

We can obtain P P , P P as follo\«: 
n&w t new 2. 

p.:! pv + p.;L px - p~- ( ' p-«- 

Cf" ^ iiCA-p: A) 


and using (2.5-4) 



p"‘ 


(2.5-6) 


Similarly, 


p-> pi 


- - ^ p'* p - p”’ p ) 


or, using (2.S-4) 


P’‘ P^ 


X i 


Ex 


(2.5-7) 


Substituting (2.5-4), (2.5-5) and (2.5-7) in (2.5-5), 


P"' C 

'r\9u n«w 


E. 

0 

0 
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-ie. 

0 

1 < 
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(2.5-8) 


The rank of the matrix on the right side in (2.5-8) above is n, and 
hence P^*'^ is nonsingular. Therefore, is nonsingular. 

This proves the theorem. 

Thus, ttie first t;ecessary condition that C must be nonsingular to 
avoid certainty of uncontrollability (see secton 2.5.1) is always 
satisfied for the spatial form of i in real space I?". This leads us 
to conclude that the second necessary condition, i.e., the requirement 
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of the linear independence of the scalar time dependent functions (t) 
of o<(t) a/er the period T is also a sufficient condition for complete 
controllability of the system. Hence, it is only necessary to test for 
the linear Independence of the ^j(t). But as it stands now it is not 
easy to conduct this test, and it is desirable to derive some conditions 
which will be necessary and sufficient to guarantee the independence of 
the «<j(t) . The next section is devoted to this objective. 

2.5.3 Mecessary and Sufficient Conditions for Complete 
Controllability for Linear Time Invariant Systems 


In this section we will derive a set of necessary and sufficient 
conditions which will guarantee the independence of over the 

period T, thereby guaranteeing complete controllability of a linear time 
invariant system. 

Fbr the spatial form of 4 in real space (see (2.4-25)) , from 
(2.4-25), 






where is the set of vector partitions 

A 

f snd ^ is a similar set of pairs of vector partitions 
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for all in wiich + (see (2.4-24)). Each 

vector partition consists of elenents (scalar ti-ne dependent 

functions) vhere is the order of the Jordan block 3^ . within each 
vector partition all the alenents are linearly independent over the 
period T provided there is no zero element present. (These elements are 
menbers of a Jordan chain; each element differs from the other in the 
highest power of t, because of the polynomial in t in the exponential 

A . 

matrix e” (see (2.4-3)).) furthermore, between any two vector 

partitions in ^ if the exponential factor of one differs 

from that of the other then these two sets are linearly independent of 
each other over the period T. Similarly, between any tvo vector 

A 

partitions in if either fj or W; of one differ from that of the 

f , ► 

other, then due to the exponential factor or the factors 

cos (Vjt., sin (J;t , these two sets are L.I. over the time T. Airo 


any vector partition of is L.I. of any vector partition of 

over the period T due to the presence of the circular functions of 

A 

in ^ . Thus, the test can be narrowed down to checking the vector 

A 

partitions in and separately, and 1) those partitions in 

A 

for v«hich are identical, and 2) those partitions in ^ 

for which both 'fj , tjj are identical. A further simplification is 
obtained by using the Jordan chain property for the elements of a vector 
partition, i.e., if the ))• th elanent of one vector partition is 
independent of the >)• th element of another vector partition then the 
tw 3 sets are independent of each other. Because the elements of a 
vector partition are a cha’n built over this th element. (An 
analogous reference is the i . Jation vector of Def.2.4. Actually the 





•J; til element is male up of the foundation >/ector.) For convenience 
will call these J^th alaments in each vector partition as "foundation” 
elenents in «»c(t). Ihus, finally the test for linear independence need 
be conducted only for the foundation elements of: 

1) the partitions in for which are idaitical, for all 

; and 

A 

2) the partitions in (f-g for vhich both , W; are identical, 
for all \ € (T*. 

A ^ 

Let (in or ^ ) denote distinct eigenvalues of the system, 

and let be the nunber of Jordan blocks (each associated 

with a \ , k»l,2,.,..,r/ r < n) , for which >u are identical and 

A 

equal to . Let 5j be sets of foundation eleneits of e<(t): 1) 

^ h A 7T 

Corresponding to >j for € C , and 2) Corresponding to 

A ^ 

for >j * • Since each Jordan block gives rise to one foundation 

A 

element each Sj contains /'j elements if is real, and 2/^j 

elements if is comple;;. For the complex case, Sj is a set of all 
elements for which both and U)^ ( identical. Let 

A 

Vj denote: 1) sets of foundation vectors (see Def.2.4) |» ^ 

A 

corresponding to all the Jordan blocks associated with X- for 

A ^ ^ I 

>j ^ , or 2) Sets of the real and imaginary parts b , b ^ of 

A ^ 

the foundation vectors corresponding to all the Jordan blocks 

associated with either one of a complex pair of ,X\^C . 

The sets y^- contain vectors (colunns) for X\ ^ C and 2/<; 

A 

vectors for € ^ . Or, 

-y,. - ([^. ■ S/.;] ^ 

[■J, ^ 
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m m 
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««<, >), are usad in Place Of , a«. to simplify 

notation. Now, usl,^ (2.4-2), (2.4-3) and (2.4-7) in (2.4-24) for &a 

fomdation elemanea of ,^(t) (in r, , 4 > »» sets S; c.,n be written 
as follows: 

r T X "I 

Sj - L® 

for € (T 

^ Ccnwjt + SinCi^t j 

A 

5 e'‘ Sm^-t - ( 

A 

«i € { C 'i^) e«*.C4*au)yt ^ j^6«w SiKitOj-fc } 

a e^‘ '' [ tfj'')^euai>>c^t - ( 3*)’8w.crj4jjt 
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for c ^ 


(2.5-10) 


Note that toe elements of S- are scalars. For % e (T ^ Sj can be 
expressed in toe matrix form 


nt' r . 


(2. 5-11) 


-^jt TJT 
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where yTau » \7^y. , etc,, and Y, is as -^iven in (2,5-9) for 

A A ^ 

C • lb express S. for e. C in a similar form, let 
U.c*ju>jt and usi'wwjt , Since each slanent of Sj is a sum of 
two components, Sj can be written as a sum of t\o sets (matrices) 
(with a slight rearrangement^for convenient algd^raic signs) as shown 
below: 


S; 
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(2,5-12) 


Rirther, if we define 


k s C > 


• /4^. Mocks. . k'] 
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-I 


then (2.5-12) can be expressed as 


S3 . <1 


BVjk 
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0 




(2.5-13) 


(2.5-14) 


tlhe order of elements in a set does not affect their 
relationship with respect to independence or dependence. 



^ m 

whare Yj is as defined in (2.5-9) for \ . 

'liie time dependent scalar functions of Sj in (2.5-10) are said to 
be linearly independent over the time inter/al [0,T] if the Gramian 
matrix associated with Sj is nonsingular, ^lore details on Gramian 
matrices can be found in [ 3 .]. The Gramians G. associated with our 
sets Sj are given by 

Gj s I Sj S; (2,5-15) 

0 


Substituting for Sj from (2.5-11) and (2.5-14) and noting that B,Yj,K 
are time invariant, (2.5-15) takes the form 


G-j = f y.) , e ^ 


vhere 


[( ^ is 1; ^ 


M = y u.uT^'^ 

0 

M = 






(2.5-lfi) 










(2.S-17) 


The matrix M is of order m for \ g ^ , ai*-' -f order 2m tor '>j e 
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Ths Gra.Tiian Zj is of order /<j for ^ j 6 ^ and of order Z/*j for 

Hence, for linear independence of all in , 

i*l,2,....n, all the Graroians Gy must be nonsingular, i.e., rank Gj 

A 

must be equal to /<y for >y e ^ , and rank Gj must be equal to 
2/ij for c From (2.5-15) and (2.5-17) we can write 


'Y'txnk (jJ’j ^ nr«\.*v k. M 


and 


‘Vo.'O k (Sly ^ vetnic ^ ^ 


'Va.'nk 'ro.nk M 

-reink ^ *Y^«k. 


e <r* 


b\ic 
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(2.5-13) 


Now, we are ready to formulate the necessary and sufficient conditions 
for complete controllability of the system. We prove this in the 

following theorem. 

Theorem 2.3: In a linear time invariant system 


X - Ak + 8u. 


(2.5-19) 

where xeR*'and ttc R*", let J be the Jordan canonical form of A and P be 
the transformation matrix such that « J. Let be the Jordan 

blocks of order . Associated with each is an eigenvalue 

A 

> , k*l,2,...,r, r^ n. Let , j«l,2,...,x denote distinct 

K J 
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eigenvalues so that jtsY. Let /<j be the nunber of Jordan blocks 
for which are identical and equal to >j . Let yj be; 1) 

A 

Sets of foundation vectors [see Def. (2.4)1 j? ^ corresponding to the 

A ^ 

Jordan blocks associated with 6 ^ » or 2) Sets of the 

A ^ A Ay 

real and imaginary parts of ^ (i.e.# and ) corresponding 

^ *3r ^ ^ 

to the Jordan blocks associated with or • 

ihen: 

1. Ihe necessary conditions to be satisfied .9o that the system may 
be completely controllable is 

A ^ ^ 

Wi ^ j (2«5“20) 


where m is the maximun available nunber of independently 
variable scalar controls u^(t) belonging to u(t) (i.e., 
dimension of u(t)) . 

A * ^ 

Note: For only one of a pair of complex conjugate A; , 

A 

)ij , need be considered. 

2. The necessary and sufficient condition to be satisfied so that 
the system will be completely controllable is 

Y<ov»ik & *■ YA.uk Yj j u 

s YAvik Yj » 2/<j V 

(2.5-21) 


Ic 


b'^.K 


for ^ C 


Note; For *>j € ^ rank Vj « and 

rank Yj ■ . ^ain for ^ only one of a pair of complex 


c*' HVjf. 



conjugate >j , need be considered. 

Proof; 1. Necessary condition - To prove the necessary condition 
(2.5*20) we will prove Lhat if this is not satisfied then it is a 
sufficient condition for uncontrollability of t^ie system. Hence, assume 
^ • f *'^ a <>•*»« ^ j « I , •..•,1 (2.5-22) 

Consider the two cases: 
a) % € (: 

The matrix B is nxm and Yj is nxyH. so that ^Yj is myHj, and 

T f 

rank B Y. $ m. But from (2.5-18) rank G- $ rank BY.. Hence, if 
ra < yM; , 

<v«\w K Gj ^ "r»L*i k b'Yj < ^ ^ /^J 


Ya««k Gj < 


(2.5-23) 


Thus, the Gramian is singular for this which indicates 

uncontrollability of the system (due to the existence of a dependent set 
S; in «<(-b;). 
b) 

The matrix Yj is of dimension nx2ykj so that sTyj is ro<^j. 

T 

From (2.5-13) K is a square matrix of order 2/ij, so B Y> I< is also of 
dimension mx2y*tj. If call 


A - 


J 


(2.5-24) 


then dimension 2nK2yMj. And rank "jf. $ 2m. But from (2.5-18) 



(2.5-25) 


5 f 


rank CSj < rank 7^. . Hence, if n <^j , i.a., 2m < 2^j, 
-y’A.nk G»; -r^f\k /y ^ 2»*v <■ 


or 




A ^ 

1. e.^ the Gramian is singular for this and this is sufficient 

condition for uncontrollability of the system. 

Ihis shovB that the condition m </^j for some '>j (real or 
complex), j»l,2, , is sufficient for uncontrollability of the 
system. Hence, we conclude that the condition 

A ^ / 

Irvt ^ V , <»'•',*•» J ^ 

is necessary so that the system may be completely controllable. This 
proves the first “a'-t of the theorem. 

2. necessary and Sufficient Condition 

a) To prove the necessary part consider the ranks of sV- 

•I • 

M^ose dimensions are mx^^. and 2m2/ij , re«?»;;ectively. Tf .ank 
B Vj A then rank B yj <^j . Similarly, if rank fj , then 

rank 7^. < 2yHj . Hence, from (2.5-13) , 
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k Gi j 

$ Y«^nk 

•'A 

> : € <r 
) ^ 

nr«OA k 

^ Aravtk 

/j < =A 

A ^ ^ 

) '•j 


Or, 



€ ^ 


k, (jij < j ^ 
•voLnlt j 



€ ^ 


(2.5-26) 


and from tho first part we know this leads to uncontrollability of the 

A 

system if (2.5-26) is true for any . 

Hence, the necessary conditions for controllability are: 


orank 6*Yj » y^j 


X; € r 



€ c:* 


(2.5-27) 


which imply 

\m •.Al 6 ^ (2.5-23) 

^ ^ A 

Also, rank y. is A for >j ^ 6 , and 2/<j for j € ^ , because y. 
are sets of left eigenvectors or their derivatives (real and imaginary 
parts; for this case Lhe proof of L.I. is similar to theorem (2.2)). 
Iherefore, (2.S*27) can also be expressed as: 


'YCKAk ® 7j * • /^j , ^ 

-fav^k 7j = Yj - V'J i C (2.5-29) 


b) To prove the sufficient condition consider the quadratic form ir 
where: 1) V e R^’'* for ^ , or 2) V € R*’'^''for Xj C 

and Gij is as given in (2.5-15). We can ther write 
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v'^^jv . f( e‘^Y;vJ H Ce Yi'') , 

(y.v)'^M (/jV) (2.W!)) 

vhere M is as defined in (2.5-17) and /. as in (2.5-24). Also, let 
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[ '' - 
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e C 
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1 



(2.5-31) 

v^ere z 6 R for 

A 

^ , and 

z frr \ € r*. 

Then 

V • 2'*‘m « 


(2.5-32) 


It is now useful to recall the concept of controllability itself. 
Mien we ask if a system is completely controllable we imply all the 
freedom, if necessary, in the choice of scalar controls «*jCVJfor a given 

nunber of them (i»l,2, ,m) , i.e., the scalar controls 

u^(t), i«l,2,o..,m are independent). y variable. Given this freedom we 
ask if the system can ever be controlled. Hence, we are at liberty to 
use all the freedom available for a given set of scalar controls. 


Now, 
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, J; * C 

-T 


U (1.5-33) 
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Me see that M is a Gramian matrix for both cases associated with a set 
£(t) of scalars, i.e., 

I ^ J ^ > -I , (2.5-34) 

^ A 

^ /T 

Since H-COcan be varied independently for all ^ over the period 

T, rhe scalars of £(t) are also independently variable over T. 
(exponential factors do not affect the independence of , and CMtvyt 
and SinJjjc being orthogonal functions preserve the independence of 
. ) . As such the Gramians associated with sets f(t) can be made 
nonsingu.' r. Hence, assixne N can be made positive definite (N > 0 }. 

Or, looking at another way. 


^ 

T * 

® T (2.5-35) 

» / W Ct) WU) 


Inhere z is a constant vector, and w(t) « f(t)z is a scalar. If 

z"^Mz ■ 0 for all time t « [0,T] then from (2.5-35) w(t) « 0 is the 
only oossible solution, or f(t)z*0. But f(t) Is a set of 

scalars v«/hlch can be varied I ndeoenden: 1 y over the oeriod T. 
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So for this equation to be true for all tine t € z 

must be zero. In other words. If for an Indeoendent 

set f(t), 

/ wctl^d.b 
0 

can alvays be made positive for a finite time f i*e., >o , 

B /o • Thus, a positive definite M can be assuned as vtrlthin our 
capability for controllability purpose. Hence, assuming M >o , we 
know that if 2*^M2eo r then o is the only solution possible. Or, 
from (2.5-31), 

A 

g « v « o , X- ^ 

2 « V - O , C* (2.5-36) 

* T“ 

Since V is x 1 -f*T if r^.ik then v«o is the 

inique solution to 8 v*o . Similarly, v is l/'j- * I for C ^ 
and if tds)k "f, then v«o is the only solution to V«o . bi'*. 

this implies v.o is the only solution to v «o This 

will require to be nonsingular. Hence, the conditions (2.5-29) for 
all X, € C,C Are sufficient to guarantee nonsingularity of all . 
Combining with part (a) the necessary and sufficient conditions for 
''omplete controllability are, therefore, given by (2.5-29). This 
completes the proof. 

In the above theorem we have derived the necessary and sufficient 
conditions for complete controllability of the system. Although the 
necessary condition (1) is implicit in the necessary and sufficient 




co'^'iiltlons in (2)^ the for'^er can serve as a quick and simole 
check to see if the system is definitely uncontrol lahlc. 

Even though the above conditions were arrived at by a rather 
lengthy process their application in testing for the complete 
controllability of the system is not very complicated, ihe first 
necessary condition is very simple (involving only costing of the 
Jordan blocks) once the -Jordan canonical fern is obtained. 

In the second rondition the rank test involves only much 
smaller matrices (mxya^,mx2^j* ), and n and /<• are {generally 
much smaller than the dimension n of the system. This v/ill 
reduce numerical errors In the comouter methoHs. Mnrenvf^r, 
the ^irst test orovides a minimum number ''or the scalar 
control? u;^ that is absolutely necessary so that the syste~. 
may be comoletely controllable. 

In this section we have been concerned with the independence of the 
scalar functions of »<(ir) in the equation for ^ in spatial form in 
real space. We derived the necessary and sufficient conditions %^ich 
would guarantee this independence. In the following ve sunroarize the 

practiced steps involved to conduct the test for complete 
controllability of the system. 

1. Obtain the Jordan canonical form J of the matrix A. Count the 
nunber of Jordan blocks associated with each distinct 
eigenvalue. If this number exceeds the dimension of the 
control vector u(t) for an/ dist..oct eigenvalue then the system 
is certainly uncontrollable. If it passes this test proceed to 
step 2. 


OS'' 



2. Set up the Yj aatrices from the foundation vectors 
corresponding to the set of Jbrdan blocks associated with each 

A A 

distinct eigenvalue . (Fbr complex Aj the real and 

imaginary parts of the fomdation vectors must'be used.) Test 
T ^ ^ 

the rank of B Y> for all ^ C , and test the rank of 

for all If their ranks are equal to the 

ranks of V- (for C, C* ^ case may be) for all 

distinct then the system is completely controllable. If 

not it is uncontrollable. 

Some implications of these conditions of controllability are worth 
noting. Consider a scalar control systen (m « 1, u « u J . From the 
necessary conditions for controllability, we have 

/Ij S M * 1 
or 

^ ^ (2.5-37) 

This implies that no two Jordan blocks of J S(«s Ad be associated 

A 

with the same eigenvalue (i.e., , k«l,2,....,r must all be 

distinct) . in order that the system may be completely controllable. 
Further results can be stated :n the following theorem. 

Theorem 2.4: If all the Jordan blocks J,, , k>l,2,...,r, of the 

A 

system (2.5-19) are associated with distinct eigenvalues , i.e., j ■ 
k, then the system is completely controllable if none of the conponents 
of the vector a<(t) is zero, i.e., 

for any i, i • l,2,...,n (2.5-38) 





Proof: From theorem 2.3 


*1 » 3 * ^ * l»2»..../r» r ^ n (2. 5“ 39) 

Mhere is the nunber of corresponding to the distinct eigenvalue 

A 

. The necessary and sufficient conditions from theorem 2.3/ give 

YS^nlc ^y. m. T«ntc y. « 1 j ^ 

ro^k /. * V; * ^ ^ ^ 

^ (2.5-40) 

\iiere 


and 

Yi 


(2.5-41) 



in which y is a single foundation vector corresponding to a Jordan block 
J. ; y^ / y^ are the real and imaginary parts of the foindation vector 
y or 7 (algebraic sign for the imaginary part ignored) corresponding 





to a pair of Jordan blocks Jj , J; associated with a pair of 

A -T 

complex conjugate eigenvalues . Hence/ using (2.5-41) 

in (2.5-40) 


rank ^Y. « rank » 1 for all A: ^ ^ 


(2.5-42) 


*Y«^nk. 


/. . «..C ."/J [» 


6 / 


Y<»nk 


all 




bV 


= :i 

a II ^ j ^ ^ 

eV 

8^a^ 




(2.5-43) 


^ T 

Fbr real , B y is an mxl vector, and hence 


V4nk 8 ^ <1 

s o 


if and only if the vector B y » 0 . But from (2.5-11) , this means 


-it T 

5?^. * € * 8^ Yj * O -frr S*"*e « <T 


(2.5-44) 


'f" ^ 





V 


And, S. is an element of o<(t). Ihus, 


'vok»>k B D 


_ n 

}*"jj < 1 »•**»« 'ij € ^ 


(2.5-45) 


if and only if a component of o<(t) is zero. For complex Aj , the two 


vectors, 2imcl, 




6^^ 


6 ^ 


6^/ 


are L. I. if ary one partition B y*^ or 5 y'*‘ is not zero, i.e., these 
are linearly dependent if and only if ^y^ • 0 (this can be 

proved easily or observed directly) . Or, 


*y^n k ■* 6 U 






^ T X J ^ # 

if and only if By ■ B y « 0 , i.e,, j. ■ 0 for some € <T . 


But from (2.5-14) this implies 




(2,5-46) 


■ * 



(2,5-47) 


And S. is a pair of elements of <K(t). Thus, 

“Train k. {, < Z »e»*»e <T 

i 

if and only if a pair of components of «<(t) is zero. 

* T 

This shows that for all % , real or complex, the ranks of B V* 
or Vj are full if there are no zero elonents in o<(t) , and full rank 
for these matrices is a sufficient condition (theorem 2.3) for complete 
controllability of the system. This proves the theorem. 

In this section we have discussed the spatial form from 
controllability point of view, and proved a few useful results. Now, we 
are ready to proceed with the construction of an a{:proximation for the 
recovery region ^ . 



2.6 Approximation for Recovery Region 


In section 2.1 we discussed the recovery region ^ and some 
approximations for 2 dimensional systems (Pig. 1#2,3). Die 
approximation was in the form of a parallelogram enclosing ^ , and 
thus forming an i^per bound for . Here, we will extend the 
parallelogram concept to the n dimensional space and obtain an n 
dimensional parallelopiped. The n directions needed to span the space 
cone f’*o(r. the colunns Cj of the C matrix in the spatial form for ^ 
(see(2.4-2S}) . Thus, a set of n vectors in these directions will form 
the semiaxes of the n-0 (D for dimaisional) parallelopiped. The 
magnitudes of these vectors are provided by the magnitudes of the cj 
vectors as well as by the components of <»<(t) . Since we seek an upper 
bound of the region for a finite time T with bounded controls, we will 
maximize each time integral component of a<(t) for the period T using 
the upper bound for all the controls u^ , i=l,2, ...,ra, i.e.,|uj *1 . 
Thus, each semiaxis of the n-D parallelopiped represents the maximun 
component a state space displacenent can have in that direction in order 
that it can be brought to zero in time T. We mentioned in section 2.3, 
for 2 dimensional case that at least one point in the recovery region 
fi lies on every side of the parallelogram. Similarly, for the n 
dimensional case at least one point in fi. will lie on every surface of 
the n-0 parallelopiped each of viiose surfaces consists of n>l 
independent vectors as edges which are parallel to n-1 of the n 
semiaxes. The semiaxis that does not ha^ ^ a parallel vector in the 
surface is the one on whose tip this surface rests (semiaxis is from the 



origin to the center of a surface). Of course, an n-D parallelepiped 
has n pairs of parallel n>l dimensional surfaces. 

Let A* re{«esent the approximation, the n-0 parallelopiped, to 
the recovery region which is to be constructed from the equation 
(spatial font) (2.4-'26)} 


Define 



T 

Ct) 4t 


0 



/ C*-) 'it' 
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( 2 . 6 - 1 ) 
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Dy • ■ ■ • ■ ■ ’^"1 

T 

j " ) 4t- ^ t « L0(T] , ) '^ J • 1 ^ 
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( 2 . 6 - 2 ) 


In terms of these notations A ^ can be denoted by a set of n vectors 
, j«l,2,...,n in matrix form, these vectors forming the semiaxes 


of the n-D parallelopiped. Or, 



(2.5-3) 


• C Dy 


U.1, /J 


To obtain an approximation to the degree of controllability P of 
the recovery region ^ we have to get the minimun of the distances from 
the origin to all the surfaces of the n-D parallelepiped. Ihe distance 
from the origin to any surface is measured along the normal to that 
surface fran the origin. Ihen the normal distance from the origin to a 
surfaceof tiie n-D parallelopiped is the component of its tip vector*^ 
along the normal to this surface. In the following theorem it is shown 
how to obtain the normal distances to the surfaces of an n-D 

parallelepiped. 

Tneoren 2.5; Let F * [ a^ a^ a^ ] prescribe an n-D 

parallelopiped in real space whose ssmiaxes are given by the linearly 
independent nxl columns a.j , j»l,2,...,n. Ihen the nonnal distances 

dj to the surfaces of this parallelooiped are ?iven by the 
reciprocals of the masnitudes of the column vectors of 

Proof; Let (F^)“‘ » [ a, ] in which a. , j»l,2, ...,n 

are nxl colunns. Since F"‘f ■ E , the init matrix, we have 


r* **< 

•■I 

[o., . . .. %] . 

m ^ 

1 



0 

V 


» 

% 

• 


J 

1 


0 

I.*- J 
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A tlo vector of a surface is the semiaxis on whose c’ 
the n-1 dimensional surface rests. This surface loes 
not have ed«?es parallel to its tip vector. 



Oy, 
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1 

2(7^ 
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'■•om which 






(2.5-5) 


where is the Kronecker delta. Equation (2.3-5) is the 

restatement of the fact that oL. is orthogonal to all a^ , k 9 * j . In 
an n dimensional parallelopiped a normal to a surface is orthogonal to 
n-l vectors of which the surface is composed of. Or, a normal is 
orthogonal to n-l of the semiaxes. Since a^ is orthogonal to n-l of 
the a^*s, it represents the normal direction to the surface whoso edges 
are parallel to a^^ for all k except k*j, Tho tip vector for this 
surface is a^ , and the normal distance d* from the origin to this 
surface is the component of a. along the normal direction given by 



(2.5-6) 








|; 

t 

and using (2.6-5} 


> 

i 
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(2.6-7) 



vhere \\ a"- H is tha magnitude (Euclidean nonn) of the colunn vector *?. 


• 

j 

of (f'")*’' . 1!iis proves the theorem. 



Itow, from (2.5-3), 
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(2.5-3) 
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Let 




CO" =[j, f, ■ ■■ 




% = Ilijil h , , 

(2.6-y) 
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Vihece g, den^to aiit vectors. Tiien 
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(2.6-10) 



Using theorem 2.5, the normal distances from the origin to the 

sur faces 



of the region are given by 





(2.G-11) 


) 

t 

i 

Tne approximation p* of the degree of controllability 

mm m ^ ' 

from 

A ~ 



the parallelepiped bound on the recovery region ^ is 


(0*^ Wh'a d- ^ j -- 1, V, , n. (2.6-12) 

5 

This approximate decree of controllability can be made arbitrarily tight 
by including additional directions in the n dimensional spare. Suppose 
e is any desired unit vector in one such direction, then maximize 

} (2.6-13) 

b 

over the period T with |u^\= 1, i=l,2,....,m. This quantity represents 
the maximan componer.t in the direction e of a state in the reco«/ery 
region A , and hence to points in ^ lie beyon^ the surface ortlio^onal 
to e and a distance given by (2 5-13). Tias is similar to ’.^i3t we 

discussed in section 2.3 or a two dimensional state space. By 

considering a 3, t of such e lirections e set of distances ( 1 . 5 - 13 ) :.=• 

A 

obtained. Let be the minimum of these dist<anccs, Tlien an improved 
estimate of the degree of controllability is (° = nun( ) , and 

f > (p can be made arbitrarily close to tiie tru; dcjree of 

controllability ^ by picking a suCfi'ient number of directions e. 

The approximate decree of controllability (O* in (2,6-12) will go 
to zero wnen Lhe system bccoaits uncontrollable as inJicatij in the 
following theorem. 


Hieoren 2.6: Suppose J 

is the Jordan 

canonical form of 

A i,. 

the 

system (2.5-19) and 

the Jordan 

blocks 

associated 

with 

the 

eigenvalues , k=l,2,..,,r 

, r 4 n. If 

the 

eigenvalues 


are 


distinct then the approximate degree of controllability is zero if 
and only if the degree of controllability is zero, i.e., based 

on the minimaa normal distance to the surface of the parallelopipad 
will be zero if and only if the systen is uncontrollable. 

Proof: The normal distances to the surfaces of the parallelopiped 

A* , from (2.6-U), are 




i 



s I 


/ h 


wnere /• are defined in (2.G-2). Since (^* « min d- , it will be 

^ s ^ 

zero if and only if some dj » 0 , out dj »■ 0 if and only if • 0 . 
Fro-n (2.0-2) 

T 

jf. c f dr , t e l«.“3 , 


and /• = 0 if and only if * 0 . Thus, p* * 0 if and only if 
some - 0 . From thee 'em 2. A the systen is unconcrollaole if and 
only if ths4& exist one or more zero components in o<(t) , i.e., some 

o<,(t) =0 . Hence, » 0 if and only if the system is 

uicontroliable. 7nis proves the fieorem. 

W\en t\r> ot more Jordai blocks are associc^.J the same 

eigenvalue some modification is necessary so that f** be zero if 

and only if trie syste:?. is uncontrollable. We adopt vhat is known as 
singular value decompositior ♦’o modify the foundation elenents 
corresponding to the set of Jordan blocks associated with one 

eigenv^.lue. We denoted such a set as S; (see (2.5-13) cod will use 



the 5eme terminology for all quantities inless otherwise mentioned. 
Similar to the approach adopted in transforming complex vectors in the 
spatial fom to real vetors, we will group the vectorssflf S that 
contain the aleoents of Sj and tranform this set into a similar set so 
Lhat the same process of replacement in the spatial fom can be adopted. 
The new elements of o((t) will give the property desired. We will 
illustrate the process for real eigenvalues, and for complex eige.T/alues 
it is exactly jinilar. 

Consider a group of vectors (a siznmation) of S associated with 

A 

S- corresponding to Lne eigenvalu? V , i.e., 

^ I ®<v* (2.G-14) 


U)ere tae «(.(t) belong to lin oolonn form it is denoted by • 

Denoting tiw sot of c^ h. C.. (nx^^ matrix? o;ii s^ostifjtinj for 
fro.a (d.S-ll) this equation becones 



<> 


(?.6-l5) 


Since the necessary and sufficient coivdition is related to the rank of 
the matrix b\, (see theorem 2.3) this is to be modified by singular 
value deco.oposition toolviique. Details of this tochniqjo con be foind 



in many taxts on Linear Algebra, •«<* » reference [4], and here ue 

will just use the results for our purpose. By this technique the matrix 
T 

YjB can be written as 


V7b.L;WjZj ^ 

in hbich Lj , Zj are orthogonal matrices of order ^ 
respectively, and W. is a /»; x ti matrix of the form 


(2.6-15) 


V)j • ["A; ol 


o 0 


(2.6-17) 


Tna A- is a diajonal matrix consisting of as many nonzero values 
called singular values as the rank of Wj . Substituting (2.5-lS) in 


(2.5-15) we ;uve 


T * T « ^ 

Si VT® j «■ *»■ * Si L; / H Zi ». 

• • 

ss C .. I*. 




(2.5-18) 


wlisre C^. » L- which is a set of L. I. vectors. Mow, let 


>* X t 

S. - e vr. 4.U. 


(2.6-19) 


: 


denote the elements «(j^(t) , i«l,2,...,^^ . With this (2.5-18) can 


be written as 



( 2 . 6 - 20 ) 


c r- 


at- 


This fonn is ready for replacement in the ^ equation, i.e., replace 
and o(.{t) by o(.(t) . If A. is of order £; (v^ich 

w j * 


Si 

indicates the rank of W. is ) then ?xpandino Vi . Z> we will get 
« 

^ ’• •• -.*<1. . ® . (2.5-21) 

This has Z';- £. zeroes for tlie oi';(t) coefficients. 8ut beinj 

Uie rank of yJo if it is less than tiie system is ancontrollable. 
In other words, if there is a zero elenent present in Sj th^n tie 

system is uicontrollable. On the other hand, if £^. is eqi’il to 

T 

Y. 0 is of full rank wnich satisfies ttje sufficient conditior. for 

j 

cor.plete controllaoil it> «*s far as this eigeiwalue is concerned. 

Hence, the coefficients ft) are independent. Thus, only a zero 

4 

indicates dependency of ®<^(t) . Ttiis process is re:. 'ated for all 

A 

distinct real ei.jenvalues . 

For complex eijetwalues the sane procedure is adopted. Tie 

equivalent of (2.6-15) is 




(2.5-22) 


where an ax2/»j matrix. Ajain is deconposed 

similar to (2.6-16) (note Y^ is a nx2^^ matrix), i.e., 

y/fij = l.vj;€ 


(2.6-23) 





Hance, (2.6>22) becomes 




<k^ 


(? - 24 ) 


r 
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where Lj , a set of 2/^- L.I. vectors, and 


A 

S. 


t W. 2. 


bj 


(2.5-25) 




denotes tne 2 /Mj elenents «<;(t), i=l,2,....,2/5; . Tie ^orm in (2.6-24) 
is ready for replacancnt in tJie ^ eqaation as discossaJ earlier, ftow, 
sup;>oso Wj is of rank i- tlun rank » Xj . And if tj < 2Aj then 

A ^ ^ T 

there will be 2^.-0. zeroes of o(-(t) in S> . Tie syste.n, of course, 
is a;icoiiirollable (rank /j < 2/<^ ). On the other nan), if Xj * 2/lj 
ttie S'.'Ste.a is cor.iplotaly control laolo (rank/^ * 2/1; ) es far > this 
ei'jen'v’alua is concerned. Hence, all 2/i. coefficients o<;(t) are 

iniapunJcT.t . Tms, only the presence of zeroes indiectoi depenJ^ney of 

j(.(t). Tiis process is repeated for all distinct complex eigenvalues 

A TT 

or . 

In tie abo/e a modification is describe) to be aaoptei in the 
spatial form of <f for the colonns in C and cie elenents in »i(t) in the 
case of more than one Jordan blocks associated with one einsnvalue. 
T.iis involves replacing a set of c^ in C by anotser L.I. sot c^ and 

A 

a correspond inj replacencnt of o<.(t) in o<(t) by Of'.(t) which correspond 
to the foundation elenents associated wi h an eigenvaluo. We showed 

A 

tnat the systen is uncontrollable if and only if an elanent o«t) 


I 



goes to zero. Combining this with theorem 2.S and following the same 

line of proof wi conclude that the approximate degree of controllability 

|| 

(• goes to zero if arvJ only if the system is uncontrollable. This can 
be sunmarized in the following theorein. 

'Riaorem 2.7j For a linear time invariant systen (2.5-19) the 
approximate degree of controllability (*** based on the parallelopiped 
region goes to zeio if and only if the systsn is uncontrollable. 

If there are more than one Jordan blocxs associated with an 
eigenvalue apply modification based on singular value decomposition as 
described above to the sets of colunns of C and foundation elements of 
of(t) associated with sets of Jordan blocks corresponding to any distinct 
eigenvalue. The rest of the proof is as given in theorem 2.^. 



2,7 Sumnary 


In this chapter attention wns directed to iencr=»l line^ir time 
invariant systens. 've startea witn tho qjcstion of how effective is a 
distribution of actuators on a very lar^o flexihle spncecraft in 
controilinj its attitude and shape? Tb comprahoni the meaning of 
"effective" sought to develop a concept of the decree of 

controllcibility. Ibis develop lent followed a ratioml nooro»ch first hy 
showing the txisuitability of certain candidates to serve as a definition 
for the degree of controllability, and then resulting in a meaningful 
riefinition wiiici ^vouli account for all tv? pertinent factors sue!', ns 
controllah.lity, total time, control effort, stability, an-* -ontrol 
objective that will have a bearing on the degree of controllability. 
Once t definition was foriulated suitable approximations had to )C 
Jevclo-xid for the recovery region and the degr^^.e of controllab' I ity so 
that Mis definition is easily applicable to real problems and 
nunerically nanajeable, tne jpproxi. cions approaching tha tru* values 
as th.j comiXitational effort is increased. The mathcnatical approach 
adopted, besides leadim us to our desired go*l, sho.v-*d the system 
equotiens in special forn Lheropy enabling us to derive some relatively 
simple tests for complete controllability of the s'ystem. 

Hius, in sliort, the emphasis in this chapter hns been developnont 
of a rational concept to obtain a moaningful definition for tv? jegreo 
of controllability and taking it to the st^ge of usefulness from an 
applications point of view. 



3. .\AU3IS or L\^ZZ FL^XI^L" 3PACEr:"l\P? 


3.1 Intro-iurtion 

■bvinj da'/olope;l the necressary “r^ncepts for the decree of 
controllability of o systf*m w? are in a position to apply those to our 
funda:»ent.al problem, the study of the effect of location of actuators on 
/ery lar^e flexible spacecraft. To do this first '.•k have to obtain the 
state spice form (2.2-2) for the system dyaa.ncs aquations of notion for 
tlie Sioacecraft. ’Hjs follo.vinj arc the asssntial st=ps ineolvai in 
obtain in ;j the stat.s space forn. 

1. Choosinj a jenaral mcrlsi for a typical eery lar; 3 c flexible 
spacecraft. 

2. Derivation of equations of motion and their linearization about 
the equilibriun state. 

3. Translatin] the generalized forces into the normalized effort 
u, and introducing vAeightirrj factors for the coordinates to 
obtain normalized state variables. 

In this chapter we will be concerned v/ith the ano. ysis of a gcncnl 
model for a large flexible spacecraft and obtaining the state smare forn 
for tha dynamic equations describing tho system. 



3.2 Typical Planar 'totion s-todel of a ^tery Large Flexible Spacecraft 


The typical planar nodel consists of ossentially a rigii boiy B , 
and a flexible body 3p , with the rigid body ns a central core 
separatini the flexible body, 7nus, the flexible body consists of tvjo 
very large appendages attached to the rigid bOvdy one on either side as 
shovjn in Fig. 4. T,ve syste.'t is considered inortially at rest, i.e., 
any disturbance to the systan results in sone .iefometion of 8p and 
some notion of such tliat the center of rnass of the systen remains 
stationary in inertial sfMce. The systen will be treated as a continuun 
con^ the h/orid coordinate approvach [5] will be used. The flexiole body 
Bp will be treated as elastic betas and the deforn.'.tion in the plane of 
rotation will be assayed nonn-nl to their undeforned axis. 







^0 


3.3 Dynamics 

The tollowing notations refer to the model shovn in Fig. 4 . 

0: Center of mass of the system stationary in inertial space. 

C: Point fixed in with vghich 0 coincides v*ien there is no 

deformation of the flexible body . 

B: Frame of cartesian body axes b, , whose origin is at 

A A 

C. The plane of motion is b,-b^ , the axis of the 

A 

ofidefonned flexible appendages is parallel to b^ , and 

A 

their deformation is in b^ direction. 

Masses of the rigid body (B^) , the flexible body (B^) , and 
the system (B^+3p) , respectively, 
dm: Differential mass element in the system. 

R: fbsition vector of a mass element in the systesn in the 

inertial frame (i.e., from 0) 

R^: I’ector dravn from the center of mass 0 to the point C. 

: Position vector of any mass elenent in the system with 

respect to C. 

r„ : Vector from C normal to the cndeformed axis of B,. . 

r«: Vector in b, direction denoting tlie offset of the flexible 

component from the base of . 

A 

r: b, component of r for mass elements of B^ . 

y(r,t): Displacement of mass elonents of B^ relative to its 

undeformed axis. 

0 : Inertial attitude of B^ (and rotation of body frame B) . 



Time derivative of (*) with respect to the ir jrtial frame 




•(*) : I 

/ 


Time derivative of {*) with respect to the body frame B, 


From ttiese we have the following relations; 




Now, let 

. y (^y t) « y(Y,t) . T « y/i, -c Y^' 
In * ^ * r W, 


(3.3-1) 


(3-3-a.) 


so that 


3 » Yn -t- Y 3 ^ k, d" C ^n ®F (-3- 


From these Lhree equations, we can write 

^o-r 8^ 
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A ' ^ 


®r/.o 


■f*Y Be 


■fwT ®n ®C 




•'r 6 « 

■(•T 8 { 


i f*T 0 p 


e-3-^) 


The expression for can be obtained by considering the center of mass 
of the system with respect to the body frame 3 as follows: 


— »n, Rt « f I'dm c f Y'4*n. f 1 ^ 4 m 


C3 - 3-5.) 



? 2 - 


where (3«3-l) has been used for^^ for 8^ . But when there is no 
deformation of * 0 (see definition of C ) , i.e., 


j y* J C **■ I ) 


Using this in (3.3-5) we obtain 


Res «^-J- (3.3*7) 


v^ere y has oeen substituted from (3.3-2). Let 


Re--— ^ Rc « Rc 

'^S ) “ 


Hence, from (3.3-1) ,(3.3-2) 


R « Rg + Y ' , Rt i- V,' t, + * r,' !>,-*•( Rt 


R -®r,' + 

(B) i-j Q V ^ 

= -© CRci-\ ) C, -- ' f K 


wlieie *r^^ » 0 from (3.3-4). Hence, 


. "'i ■ (S'L-CnoV W)V » 

t ^ ^ (*V ) + ( ® Rc ^ (.3- 3- '0 



3,3.1 Kinetic Energy 


Ihe kinetic energy, T, , of the system is 

B i I * **' A ('3.3-^^) 

v^uch on substitution of (3,3-11} yields 

Xy*fcf« ' *yifet«v 

( cv)’’**"! t®\ (V-*- 

^ '’jyiN-v 'iyrW*^ 

(,3. 3->3; 

Since r/ , r^ are different for B and 8 (see (3.3-4)) the integrals 

I X ^ It 

over tne systen that contain these have to be separated into the 
integrals over aod 8^ . Hence, substituting for r^^ ,r^ , etc, from 
(3.3-4) these integrals can be written as; 

sy»te^ ^«P Bp 

1' 


(«v 

j Tr J » 


^YSbtw>. 

OiD 

} ^'C*V 


*yikB. 

(v; 

/ * 


sri^ 


®P 




•rytle 

• (vtl f V ) * f 9 






Where ™p * <** » ^ * (r')^+ (r^)^ has been substituted for . 



Define the following system cx>nstants: 




Sf, * 


J 


/ 

• J -V • j 


C 3 - 3~»5 ) 


jyite»s 


in wich represents the system moment of inertio about point C vjtjen 
there is no aefornation of Fp (y * 0). Subfft'tuting (3.3-1**) in 
(3.3-13) 'JT»n 'jsing (3.3-15) we obtain 


X * e 


^ * 6p Sp 

® [<=.<' «J + j " i -'j ^ i f 


3c 


ivc 


(3- 3-ifc) 


v<here = (da . 


tysiimwm, 

At this point modal analysis *;ill be adopted for the appendages. 


Define 


3t.V,t-) . $Cv)>^(V> 

tlw) =[(*,(-), 

. [ncv), , 


(3.3-t?; 


where ^ , *1 , j»l»2,...,N, ore scalars, which are functions of space 

J 

an! time, respectively. With this definition we have 



fr 


j ' [ i f 1 lk> ; f - {(«[ / (,)4„] ^ IV) 

*» •' H % 

’ [/#Ct) 4.»] ^(t) ; 1 5 V . ^ Tt)f f fi;,f W4-J ^(t) 

jri*^ .[[rf{Y)d«] ijtt) 

•» •, 


(33-«0 


and usin^ (3,3-3) 


^t* --^ [ Jftr)4,*j.{(tr) 


»->V 


Also, define 


• C /^l , > /**N ] • / ^ C^) 

•f 

®F 

• [ A#, , , »A#^] • /^ 

^ (3-3-ao) 

Substituting (3.3-18) and (3.3-19) in (3.3-15) and usi,^ (3.3-20), the 
expression for the kinetic energy of the system becomes 

■^s-6 (i 


«4iere 


(3- 3-^0 







M « A — — ^ 

‘i»*^'“-^ (3 !-»>•) 

3.3.2 Potential Energy 

Since tne system is inertially at rest the potential energy in the 
systen is solely due to the elastic deformation of the flexible body 
. Let be the ptential energy in the system, and for elastic 
bea:ns we can write 

V, .J jCeiJ) (.’-s-i’) 

®F 1 

•f 

.j ((ex; {») [ 

= if j 

“ i ’T" 1 

where y is substituted from (3.3-17), denotes cue second 

derivative of ^ witn respect to r, (SI) is the stiffness of Lne body 
, and 

k . j (,EX) 

®F 





t; 




I 


(• 


t 


t 


t 


3.3.3 Generalized Fbrees 

lb derive the <jeneralized forces we will assume that external 
forces and couples act at discrete locations throughout the system. 
While forces will be assuned to act at points couples will be assumed to 
act on infinitesimal elements so that those can have rotation. Let 

denote the inertial velocity of a point at jth location in the 
system .here a force Fj is applied, and let denote the inertial 

ar^ular velocity of an infinitesimal element at jth location in the 
system where a pure couple is applied. If denotes the 

generalized force corresponding to the generalized coordinate q(t) , then 



(3.3-26) 


where ^ stands for a class of forces and ^ for a class of pure 
couples. The coordinates q(t) are and in our case. From 
(3.3**10), can write 







in which we note that as the location j varies r' , r^ and their 
derivatives are the only quantities that are affected. The inertial 
angular velocity can be written as 



(3*3-2t) 


where is the slope at jth location on relative to the uvileformed 
axis of . We have 





(s. 


t«o(^ . 




/ 

'*^*®*^* * (tj) denotes the first derivative of ^ with respect to r at 
r • r. . For anall deformations tanof.* •<. , and hence 

J <1 

s \W (3.3>30) 
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vhich gives 


=1 ® y 


4^ 0» 


^ «p (3.W0 

Now, from (3.3-19) and (3.3-20) 

i 'i {.33-31) 


and fran (3.3-4) and (3.3-17) 
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-r,i Ur a, 
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4«r B, 




(3. 3- *3) 


With these equations, from (3.3-27) and (3.3-31), for q -® and q « »|^ , 
we can write 


• j ( lir/tl S ^ 


81 


ut 


M 


f -< 


iS' . I . . f o •*' *• 




fj 

1} 


Lj ' ^ ti K 


f> 


(j.3-35) 


Ihe forces and ttie couples are applied in the plane of motion. Ihe 
equation for the generalized forces in (3.3»26) then becomes 
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^ ""i/' (tw 
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1^*3- ^ ‘ (3.J-3 


We see that 0^ has terms containing and we desire the generalized 
forces to be independent of the co<ordinates. Ihe nature of the problem 

A 

is such that the deformation is only in direction, and hence it is 
reasonable to assume applied forces to be in b^^ direction only. 



C3.3-5») 


Therefore, we will set L;« 0 so that 

jR • K 


With this the final expression for the generalised forces are: 
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• Z C ^*1 Vj + 1 i 

* *'^C (3-3-3^^ 


where 
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^ 4 fH Bp 

(3. 3 “40) 


3.3.4 Squations of Motion 

We have derived the kinetic energy , the potential energy 
and the generalized forces Qy* Ignoring the damping (this is 
reasonable for spacecraft ) , the Euler>Lagrangian equations can be 
written as 

and L • Tj “ . From (3.3-21) and (3.3-24) the Lagrangian 

L .T,_ V, . t’mI f'u/'Tl) i 

(3- 3-41) 


from which the equations of motion for the spacecraft are: 



+ ^ C^'V**'^»J (3-3-4j) 

Where M, K are synmetric and 3-, ^ are as given in (j.3-39). 

fy \ 

3.3.5 Linearized Gigaations of notion 

In deriving the final form for Lne generalized forces we already 

made the assunption that the defamations are very small. 'Ve will 

consider arbitrarily small deComations in the appendages so that 

nonlinear terns in the coordinates '[j of can be effectively ignored. 

Also, if the -efornations are small toe perturbations in the attitude® 

are also anall. Hence, nonlinear terns in the attitude perturbations as 

well as terns involving products of deformation coordinates and attitude 

perturbations can he neglected. Further, since our concern is mainly in 

returning the systan from a disturbed state to its noninal (eq-jil ibriun) 

state we can, without any loss of generality, assune the nominal state 

to be at rest, and will consider attitude perturbations about a mean 

value. Therefore, let = a 4# 

; ®*0*O (S-3-4®) 

v^ere Q is the mean value of the attitude and 0 its perturbation about 
this mean. From this 

. g«a ; 

Substituting toesa in (3.3-43) and neglecting nonlinear terms in the 
coordinates we get the following linearized equations of motion for 
the aystem: 



(3- S--40 


‘I**’/' ^ 


Or, in matrix form 
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JLM U 


O o 


6 


9. 


(3.3-<^7) 


Where 1 is an Nxl zero vector and 0^ Is the transoose of fi. 


If define the following matrices 
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X, /X’ 

/*. H 


k « 


o o 

0 k 


) 9* 


9 
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W J 


(3. S-41) 


then the equations of motion in (3.3-47) can be expressed as 
^ IL ^ ^ ^ 

H k « Q 


(^3- 3-4^) 


or 


R-*? 




v^iich destribe the system in the hybrid coordinates 9 » We could now 
proceed to (xit these equations in state space form. But, since these 
equations have some special property we will investigate an alternative 
form for these equations. ttote tnat tiie matrices K are real 
synmetric and ? is positive definite (associated witn kinetic energy) . 
Such a pair of matrices can be diagonaHzed simultaneously by a 
principal matrix. Often it is desirable to work in these normal 
coordinates so we will derive the system equations in orthogonal form. 





V* • 


Let V » [ V, 1 be the principel matrix whose colunns 

are of 'dimension N+1 such that 

(3 3-si; 

where E is the uiit matrix of or-ler N+1 and ^ is a real diagonal matrix 
of order N+1, Detcails of such transfornation can be found in texts on 
matrix theory, e.?. , [*?]. Hence, a new set of coordinates 

T * *^1 obtained by the transformation 

(33-Sx) 


Substituting this in (3.3-49) and premultiplying by V 
(3.3-51) we hoy/e for tne system equations 

coordinates in this case are often referred to as vehicle normal 


and using 

(3.3-S3) 


modes. 


Tho advanta*’? of (3.3-531 over (3.3-50) from the al<>:ebraic 
ooint of v'lev/ is aooar«nt. Before we discuss state snacs for+i 
for the system equations we will translate th*» generalized 
forces Q into normalized control effort. 


3.3.5 Normalized Control Effort 

In chapter two it «^s pointed out that the control effort r.u.'st bo 
normalized. The control effort is imbedded in the generalized forces ^ 
we derived in (3.3-39). The generalized forces are due to the control 
actuators. We will consider two types of actuators; 1) Force type, and 
2) Torque type. Hie output of these actuators will be forces and pure 
couples, respectively. Any real actuator lias a maximun output possible 
which we will term "actuator strenjtn", and it is this maximun which 



will be used to normalize the generalized forces. Referring to 


(3.3-39), let 


Vj ■ 


(3- 3- S'*) 


Where Y. 




» K" 


are the strength of the force actuetor and the 


torqpje actuator, respectively, and j is any location for these actuators. 
The normalize:! control effort is 'Jj(t). Nbte that Uj(t) is different 
for each type of actuator (i.e., j belongs to different sets '3r * C ) • 
Substituting (3.3-54) in (3.3-39) 

fre 

wiiere o^; is zero for locations on and unity for loc.ntions on Bp . 
Define the following matrices: 






“m 


n 
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* 

u.* 


. • 1 

1 



«:-i 


H.] i » 

>^.1 ; ('* i ) • ‘ 

( 3 ..‘ sO 


t 

#N. A 

In L'lis cas.j H- is a row vector an-i H. , H.^ ara M:<1 oolunns. with 

9 

thosj notations (3,3-55) can a written as: 

t 

<?, • H, 

n* 

4 

and fro:r\ (3.3-43) 

§ * HU, (3.3-js) 

T!ie generalize'’, forces are now in Lhe for.n of nornalized control effort 
u. Since we have the total nunber of controllers equal to -n, s * m 
in the above equations. Ihe -natrices are of dimensions l>on and 

respectively, so that H is of dimension (J+Uxo, And u is of 
dimension :nxl. 


3.3,7 Normalized State Space Form for System Equations 

will consider the state space form for the equations of motion 
(3.3-5^) and (3.3-53) in the two systems of coordinates \ and ^ , 
rasp-ectively. 

1) sytem of coordinates 
Defining 


I 



(3-3-51) 





t • ^ a, n, 


\ \ 


the equation 

•- ^ I N St 

+ ic)^ " M « 


M V4U 


ran bs put in tlie forn 


A** 


(3. 3-60) 


where is an augmented matrix whose elements are elanents of "m'*'k 
anJ zeroes and mity at appropriate slots, B* is an augmented matrix 
whose elements are elanents of ^*'h and zeroes, such that (3,3-''')) 
reprosents the system equations (3.3-50) and identity relations for all 

9 

the elanents of . 

2) ^ systan of coordinates 
Define 

A*''*Lt t i, ^ »-4') 

Ihen tl'^r equation 

+ A ^ ^ 

can oe put in the form 

X* • A*X* + e'W' (3-3-<»l) 

where A* , B* are augmented matrices, similar to (1) above, with 
elements of A and V**H, respectively, such that this equation represents 
the system equations (3.3-53) and identity relations for all the 





elements of ^ . 

lb obtain the final state space form oorresponiing to either of the 

^ - 

two systems of coordinates, ^ , normalization of the state >^ariables 
X* roust be effected. This is done by assigning weighting factors to 
these >/ariables so as to reflect the relative importance we attach in 
controlling taeso. Let N- , N* , i»0,l,2....N, be a set of chosen 

>eighti.*i-g factoris so tiat 

ft * 

. * 




N- 


7J, 


, t so, ',1, - 


U 3-‘0 


f . 


vhere X * , :<• , j»l,2,....2'J+2, are elements of x*, x, respactively. 
For toe ^ system 


<.» ; *»’• 4 1; 


*< 


and for the systen 

X,. « €. • X,* * t. » • ®/ 


»,i, / w 




(3*3- ***) 


From (3.3-33) we can write 


where 







( 33 - 

44 ) 


N, N. N, 

N, - - - 

■ 

J 

(3. 3 - 

47) 

for X * 

from (3. 

3-55) in 

(3.3-33) 

or (3. 

,3-52' 

and 

r. by d;* 

we obtain 

the final 

state space 

1 form 



X • Ax 

+ 6 n 




(3.3- 

49) 


where 





^ 1 ? 


A = oI’a* d 


e 


-4 


Oa'B 


Thus, the e^'^tions of nnotion in either of the two systems of 
coordinates ^ ^ can be put in the normalized state ^?e form 

(3.3-S3). It should/ however, be noted that the factors , 13^ are 
not the same for both the systems t ^ f • The coordinates *l are 
physically more meaningful than ^ which are abstract (a combination of 
hybrid ooriinates $ , 1) and difficult to interpret in physical terms. 
As such it is easier to assign the %«ighting factors for the ^ 

systian cnan for the ^ system. Nevertheless, it may be possible to 
translate a set N. , in ^ system to an equivalent set ‘A^ , N; in 
^ system using some understanding of the system behavior, etc.. 

4 

However, in tiie present analysis we will assune that it is possible to 
obtain a set N. , N- in either system of coordinates. Mso, note that 
any one of the .>!• can be set equal to isiity since the weighting factors 
reflect only the relative importance among the coordinates. 

Before ending discussion on equations of .motion we will derive 
expressions for the mode shape integr:.is based on this model. 


3.3.3 iode Shape Integrals 

K ti f 

The mode shape integrals A; 'A ‘-1^2 H, 

3*1,2,. ..,N in (3.3-23) and (3.3-25) need to be evaluated to obtain the 
consta.nts of the equations of .motion. These integrals are taken over 
the flexible body which is considered as cne body, and hence any 
mode applies to the entire body. Because B^ is split about a rigid 
body these integrals are evaluated as explained in the following. 


V 





z 


' I. 


\ • 


t 



A mode shape is chosen separately for each appendage whose fixed 
base is tlie rigid body, and these two mode shapes for both the 
appendages togeLher describe a mode shape for the flexible body . 

Thus, there is a discontinuity in a mode shape of Bp due to the rigid 
body. Let f. (r) denote a mode shape of the flexible body 3p , and 
define 

I Si r « 1 

(3.3-70) I 

T * A : 

vtsicre ^.(x) applies to the appendage in -b, direction (left) an! ! 

A 

applies to the appendage in +b, direction (right); r is as defined 

A 

before anj x is neosured from the rigid body base along b^ direction; 

.d • 

£ , L are tne lengths of the appendages (from the rigid body basa) on 

* 

the positive (right) and negative (left) sides of the b, direction, and 
r^ , r^ are the offsets of these appendages about the btse of r^ in 

A A 

tbj,-b| directions, respectively. 

The mode shape integrals from (3.3-20) and (3.3-25) can now be 
expressed as: 

O) = I -h I 

' S V 

W A(i . 

A ' < 


f 



(CO 


^ ^ C.M,... **, ^>1, 

where b' , are used to denote the left and right appendages of , 
respectively, and (SI)*, (EI)^ are used similarly. Let 


/ ®p 

lc/Cv)4r 4^ Sf 


(3.3-72) 


denote a differential mass, where c” (r) , (r) ere linear mass 

9 W 

densities for B* , B^ , respectively. Note that r can be expressed in 

terms of x for 3* , b 1 as follows: 

P F 


y- (x- V , ^ s- 

x**> V , ^ r ^ 




fro.n wliich 


Cnr + V 

t X- *»« 


• X s» s. 0 


, »F 


0 S X S A , ®F 


With these we can write 

-Tie 


j(f")aw, . I (^-;c;cr)av * I (0<i 

Hi>v) Jc 

i\rt 


.'(T)ax 

r 


. J • J*(-f 


(3.3-73) 


(3.3-74) 


(3.3-75) 
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o^ere f " , f ^ are any functions corre^ndin^ to b;, 
respectively. integration of the mode shapes can be put in a 

nondimensionalized form as follows. 

Define 

Y • h , f *X* X -frr 

^;CY)« €; ^ 




^ «P 


* ClTe ^ 


Ye ^ 


(3.3-7-;) 


vhere is some constant having a dimension of length and 

are nondimensional fanctions. Tnese lead to, from (3.3-73), (3.3-70) 

and (3.3-72), 


*Y* * 


7C*-Cy,’) V ^ ^ 4w « fv 

7t**vCV/ 0r ' |c(r)raY^ 

C*) ■ 6j C«*) j - 1 i » S. 0 

^^.Cx) . ^/CxN , 0 * ft » 

jt X Jl f \ ^ ^ * C^*) 


(3.3-77) 


4 : f 

V 1 



In nondimensional fora (3.3-75) can be written as: 

/ O 

I C^")***^ = f i 'f') Cj(.y) ix = [ Cdx* 

^0 0 


Jsing (3.3-75), (3.3-77), 

appropriately in (3.3-71) 
mode sliape integrals. 


(3.3-75) 

(3.3-73) an-J substituting for f , f"** 

WB obtain tiie following expressions for the 


+ ( c;^/u*)Ax‘ 


- » 


(iO 


(fii) A, 


O U 

• f «r - fcj 

+ I c/ tV^*(x')d»* + J 

V . 
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(3.3-79) 
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where 


Cf * 6. C'it)l 

c^ * €j (j<) i 

4 ^ Z* . ^/"« 
rj e Tj y 


ax- 


(3.3-30) 


The expressions (3.3-79) will be used to evaluate the mode shape 
integrels for any chosen mode shapes. In the next chapter we will study 
a s^cific model of a larje flexible spacecraft and obtain nvmerical 
results. 



3.4 Sicnmary 

In this chapter we analyzed the dynamics of a tyoica] 'node] of a 
very large flexible ^cecraft inertially at rest, and obtained 
equations of motion, linearized these equations about an equilibriun 
state and put them in state space form to be studied from the 
controllability point of view. The equations were derived for planar 
motion, and the flexible body was assumed to be elastic beams undergoing 
transverse deformation. Linearization was based on the assunption that 
the flexible body undergoes small deformations. 



4. AFPLIOVTIQNS 10 A SPeCIFIC MOOBL OF A LARGE! PC£XIBLE SPACECRAFT 


In the last chapter a general model of a large flexible spacecraft 
was analyzed and equations of motion were derived. In tiiis chapter we 
will choose a specific model in order to obtain nueerical results and 
apply the c^^ncepts of the degree of controllability developed in chapter 
to study the effect of actuator locations on the controllability of 


the system 



4. 1 Model Description 


Figure 5 shows a specific model of a large flexible spacecraft. It 
consists of a cylindrical rigid body and two idantical appendages, i.e., 
the length, linear mass density and stiffness of both the appendages are 
identical, the latter two properties being uniform throughout the length 

A 

of these appendages. Ihe offsets of these appendages about the axis 
(in the rigid body) are equal. Each appendage is treated individually 
as a cantilever bean with its root in the rigid body as the fixed end. 
Any cantilever mode considered for the tvo appendages taken to-gether 

constitutes a single node for the entire flexible body 3p. ^or 
any cantilever node of an aopendats two tynes of nodes for the 
flexible body 8p will be considered: 

1) Symmetric mode — which a cantilever mode of an appendage is 
imposed symmetrically about the rigid body for the two 
appendages, i.e., the appendages execute a symmetric mode of 
motion as in Pig. 6. 

2) Antisymmetric mode — in which a cantilever mode of an 

appendage is imposed antisytnnetr ically about the rigid body for 
the tvo appendages, i.e., the appendages execute an 

antisymmetric mode of motion as shown in Fig. 7 . in this 
analysis ve will treat separately the symnetric and 
antisymmetric ;nodes of motion of the flexible body . 


4.2 Degree of Controllability 


4.2,1 Hode Shape Integrals 


Fbr the model described above, we can write, in the notations of 
section 3.3.8, 

f. I’-l •, V* Y," 

m ir) » C, « 

Wp »2Vl 

^.Cv)s(JkCx) ^ Y ^ 

^ V ^ A . r« 

^.U) -iix*40 

* * ^ 4 > 

f?.cx) s c. , 0 * x*$ 1 

C. « » Cj * 4« 

’* ' ' *• *^ * mp 

^ • i w 


. (4.2-1) 


vhere 4^ , apply to the left and right appendages, respectively, and 
the asterisked quantities are nond i mens Iona 1 . Usin; these in 
(3.3-73) the exoressions for the mode shape Integrals become 

« <^i£ / -I- I 

H • 



^ 0 I 

-c.Ay/[ I - I f*i^*)^^**] 

-I 0 

(m“0 a., s j j cIhJ 


*<{ 


_ <.Ex; g. 




C, <1 


( 4 . 2 - 2 ) 


The cantilever mode shapes [^] for an appendage can be written in 
nond imensicnal form cis 

» (ce*ft.pjl X*— X ) (si‘nC|?.4x— SJn.^lx*} 


S X* eg 1 




( 4 . 2 - 3 ) 


khere f. , "/• are constants, and the mode shapes are written for the 
right appendago (B^ ) . The modes for the entire flexible body are 
obtained by choosing 

S <T , -> * X* AO 

J J ' 

S 1 «y»wine+w*t t^tJe ef ®f 

1 Mymmrnkrtt, t*»o4« Of Bp 


( 4 . 2 - 4 ) 


•J*' vryiiw; 


' f 


The mode shapes together describe the mode for . 

Making a change of variable, x* « -i*, in the integrals involving 
^•*{x* ) in (4.2-2) these integrals can be expressed in terms of the 

A 

integrals involving (x*). And the mode shape integrals in (4.2-2) 
can be expressed as: 

I 

o’) s / 9*i^*^^"*-* 

w Av,. • + 


iiV.) K 




(i+O / $r (X*) A‘cx*; <tx* 

i* ; ' 


•■i J 


(4.2-5) 


The cantilever mode shapes (x) (or (x)) are orthogonal 
functions with respect to the mass density c^ , and hence the following 
properties (see [5] for details) can be obtained: 


id 


J CX) rf. CX) C^Jx . «»<j 










(4.2-5) 


liMre is a normalizing constant to be chosen later, toj is the 
frequency for the j cantilever mode, and is the Kronecker delta. 
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\\0 


These properties are tr'je for (x) alio. In nondimensional form 
(1.2-6) yields 


o(x* « ®<; 


(4.2-7) 


Now, define 


e. . / x*^/cx*;«tx* 

J / J 

e 

9j • 


(4.2-3) 


vhich jy the substitution of (4.2-3), can be obtained from the 


following : 


cc. « L- 


^ |1.i 

. j 

e . — L-- 

j ^ . 

b 

ft4 

■t 

-1- 

1*1+- 

1 1 

jl.l ' 


+ ) Sint jm + /j 


(4.2-9) 
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Llsiivj (4,2-7) and (1.2-y) in (4,2-5) tha Tiode shnno intacjrals can he 
expressad as; 


(.1) 


- C; (.1,4; CXj 


i'i; 


= CjHi- .') ( ®i + V o-;) 


{iff) 


= acj 6j t; •'ii 


(I’v) 



(4.2-10) 


Uiere and ^ * 1 have been substituted. 

AlthoLTjh one could consider many mode shapes will restrict hare 
to two (j*l,2) cantilever mode shapes of an appendage, the first and the 
second, Tiiis gives four inodes for 3p , two symnetric and two 

antisymmetric^ which are treated separately. iv}uations (4. 2-lU) 
represent t\vj sets of mode shape integrals for Bp , one for the 
symnetric mode (^= 1) and the other for the antisymmetric mode 
(^= -1). 

Table - 1 shows the various constants related to t.ae first and 
second cantilever modes given by (1.2-3) for j - 1,2. Tie nodt-s for a 
mode are denoted oy :<^j in the table. 


4.2.2 Computation of Inertia ami Stiffness '^latrices 
From (3.3-22), 


M = - 

s/’-r 


^ /V 

rhs 



(4.2-11) 
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If , /^ii denote the eleoients of H, /i , reepectivelyr then 
substituting for A;j , from (4.2-10), and Sj from (4.2-1) 


we obtain 


X. ^ 

Mii * - -^=v- (.<■**) ac'X; 

»np * fhs 

ic,, = 


(4.2-12) 


We will now use the freedom to choose the normalizing constant C. (see 


(4.2-3)) so as to obtain some algebraic advantage. If w» let 

*^04 * , j - a. 


* 




(4.2-13) 




j * I, a> 


(4.2-14) 


Substituting for Cj from this in (4.2-12) 




( i-ei J <lc 

“srsT^ 


k.. « T 4 J!»i a / 

Pj 


(4.2-15) 




The data for the model of Fig. 2 are given in table 2, and Table 
3 shows the required input data derived baaed on the assiznptions for 
this andel. Table 4 shows the frequencies m. of the cantilever modes 
and the constants 0* for both synaetric and antisymnetric inodes of . 
Using Tables 1, 3 and 4 the elments Mjj , /ij , K- in (4.3-5) 

are evaluated for synaetric (^*1) and antisymmetric (i« -1) modes of 
and result in the following matrices / K« M, K: 

1) Symmetric mode of motion of B^ (/• 1} 


M s I. 


l- ee o 


-3*IX4U« to 
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^"■=1, [o O] 
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S-i, 


-1 


3 . 1440 s* 4 10 ” O 

0 l-1^33t4 
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0 1.144eSl««o 
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(4.2-16) 



2) Antisynm«tric mode of motion of Bp -I) 
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whore M, H are the inertia and the stiffness matrices and are defined in 
(3.3-48). 


4.2.3 Computation of Matrices for Generalized Fbrces 


From (3.3-55) 
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in v4iich i refers to the actuators and is defined in (3.3-40). 


Also, define 
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vuhere k^. , are some constants related to the actuator strengths. 
ctow, 
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Substituting (4.2-19) and (4.2-20) in (4.2-18) «id using (4.2-10), 
(4.2-1) and (4.2-14) for /^^,t » Cj , respectively, we obtain: 
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^i^*CS " ^ «.2-22) 

*.r. Chi valuilof & D, , Oj tor our modil an jlvan in 

Tables 3 and 4. The distances r*^* , r^.* represent the variable 
actuators in B|^ , , respectively, measured in bj 

direction. Note that r'^ , rj* can be considered as one variable and 
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the factor takes care of the situation of whether the actuator is 

on Bji or , Substituting the values from Tables 3 and 4, 

(4.2-22) can be written as: 
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2) Antisynmetric mode of Bp {S» -l) 
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where (r * ) are to be appropriately substituted according 

to (4.2-21)^ and Q is defined in (3. 3*^8). The vectors , 

are completely determined for any r* . 


4.2.4 State Space Fbrm for System Equations 

Having obtained the necessary matrices we can now write the 
equations of motion (3. >-49) for synmetric and antisymmetric modes of 


Bp as 
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(4.2-25) 


\4iere K, ^ are obtained from (4.2-16) and (4.2-23) for the synmetric 
mode of Bp , and from (4.2-17) and (4.2-24) for the antisymmetric mode 
of Bp . Tb obtain the ortho?onal for*n of these equations (vehicle 
modes) we diagonalize M, K by the principal matrix V as explained in 
section 3.3.5. The principal matrix is given by 
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The principal matrix for the symmetric mode of Bp is taken as the unit 
matrix within the 'working accuracy of the nunerical terms. This means M 
in (4.2*16) for synnetric mode is appr ximated as a diagonal matrix (the 
off diagonal elements are much smaller than the diagonal elanmts) . Or, 
m 'X. E -for *f (4.2-27) 

< 

where E is the unit matrix. Hie orthogonal form of the equations, from 
(3.3-53), is 
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(4.2-28) 


where 

V ^ J 

A • v’^K ^ t \ '»] 

* (4.2-29) 


The coordinates ^ are vehicle normal modes. For the synmetric mode of 

Bp these are the same as If coordinates. Table 5 sives the 
values of for both symmetric and antisymmetric modes of Bp. 
Now, ) e t 



(4.2-30) 


where H is given by (4.2-23) or (4.2-24) for symmetric or antisymmetric 
mode of Bp . Then 
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where 
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o-f 3p 

* X ^ji f*'' 

■'•' (4.2-32) 

and v.j^ are elanents of V, ,'',lc»l,2,3, tal,2,...,m. From (4.2-32), 
(4.2-23), (4.2-24) 
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(4.2-33) 

where for syrometric modes of V;„* , and for , k* 1»2,3, the 


first or the second set can he used for any c (l,2,...,m) depending on 
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the type (force or torque) of an actuator. Of course, appropriate 
vectors , Z; must be used for s/imetric or antis^etric modes of 
Bp . Substitutirvj for from these we obtain: 


i s kpc (v,*Vc'+ v„wl‘V 


for force type of actuators for any i , t »l,2,..m 

iic - ( v,x + ''3x*a‘^) 

» kcc ( *► •» ** ) 

for torque type of actuators for any C , 


(4.2-34) 


where y, W; , Z; are appropriately substituted from (4.2-2S), 
(4.2-23), (4.2-24) for symmetric or anti symmetric mode of motion of . 
Also, note that a combination of force and torque types of actuators can 
be used (substitute ap(X‘opriate set for that c"'nunber; the total mmber 
of actuators is m) . 

Fbr the syinnetric case these reduce to: 
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(4.2-35) 
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Ttie state space form of the equations (4.2-2S) can now be obtained 
in the manner described in (3.3-61), (3.3-52), and for the system of 
vehicle normal modes ^ this is given by 

X* * A^x*-)- B^U 

(4,2-36) 
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Fbr symmetric and antisyinnetric modes of appropriate ^ and h^> 
are to be substituted from Table 4.5 and equation (4.2-34). ttormalizing 
the state variables as described in section 3.3.7, we obtain the final 
state space form corresponding to the orthogonal modes ^ as 
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S-4 


. «H b.i,^ the =«sen lighting factors. The value of 1 ,, i. „„ 
for both types of ^es of 8 ^ and g, .te to he taten fro. ibhle 
4.5 for the apfcopriat. node of B, . The eluents h,. , k-l 2 3 

^ *foUarly stf„titut«. fron ,4.2-34, usil, ,4.2-26,' 
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(4.2-23) arxJ (4.2-24). 

Note: This orthogonal form leadc to simpler A matrix consisting 
o£ 2x2 diagonal blocks which simplifies the alg^ra (the system can be 
separated into N+1 second order subsystems for computational purposes) . 
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4.3 proximate Degree of controllability 
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In chapter T*io, section 2.4.3^ • outlined tiio procedarc to obtain 
the spatial fonn (Oe£. 2.3) in real space for the state displa.jr.ait 

i , Following this procedure we will first obtain the Jordan canonical 
form J of the matrix A in (4.2-35), and also the transformation nutrix P 
and P**. Ihose are as follows: 

T s ^’ou^ ^ T, “* **r^v '-i^3 — '■J^3 



(4.3-1) 

vhere a , a.. , a,. , etc. are defined in (4.2-33). Ttje eigenvalues 
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are 0, 0, ±i JX, , li Ja, and the Jordan blocks are J, , J^ of 

order equal to 2, 1, 1, I, 1, respecti>/-*''v. *I!« spatial form is 
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'tow, replace the complex pair of partitions in (4.3-2) by real 
partitions from (2.4-25) as wcplained in section 2.4.2, to obtain ^ in 
real space, i.e., 
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(4,3-5) 


.Mow« since all the JOrd-'A blocks jr^ are associated with distinct 
eigenvalues, ■ 1, j»k*l ,2, ... ,5, w»uch fron theorem 2.3 gives the 
necessary condition m ^ I so thit the system may be completely 
controllable (in is the dimension of u(t)) . Ihere is no need to consider 
the rank test (necessary and sufficient condition) for complete 
controllability because we are going to proceed with the derivation of 
the approximate degree of controllability ^ , and from theorem 2.r> we 
know that the system is uncontrollable if and only if goes to zero. 

Die r.h.s. of o((t) in (4.3-S) can be written as 
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these, we 
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v*iere a,^^, , c, , hav/e been substituted from (-1.2-33) and 


(4.3-1). 


The aooroximate recovery region <*(6 dimensional oarallel 
ooioed) is ?iven, accordin'; to (2.6-3), as 
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where, as in (2.6-2), 
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From (4.3-5), evaluating (C^)” we obtain. 
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Applying the maxlnization to <j(t) in M.3-9), we con obtain: 
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rns Toraal discanaaa to tha surf’iaas of tho parallelopiped region 
are obtained usinj (2. ‘5-11), i.e.. 
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v^ere / are given in (4,3-14), (4,3-13), respectively, and 

A 

h^. are given by (4,2-34). 

Assuming that T is long compared with t!te period of the sine wave 
(cosine wave) (so that the effect of partial completion of the final 
(A»riod of oscillation of the sine wave is negligible) , the absolute 
value of the sine or cosine wave can be replaced by 2/jr * a verage 
over a period. With this substitution in (4.3-14) and substituting for 
IlgJI from (4.3-13), (4.3-15) becomes 
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in vjhich (4,3-1) and (4.2-33) are used for c, , z^, a,^ , a^^, and a^^ . 

“ria ajjproxinate decree of controllability based on the 
parallelooioed bound can now be written# accordin*? to (2.6-12)# 

dS 

* min j , J ‘ 

j (4.3-17) 

where the d. are given by (4.3-16). lb evaluate we need to compute 

A 

Ih^J, K»l#2,3, i«l,2#...,m# wlUch are given in (4.2-34). Hie 

information an the type and stren^gtii of actuators and their locations 

A 

is contained in hj^^ , TO put (4.3-16) in a more elegant form, define, 
in (4,2-34), the following: 

V.K 

r (Vi*) . V,. + 

V »i / IK a.*t I jw i (4.3-18) 

K « 1 C • •# 


where i refers to a particular actuator. If we pick a certain location 


t* for some actuator i»l then the right side of (4.3-18) is a number. 
Although we will use a finite nunber of actuators at discrete locations 
we can piclc any value for r * within the available range since our model 
is a continuun. Hence, the index i can be dropped in the above and r 
considered a continuous variable. Also, as mentioned in section 4.2, 
r'* and r.* can be considered as one variable since the change of 
location from the rigid body to the flexible body is taken care of by 
0^^ in 'A* and . Thus, there is one continuous variable r in bj 
direction. The functions in (4.3-lS) then can be written as 

CK ^ tk Xfe « I ( 4 . 3 - 19 ) 
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In the above the index i vnich refers to a particular actuator and its 
location along r* has been dropped. The functions correspond to 
force type actuators and the functions Q correspond to torque type 
actuators. Let 
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Then from (4.2-34) and (4.3-18) 
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We will take the constants k^^ , to be positive (see (4.2-19)) 

viiiich than leads to 
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Tnese are to be substituted in (4.3- 

•lo) to compute the ij .i^iose .minimun 

value, according 

to (4.3-17), 

is 

the approx imate 

degree of 


controllability f of the syste.n. 


In the next chaoter we will discuss 
examine the effect of actuator locations 
controllability of the system. 


the results and 
on the decree of 



5. RESULTS .VND DISCUSS TUN’S 


5. 1 IntraJuction 

In chapter Four derived expressions for the approxiin«ta decree 

of controlleibility based on a S[> 2 cific .nodul of o larje flexiola 

spacecraft. In this chanter we will disruss tne results for single and 
.•nultiple actuator distributions and exanine their effect on the 
approxioot-e degree of controllability 
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5.2 Influsnce Carves for Aotuator Locations 

The exprassions for tha normal distances in (4.3-1') can be 
written essentially as products of tlaree types of functions f , / » 

and fj where f, is a function of total time T, f^ is a function of 
weightinj factor and is a function of actuator strength and the 

effect due to its location in the system. Or, a typical normal distance 
d con be written as 


= ^ ^ 
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It is not important how we define only significanca 

is that T belongs to f^ , the ;«ighting factor belongs to f , and Lhe 
terms rol're.-J to the actuators oslong to f^ . Hence, in (4.3-15), 
noting (4.3-24) cluooa 

,11c, ifi', 'lO 
- 1 1^:1 

'*3- '*(T 

k* 1, 2, 3 (5.2-3) 
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Then, the nonnal distances d . in (4.3-lS) become 
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Ic» 1, 2, 3 (5.2-4) 

viiiere is defined in (5.2-3) and the rest need not be defined 
precisely at this stage. But ve will include the \ term in f^ so that 
the dimensions of f^ will be t!ie satne in either form (see (5.2-2)). 
The factor ^ is !<novr» because it is a systen constant. Hence, the 
unknowns in dj are total time T, weighting factors M, nornalited 
actuator strengths k^ , k^ , anfi the influence factors If^i, If^l due to 
tlie locations of actuators. All of these must be knowi before can 

be computed. 

If we are given a distribution of actuators we can presume tha type 
of actactor (force type or torque type) at any given location and also 
its strength (k^ or k^ ) . From the knowledge of only tha type of 
actiator at any grven location r* tlie I I (Ir^l or ll\l) factors are 
determined, with the additional knowledge of t'la strength of these 
actuators the f-inction f^ is deternined. It is reasonable to preside 
that we may not v«nt to change the weighting factors often. Hence, 
could also be detetmined if we decide on a set of weighting factors for 
the systesn states. This leaves total time T as the only remaining 
variable of d^ , so d- can be plotted as functions of total time T, 
And hence, ^ can be evaluated for any T. 
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The procedure to detemine th<» approximate ie^ree of 

controllability can be sontnarized in the follow' nj steps: 

1) Choose tine locations r,^, t* c* where the actuators are 

to be locate!. Classify them into the two groups, me for the 

force type ac^■Btors and the other for the torque typr 
actuators. 

2) Tabulate i fpit' 3nd I values for these locations 

corres Jordinj to the appropriate jroup ani for the type of mode 
of motion of Bp (symmetric or antisymmetric) . 

3) Compute the contribution of each actuator by taking the product 

of its strength (kp , k^ ) and the appropriate 1 f* I value 

tabulated in step 2. 

4) Take th* sum of all the contributions in step 3 to obtain the 
function 
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5) Choose a set of weighting factors 
functions 


ind compute the 


C • O; I, V 


in the expressions for Jj , j*l,2,...,^» 

j i 

Gi>/en a total tine T co.mpJte » f ^(r) in the expressions for 
o* f 
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7) Prom tiiese evaluate the normal distances as follo-s: 
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3) The approximate degree of controllability ^ of the system 
based on the given distribution of actuators is the minimum of 
all the normal distances d; evaluated in steo 7. 

j 

rius, by the above procedure (*^ can be evaluated for various 
distributions of actuators (type, strength and locations taiovn) , and 
tfiese distribution patterns can be r inked in descending order of 

the higher tne the better the corresporeiing distribution 

pattern, 'fhe decision as to *^ich of a jiven set of distributions is 
Detter may or may not depend on T. If a limited range for T is 
prescribed for control purposes (this may be the case in spacecraft 
control) then it is possible to decide from tiie curves the 

distriDution that is best on the average over this range of 7. 

At this point we can say our objective of ranking various 
distributions is fulfilled. However, we could gain a little more 
insight into the behavior of from the point of view of location of 
actuators from what we will call influence carves (I f* I curves) for 


actuator locations. 

By 1 r* I curves ws mean the two sets of curves If^l acid I fjl. The 
1(11 curves belong to the force type of actuators and che IJ^ I curves 

r • 

belong to tdie torque type of actuators. (Mots: these ace separate sets 

of carves for the symmetric and antisymmetric modes of motion of Bp , 
these two types of modes being treated separately.) The I I curves are 



functions of only th« ac uator location r*. A I f* I (I f^l or |f^ I) sat 
consists of n/2 curves If^l (IfJ^I or K*l,2, . . . ,n/2, w^ere n is 

the order of tho systcn, i.e., there arc as .nany component curves of If^l 
as there iue position /ariaPles (as opposaJ to velocity variables) in 
the original state v«actor (not nomalicec)) . (For our model n-^.' ^ach 
irjjl curve corresponds to a particular mods (in this case vehicle nornal 
•node) . 3ut tiie distances need not in -general correspond to any mode 
in thi original system, because the dj are aloni normals (t*> thv? 
surfaces of a parallelepiped region whose orientation in the state 
space depends on the s)iape of 

Fijures 8 through), show tiieso 11^1 curves. Figures 8 and )0 
represent wo sets of l/^l curves for force typo actuators for the 
symmetric .modo and antisyiometrlc .node of motion of , respectively. 
Figures 9 3.nd 1) rapr.sent the &»ts of !f^l curves for torque ty^ 
actuators for thj syanetric nolo and antis>nnetric .noio of motion of 
Bfi . These curves .are plotted as functions of tho .absolute valu-* I r*l 
of the actuator location. Die absolute value is used because these 
curves are dunlicateci for i r* (the appendages are identical and we 
have a symmetric physical systen' . Tf>e location r* is measured from 
the point C (see Fig. 5 ) along bj direction. For locations 

0 ^ lr*l 4 r* » ■ 3.03-13'’133, the actuator is on the rigid holy 

and for locations r* < lr*l ^ (I + r^ ) Lh? actuator is on the flexible 

body 3p on either of the two appendages. Tw absolute value lr*l is 

plotted along the x-axis and th-* I I or I (1*1 are plotted along the 
y-sKis. 

From (4.3-23) it can be seen that the contribution dui to each 
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Actuator is added up measured in terms of absolute values to obtain the 
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dj . liareo/cr, the contribution of e/ery actuator is a product of two 
terms# one involving only the strength of tho actuator (k^ » ^ 

the other involving only tne effect due to its location in the system. 
Ttia influence curves (Ip I curves) are very useful in computing the 
contribution to the degree of controllability of each actuator. 

Because these curves, being independent of actuator strengths# can be 
used for different actuators at the same locations# ervd the total 
contribution of a new set of actuators is obtained simply by adding up 
the products of tJvo new actuator strengths and the same influance 
factors for the old set. 

In the next several sections w» will e<a;une various options of the 
control actuators and their bearing on the .approxinate degree of 
controllability of the system. These ootions will be exanined for 

both types of nodes of notion of 9p , synnetrir and antisymnetr ic. 
'Phese options arc as follows: 

1) Single force actuator 

2) Single torque actuator 

3) 'Multiple Actuators 







S.3 Single Fbrce Actuator 


5.3.1 Intro>ioction 

The normal distfincas , from {4.3-1^) (m ■ I in this case) , are 
given by 
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where (4.3~24) has been use.i. For any chosen set of weighting factors 
total time T and for any strength of the force actuator 
each Jj is proportional to the height of one of the three curves of 1(^1 
(Fig. S for synnetric mode and Fig. 13 for antisymmetric noue) . 'Thus, 
if d. were to be plotted as a function of location I r*l we would have 
six curves, a pair of Lhem proprtional to each I Because of this 

proportionality the maxiituin, minlmun and taro values of any dj curve 
correspond to the maximun, minimun and zero values of its corresponding 
I I curve, i.e., the behavior of a curve is the same as the 



behavior of its corresponding pair of dj curves with respect to these 
values. Thareforc, we can examine the I fJ,jl curves to gain some insight 
into the behavior of the approximate degree of controllability as a 
function of location lr*l of the actuator. 

5.3.2 Symmetric Mode of 

Figure 3 shows the set of if^l curves for tha symmetric mode of 
motion of . There are four locations of I r*I at v#iich (•* goes to 
zero indicating the system is uncontrollable at these locations. 'Pie 
1 r^, I and irpi^l curves have one such location each and the 1 I curve 
has two such locations. These locations are tabulated in Table 5. 
Recall tiiat for the synraetric case the venicle normal modes are the same 
as the hybrid modes ( B . •[ ) # and hence the zero of any 

indicates that the correspond ir*g hybrid node is uncontrollable. ^t 
lr*l * D (point C) the rigid body mode ® is uncontrollable but the 
appendage noides are controllable. At this location the 

actuator can generate no torque on the vehicle in the inertial system, 
so 0 is uncontrollable. The symmetric motion of Bp (Fig. 3) produces 
only translatory motion of anrl by controlling the modes this can 
be killed. But if the disturbance includes 0 this cannot be 
controllad. The locations on at which I Ifpjl are zero arc 

i, , 

inertial nodes (obtained by putting ^ » 0) for the corresponding modes 
(first and second modes) of Bp . At an inertial node the effect of the 
force of the actuator is not fait by tha corresponding mode, and hence 
its shape cannot be controlled. 



Having ioaricad che zeroes of the irVi curves let us examine how they 
behave over the range of |r*l. For any location lr*l on the rigii body 
Bn (0 4 Ir^ < r* ) I and 1/^,1 are constant (Fig, 3), hence no 
improvement is attained in the controllability of the two modes and 
On the flexible body 3^ for locations lr*l beyond the zero of 
* '>.05) up to the crest of lf^,l (lr*l« 0.5) all three 
If;. I curves (see Pig. 8) increase monotonically. Likewise for 
locations beyo.vd the seconi zero of If^^l (lr*l c h.**?) up to the tip of 
the appendages all the three If;^ I curves increase monotonically. 
Hence, o/er these two ranges, 0.05 4 lr*l 4 0.5 and Ir^l >0.32 the 
degree of controllability (*^ of the systen improves as lr*l increases, 
'toreover, the individual maxiraun of each If^nl curve occurs at che tip 
of che appendages, and hence is maximun for the location at the tip 
o£ Cite appendages. 

In concltjsion we can state the following: 

Fbr the symmetric mode of motion of B|t , if a single force 
actuator is Co be used for attitude and ^ape control the best location 
for the actuator, from the controllability point of view, is the tip of 
the appendages regardless of the strength of the actuator, the total 
time T allowed for control, .and any set of v«ighting factor? one might 
wish to choose for tha normal modal coordinates. 

Othar salient features for the force actuator for symmetric motion 
of Bp arc: 

1) Locations in the neighborhood of the zeroves of the influence 
curves (Ifjil curves): Since at the zeroes of 1(^1 curves che 
degree of controllability ^ of the system goes to zero, 
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locations in the nsi^ghborhaod of thoso zeroes exhibit a 7cry 
poor quality of controllability of the syst.^nj an'^ shouli be 
avoided as far as possible. 

2) Two broad distinct ranges of I r*| on , 0.05 ^ lr*I ^ 0.5 

an-.! I r**l 0..02: 0/er these ranges controllability of the 

system improves steadily ( increases) as I r*! increases. 

3) Controllability of the B mode (rigid body mo^de) improves 
steadily as I r*I increases over its entire range in the system 
(0 4 ir*i a + «J). 

4) Tho quality of controllability for the two modes and 

r 2 r.K*.ins unchango:! for locations on trie rigid body 3^^ 
(0 ^ I r*I 4 rj) . 

Tiese features can be usefully applied wven searching for other suitable 
locations for the force actuator thin cne tip of the appendages. Suth 
occasions iiighc arise due to p<v/sical restrictions and other pr.ectic^l 
considerations. In practice one would ha^e to consider e sufficient 
nuibur of modes in the model to give a goof) representation of the system 
dynamic baliavior before the actuator location decision can be mafio. 

Tius, for the case of symmetric mode of motion of Bp -ve have some 
valuable information regarding the quality of locations I r*l in the 
systern when a force actuator is to be used. 

5.3.3 Viti symmetric dodc of Bp 

Figure 10 shows the set of Ifjl curves for the antisymetric mode 
of motion of Bp . The zeroes of the Ifj^l curves are tabulated in Table 
5. For tha antisy.ninatric modes of Bp the vehicle normal modes are 
different fro:o the hybrid modes 1 (• ,^ ). But is similar to the 
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rigid body mode ^ and and are similar to the appendage modes 
1^, » There are four distinct locations tr*l at which at least one 

M^t goes Co zero indicating system uncontrollability. At tr*l « 0 
(point Q all ir^l curves are zero which indicates none of the three 
modes is controllable. As stated in section 5.3.2, at this location the 
force actuator cannot generate any torque on the vehicle in the inertial 
system, and hence B is uncontrollable. The antisynroetric motion of 
(Pig. 7) produces only rotary motion of , and since this is 
incontrollable at lr*l « 0 the appendage modes Also 

uncontrollable. The locations on B^ at which I are zero are 

inertial nodal points (obtained by putting • 0 ) for the 

corresponding vehicle normal modes (the location lr*l » 0 is also 
inertially at rest) . At the inertial nodes the effect of the actuator 
force is not felt by the corresponding normal modes (hence. 

It is uncontrollable. 

As in section 5.3.2 we will now examine the nature of the If^l 
curves. For locations in the range 0 Ir*l 4 : ^ (in the rigid body 
B^ ) the Ifp I curves are all linear which indicates improvement of 
controllability for all the modes away from |r*l • 3 up to the root of 
the appendages. On the flexible body B|g for locations beyond the. zero 
of ir^^l (lr*l r 0.25) up to the crest of If^^l (It^l • 0.3), and beyond 
the zero of If^^l (lr*l 0.32) up to the tip of the appendages all the 
three ijpul curves increase monotpnically. Hence, over these ranges, 
0.26 ^ Ir^l ^ 0.5 and lr*l ^0.82, the approximate degree of 

controllability of the system improves as lr*l increases. >g3ln, as 
in section 5.3.2, the individual maximon of each ifp^l curve occurs at 


ni.f 


the tip of the appeniages. Therefore, P* is maximun for the location 
at the tip of the appendciges. 

The conclusion here, for the antisymnetric case, is the same as 
that for the case of symmetric mode of motion of Bp, in section 5.3.2, 
The tip of the appendages is the best location for the actuator for 
attitude and shape control regardless of the actuator strength, total 
time T and weighting factors for the normal modes. 

Other salient features for the force actuator for antisymmetric 
motion of are: 

1) Locations in the neighborhood of the zeroes of the influence 

curves (If^l curves): Same comments as for the case of 

symmetric motion of Bp, in section 5.3.2, apply. 

2) Two broad distinct ranges of I r*l on S^, 0.23 4 lr*l 4 0.5 and 

lr*| > 0.82: Sams comments as for the case of symmetric motion 

of Bp apply. 

3) Controllability of all the ->onnal modes improves steadily as 
lr*| increases over the range of locations on the rigid body 

(0 4 I r*l < ) . 

As mentioned before these features can be usefully applied in the search 
for other suitable locations in the systan for the force actuator than 
the tip of the appendages.' 

Thus, here also for the case of antisymmetric mode of motion of 
Bp , some valuable information regarding the quality of locations I r*' 
in the system is obtained with on 1 y the knowledge that a force actuator 
is to be used. 



i; 


SuiTinarizirv 3 ,for a single force actuator control we can say that the 
tip of the appendages is identified as the best location for attitude 
and shape control for symmetric and antisywaetric nodes of motion of 
Bp . This conclusion is independent of tho strength of the actuator to 
be used, tlie total titae T allowed for control, and the weighting factors 
prescribed for tho modes. Additional information is obtained which can 
be useful in searchin-g for other suitable locations in the syston. The 
information on regions ihars an actuator siiould not bo located, and the 
information on distinct ranges for locations lr*l where improvement of 
controllability is guarcanteed as I r*l increases arc very useful. For 
the sy.imetric case, some inder standing of the quality of controllability 
is gained for the rigid tx>dy mode 9 and the ap^oendage modes *1^ , • 

Tie information for both the synnetric arKl antisymmetric modes of motion 
of 3p taken together provides a better understanding of the system 
controllability, and identifies certain narrow ranges for locations I r*| 
that will bo suitable as altornati»/es for the force actuator in order to 
effectively control both attitude and shape in both types of motion of 
Bp. 
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5.4 Single Ibrqua Actuator 




5.4.1 Introduction 


The noraal distynces d» , fron (4.3-lS) (m = 1 in this 


case) , are 


gii/en by 
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v^ere (4.3-24) has been used. Hiis is similar to the case of a single 
force actuator in section 5.3.1. The d* are proportional to the 1/^1 
curves. The rest of the argunent is exactly the same as for a single 
force actuator except replace by (the strength of the actuator) 
and I by I r\„l in section 5.3.1, and refer to Fig. 9 for symmetric 
motion of aiod Fig. 11 for antisymmetric motion of . 

5,4.2 Symmetric 'lode of 3p 





Figure 9 shows the set of If^l curves for the symmetric mode of 
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motion of . TTie zeroes of these curves are tabulated in Table 'i. 
For the symmetric motion of the vehicle normal modes are the same as 
the hybrid modes | *^ ) / and hence t)ie zeroes of the I I curves 

indicate that the corresponding hybrid modes are uncontrollable at those 
locations. One feature that is different from the force actuator care 
is that the two curves, IfJ^I arvi 1 1^1 are zero over the entire rigid 
body range of locations, which indicates the appendage modes t, * 1j^ 
are uicontrollable with a torque actuator located in the rigid body . 
The S mode is controllable. As before, it is to be remembered that the 
symmetric motion of B^ (Fig, 6) produces only translatory motion of 
Bj^ , hence a torque on is useless in killing the motion of B,, . At 
the locations on Bp at which I fljl are zero the inertial slopes 

are constant for the corresponding modes (these locations are obtained 
by putting the inertial angular velocity ^ of an elenent on Bp equ’l 
to zero). In tliis case, the roots of the appendages happen to be such 
locations of constant inertial slopes for both I and lil^l curves. 
And these locations are included in the range of locations over . 
There is only one other location on Bp for a zero of the I "urves, 
an<3 that is for the lfc»l curve. At the locations of constant inertial 
slopes the effect of the torque of the actuator is not felt by the 
corresponding modes, and hence their shape cannot be controlled. 

The behavior of these curves is more uniform than those for the 
force actuator. The I fl, I curve is constant throughout the range 
lr*l , so no improvement is attained in the controllability of the 9 
mode. The appendage modes are viicontrollable over the entire rigid body 
range. For locations Ir^ from the root (r*) of an appendagf. ip to the 
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crest of I I (Ir^l « 0.21), and beyond the zero of 1 ^1*1 (lr*l * 0.51) 
up to the tip of the appendages, the twc curves iiUi. i<« I increase 
monotonically v^iile the curve I is unaffected. Hence, over the 

ranges r*^ |r*l 0.21 and lr*l ^0.51, the approx inate degree of 
controllability of the system either steadily improves or remains 

constant as lr*l increases. Also, the individual maxiitum of the two 
curves l^csl occurs at the tip of the appendages. Therefore, 

if not decided by I I , reaches a maxinira for the location at the 
tip of the appendages. 

In conclusion we can state the following: For the symmetric mode 

of .notion of Bp , if a single torque actuator is to be used attitude 
and sh.apo control the best location for the actuator, from ts . 
controllability point of view, is the tip of the appendages regardless 
of the strength of the actuator, tho total time T allowed for control, 
and any set of weighting factors one might wish to choose for the normal 
modal coordinates. (The same conclusion as that for a force actuator.) 

Other salient features for the torque actuator for symmetric motion 
of Bp are; 

1) Bad locations in the neighborhood of the zeroes of the 
influence curves ( I I curves): Same comments as for the case 

of a force actuator, in section 5.3, apply. 

2' Two broad distinct ranges of I r*l on 3^ , c^4 1^*1 4 '^•21 and 
|r*| >0.51: Over these langss controllability of the system 

either improves or remains constant as I r*l increases. Over 
those ranges the quality of controllability of the S mode 



(attitude) remains unaffected by change of location of the 
actu’Jtor, but the quality of controllability of tho apoendcage 
modes ''Ij , (shajae) improves steadily as ! r*| increases. 

3) The quality of controllability af tlie ® nods (rigid body -node) 
is constant for -nil locations in the system. 

4) Shape control is not possible > locations on the* rigid body 


'.Hissu features can be usefully applied as described before (under fores 
actuator) in the search for other suitable locations. 

Thus, for the case of the synmetric moc)e of motion of '"ip we .av - 
some valuable itifornatior on the qumlity of locations I r*! in the system 
even witii only the knowledge thc<t a torque actuator is to be used. 


5.4.3 Ymtisynnotric doJe of 


Figure 11 sho'ws the Swt of If^ I curves for t) antisymmetric mode 
of motion of Sp . The zeroes of these curves are tabulated in "daol-e 5. 
For tile anti synmetric nodes of Bp the vehicle normal modes are different 
fron tne hybrid modes ^ (. 9 , 1 ), But is similar to the rigid body 
.node & and are similar to the appendage mooes 1j , 7^ • Hiere 

aro tnree zeroes in all, one for ifl^l and two for Ifcjl* ! r^l •» h 

the system is completely controllable although its degree may be very 
low. Recall tint antisymmetric motion of Bp produces only rotary 
motion of and this can be killai by a torque at I r*l ■* h thereby 
eliminating the motion of . At tho location? on Bp. at wnich 
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are zero the inertial slopes are constant for the corresponding 
normal modes (these locations are obtained by putting the inertial 
angular velocity -£k of an element on equal to zero) . At these 

locations the effect of the torque is not felt by the respective modeS/ 
and hence they are uacontrollable. 

Ihe behavior or the 1 1^1 curves is different compared to those for 
the syiTtnetric case for locations closer to the rigid body 3^^ . O/er the 
rigid body range of locations all the If^l curves are constant 
indicating no change in the degree of controllability of the system (or 
the quality of controllability of any mode) is obtained by moving the 
actuator (torquer) . The systen is canplevely controllable at any 
location on . For locations beyond the zero of IfJ^I ( lr*l e 0.13) 
up to the crest of 1 { lr*l ft 0.21), and be^^nd the second zero of 
I I ^ I r*l » 0.51) up to the tip of the appendages the two curves 

increasa monotonically -4iile the curve I is unaffected. 
Hence, over the ranges 0.13 ^ lr*l < 0.21 and lr*l J 0.51, the 
approximate degree of controllability ^*of the system either steadily 
improves or remains constant as |r*| increases Tie individual maximun 
of the tvn curves lies I occurs at the tip of the appendages. 

Tierefore^ (9*, if not decided by reaches a maximun value for the 

location at tie tip of the appendages. 

Ihe conclusion here, for the antisymmetric case, is the same as 
that for the case of symmetric mode of motion ''f in section S.4.2. 

The tip of the appendages is the best location for the actuator for 
attitude and shap<^ control regardless of the wictuator strength, total 
time T and ia»ighting factors for the normal modes. 



Other salient features far the tor^e actuator for antis/nmetric 
.notion of are: 

1) BaJ lo:ations in the nei^hhorhno'J of the zeroes of the 
influence curves (ir^l curves). Sane comments ns for all the 
precedin'] coses apply. 

2) T/io Proe.i distinct ran- 3 es of lr*l on 3p, 0.13 4 lr*l ^ 0.21 nni 
Ir*l >0.51. Dver these ranges control lability of the system 
either improves or renains constant as lr*i increases. 

3) O'Jslity of controllability of all t)w nomel nodes is constant 

for any location on the rijii body is constant for the 

systan Cor these locations) . 

Tiese features arc useful in the se.arch for other suitnblc Iccations. 

ilius, acre also for the case of the antisynmctric !aode of motion of 
Bp , some valuable infornation on tne quality of locations I r“l in the 
systten is obtained with only th-e knowleijc that a t.orqj? ectuotor i? to 
be used. 

Sa'anarizi.nq for a siii 3 le torque actuator control we can say that 
the tip of the appenJaqes is idantifie.i as the be.st location for 
attitud® and shape control tor syrenetric and antisymmetric modes of 
motion of Bp. for this model of a spnccxrrcft v/aich can only dehorn in 
the man.ner described by the t.« noiles of (obtained fro-n the first 
and the second cantilever modes of the appen-^.a^es) considered. 
conclusion is indepe.hdent of the strength of tho actuat-jr to be used, 
the total time T allowfd for control, and the wei-jhtin-j factors 
prescribed for the modes. Information .on re] ions of poor 
controllability is useful to check that an actuator is not located close 



to any of Infomation on distinct ranges of I r*l where 

improvement in at least shape control is possible is also very useful in 
the search for altarnative locations. For t!je synnatrir case, some 
understand in^ of tJje quality of controllability is qcined for the riqid 
body 'oodo 9 an! che apoendajc modes 1^. Tie infon*tion for both 
the sytnetric and antisyxietric .oodes of motion of taken toqet’ier 
pro/idcs a better understanding of the systaa control! rbility, and 
identifies certain narrow ranies for locations Ir^l that will bo 
suitable as alternatives for the torque actuator in order to effectively 
control both attitude and shape in both types of notion of . 

Tiere is a striking difference between the force actuator control 
and the torque actuator control. Unlike tire case of force actuator 
control, in torque actuator control there is one .node, a tytoe of riqid 
body node (attitude in t.ho c-osc of symmetric notion), .Viose polity 
of controllability is unuffcctei for any location I r*l in the systai 
( I f^, I is constant) . 
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5.5 SunKiary for a Single Actuator 

Fbr a single actoator, force typ-J or torque tyt>e, v-e gained 
valuable insight into the quality of locations I r*l with respect to the 
approximate degree of controllability of the system. Ibe 'Tinjor 

finding was that the tip of the aRjendages is unconditionally th-^ 'oest 
location from fie point of view of attitjdt-* and siiape control for either 
type of actuator. In other words, for our flexible system there exists 
no other location lr*l in the system at wiiich an actuator of either type 
and of any strength will give a better controllability over any period 
of time T for any type of (sy.xnetric or antisysnetr ic) inotion of Bp 
under any arbitrary ’weighting of the .naies, Ihis is a vary general 
result. Besides this conclusion, we also observed other significant 
features in the analysis for a single actuator. One of thase ^/as the 
detemination of bad locations for actuators. The other vjns some 

distinct ranges of locations I r*l , over *.hich the degree o* 

controllability of ch.? system inproves stea^Jily as I r*l increases, 
meaning the quality of control of all tiie :.k improves. Over other 
ranges the quality of shape control definitely improves os I r*! 
inersoses but the degree of controllability of the systea never 
Jctericrctes moaning toe quality of control of the type of rigid body 
inorJe remains constant. Tiese features on the quality of locations I r*| 
are extrenely useful in searc.hing for alternative locations that will 
give reasornb?'/ good controll’bility for tie system. 
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5.5 Multiple Actuators 


I S“ 


In section 5.3 5.4 va analyzed the situation in v/iich n single 

actuator is usa'3. Tia anv».lysis providad us with very valunbls insight 
into tho ioehavior of tho d[^ro'Ciinate dagree of controllability as 

the location |-y*I of the actuator is varied in the syst?n. will 

now study the situation in winch 'nore than one actuotor is used. 7ia 
ain is merely to extract guide lines for preliminary design of 
•distribution patterns for the actuators. 

As in (5.3-1) and (5.4-1) tne normal distonces -i • can be varitten as 



wiiere as b-ufore c, , c. ,...,c, 
weighting factors. Let us 
(5.S-1) as 


are constants 
denote the 


depending on T and tho 


par inthotic*! 


z'l expressions in 



( 5 . 1 - 2 ) 










Ob'serve that each ter.n of is nonneiative and represents an 

independent contribution of an actuator to the nornol distances d^ . 

A A 

Suppose all lip^l curves increase fron If, I to Itj^l then all 
tiw dj increase simultaneously vhen a force actuator is moved from 

^ i|lk 

l'f,l to Kj^l , Continuing in this manner ve place one force 
actuator at the best possible location L,. >»ote that our system is a 
continuum, but actuators are discrete and occupy only discrete 
locations. Hence, our analysis may suggest piecing actuators at the 
sauie locations. From physical considerations ■-« will assume that only 
one actuator will be placed at any location. This '-gould mean spread in-j 
the actuators as a chain. A^ain from practical considerations will 
assuns tney are separated by some arbitrary distance, and hence discrete 
locations in a continuous ran-je will be taken in this context. Nbw, 
consider t.he second force actuator and place it at the second best 
location L^, and so on, where these discrete locations are to be 
interpreted in the context described above. Tiius, for each .actuator m 
use tne analysis of a single actuator case, w‘uch amoirts to ranking the 
locations IY*I in descending order of 'Quality for a single force 
actu.ator. ^bw, fill t.hese locations vdth the force actuators. It is 
reasonable to distribute tha actuators as evenly as possible betvaien the 
two appends^es (tiiero is al'Mays a p'lir of locations, one on each 
appendage corresponding to any It^l ) . One other important observation 
should be made regarding v^aich force .jctuator should occp^sy which 



location. From ( 5 , 5 - 2 ) it can be seen that if we place the strongest 
actuator at t)i3 best location the contribution to the *Jj is further 
increased. Thus, if oeer a range of K*l , |y*| 4 Iy^I 4 , all 

die lf^„ I curves increase monotonically, and ^ rank possible 
Iccetions Lj , L^, in iescending order of quality then the 

actuators should be placad at , 

respectively, where K_ X k_ > k_. . 

F p 

i<ep20t this wuole process independently for t'.ic torque actuators 
using tn*.? curves for the case of a single torque actuator. Note 

that for tliis case the Ifci I curve is constant througfiout, and hence 
can Oe ignored in ranking the locations assuning that ’ve are interested 
in controlling both attitude and shape. Tie other tvjo curves 

can increase monotonically over certain ranges of ly*! . As in the 

case of force actuotors, the torqua actuators should placed such that 
tie strengti of an actujcor corresponds to the -juality of its location 
in the set of ranked locations (the stron^jest actuator at the best 
i»ssiole location) . 

:tere 'Mj Iiave analyzed the sets of forca and torque actuators 
iniepindontly. It is possible that this could ask us to place a force 
actuator and a torqua actuator at the sa.no location. As explained 
earlier the actuators could be separated by some arbitrary distance so 
that adjustments c.an be made to accommoilate an actuator from each group. 

Since we have used the analysis for a single actuator case to study 
the multiple actuator case tho method outlir,d here is applicable to all 
types of motion of the flexible bady Bp that are considered for the 
single actuacor case. lb consider both synnetric and antisymmetric 
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.'notion of it is only nacessary to take the appropriate ranges of 
fran the single actuator case, in this -nultiple actuator case we 
have been concerned only with the ranges of lr*l over which the If' I 
curves exhibit -nonotonic behavior. In ot.her ranges of Iy^I the actual 
behavior of (® is not apparent, ami actual corjputotion of (•* becones 
essential including tJie now important influence of the weighting factors 
(and the time T, if nonlinear) . Tho locations close to the tip of the 
appendages are definitely superior in terms of the resulting degree of 
controllability for a spacecraft whose flexible body motion is 
restricted to the two appendage modes considered. Finally, it must be 
laantionej ttiot tlie system becomes uncontrollable only if all tlie force 
actuators are located at the zeroes of a single I curve and all 

the torque actuators ace located at the zeroes of the corresponding 
curve. Which will me .e a pair of dj equal to zero. 

In tlw approach 'wj hav outlined for obtaining suitable locations 
for the case of m’jltiple actuators the following are the T«ain steps; 

1) Consider the force actuators and torque actuators separately. 

2) Use the correspond ing single actuator analysis indepandently 
for each type of actuators. Corresponding to the tyne of 
.motion of 3p to be controlled mark out the ranges of lV*l 
bver vvuch the iP* . curves exhibit monotonic behavior. 

3) In these ranges of \x*\ rank locations in descending order of 
quality. 

4) Place the force actuators at their corresponding ranked 
locations in descemUng order of their strength. Similarly, 
place the torque actuators at their corresponfHng ranked 



locations. In ths case of conflict of locations make 

adjustments to accorjinodate the conpeting actuators in any 
aj^ropriate manner. 

We can consider triis approach a logical way to set up preliminary 
distributions for a sat of given actuators. Further modifications and 
adjustments in t!ie placement of actuators can be carried out based on 

other practical considerations. If these considerations lead to 

|| 

locations in other regions Wiere the beha>/ior of (• is not obvious then 

/•* 

computation of is necessary with the inclusion of the weighting 
factors and possibly time T. 

Let us now consider a two actuator case as an example of a multiple 
actuator case. Following the steps vibo/e ve can obtain the tip of the 
appendages as the best location in the system. Since ws have two 
appendages the best distribution for two actuators is one actuator at 
the tip of each appendage. In a multiple actuator case vnth more than 
two actuators two of tham can be placed at thesa tips. Tie tip being 
the best location ^ould be the place for the strongest actuator among 
the force actuators or the torque actuators. 

As a remark, observe that the 11^1 curx/es for the torque actuators 
tend to flatten near the tip ot the ap.aendagos indicating tfiat moving a 
torque actuator near the tip does not significantly affect the behavior 
of nils fact may be used in making adjustments in the locations in 

case of conflicts between force actuators and torque actuators. 
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5.7 Suinmary 

Ihe procedure vas described to determine the approximate decree of 
controllability ^ for a 9iven distribution of actuotors, weighting of 
the nodes, and total time T, ‘P'isn we proceeded further to analyse the 
influence curves (iP I curves) to gain some insight into the behavior of 
f , wf.ile maintaining complete freedom in the choice of the weighting 
factors and the total time T. A single force actuator and control 
and a single torque actuator control were examined for both 
symnietric and antisymmetric modes of motion of Bp • We 
gained some valuable information from this analysis. One 
conclusion was that the tip of the appendages is the best 
location in the ays turn for attitude and shape control. 

This conclusion does not depend on the time T or the weight- 
ing factors. We also Identified regions near zones of the 

curves# which exhibit poor degree of controllability of 
the system. In all cases# there were two broad distinct 
ranges of It 1 over which ^ behaves monotonically so that 
locations |-y*) can be judged qualitatively. 

Having obtained some insight into the behavior of 

for a single actuator we next examined the multiple actuator 
case. We showed that we could essentially treat the multiple 
actuator case as an extension of the single actuator case 
provided all the rctuators are of the same type. We con- 
fined ourselves to some special regions of locations 
in which the behavior of is simple. A ranking system 





was set up for the locations |V) • Also, the relationship 
of the strength of the actuator to its location was indi- 
cated which would give a higher degree of controllability* 
We extended the analysis to include systems which have 
multiple actuators of both force and torque types. We 
showed these could be analyzed independently for each type, 
then adjustments in the locations of actuators made in case 
of conflicts. 

As an example of multiple actuator systems we 
considered a two actuator system. This gave the same re- 
sults as those for a single actuator located on each appen- 
dage. 

Thus, essentially the analysis of a single actuator 
case provides valuable information on the quality of loca- 
tions from the controllability point of view, and this 
information can be used for multiple actuator systems to 
obtain in a rational way some reasonably good distribution 
patterns for the actuators. 



6. SUMMARY AMD CONCLUSIONS 
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We started with the problem of attitude and shape 
control of very large flexible spacecraft, in the last four 
years a large nundber of potential future spacecraft projects 
have been identified which require spacecraft of wiprece- 
dented site, and hence unprecedented flexibility The 
problem is no longer confined to the attitude control of a 
rigid spacecraft, but both attitude and shape control for 
the entire vehicle. In order to achieve attitude and shape 
control actuators have to be distributed over the entire 
vehicle. How should the locations of the actuators be 
chosen in order to best control the flexible vehicle? This 
requires significant advances in the state of the art, and 
little has appeared in the literature toward this direc- 
tion. 

In this work we started by showing the necessity of 
a concept of the degree of controllability of a system so 
that it is possible to cong>are different distribution 
patterns of actuators based on this measure. Since the 
spacecraft dynamic equations can be put into the linear 
time invariant state space form after linearisation about 
a nominal state, we considered general linear s>’s terns for 
analysis. Then we conducted a search for a meaningful 
definition of this concept of degree of controllability. 
Several candidate definitions were scrutinised and found 
unsuitable, because they did not include all the pertinent 
factors such as controllability, total time, control 


\'*'T , 


3 > 




effoct, stability, and control objective, which will have a 
bearing on the degree of controllability. Hence, a defini- 
tion bad 1.0 be sought from fundamental physical considera- 
tions. The approach led to a meaningful definition for the 
deg.'ee of controllability. 

Having formulated this definition we sought for a good 
approximation so that this is applicable to real problems 
and numerically manageable, the approximation approaching 
the true value as the computational effort is increased. 

The mathematical approach adopted, besides leading us tu 
our goal,^ showed the system equations in a special form 
thereby enabling us to derive some relatively simple tests 
for complete controlledsility of the system. Thus, we 
developed an algorithm for computing the degree of con- 
trollability of a system. 

To apply these concepts to our problem we had to 
first obtain the state space form for the dynamic equations 
of spacecraft motion. We analyzed a typical model of a 
large flexible spacecraft and derived all relevant equations. 
Then, to carry out a numerical example a specific model of 
a large flexible spacecraft was examined. The effect of 
actuator locations on the approximate degree of control- 
lability of the system was analyzed for single and multiple 
actuator distributions* The single actuator case, for 
both force type and torque type, yielded valuable infor- 
mation on the quality of locations in the system of 
actuators with respect to the degree cf controllability. 

We identified the tip of the appendages as the best 
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location for a single actuator of a force type or torque 
type. This is for the model described by t»ro appendage 
modes. This conclusion is very 91 ^<ral in ^he sense that 
it is independent of the tine T and the weighting factors 
ar.sociatea with the modes. The weighting factors reflect 
the relative inportance in controlling the different inodes. 
There wer other salient features of the single actuator 
analysis. These include: 1 ) the regions of poor control- 

lability near the zeroes of the influence cur'/es which 
show the effect of a single actuitor location on the degree 
of controllability, and 2) the regions in which the degree 
of controllability either steadily improves or never 
deteriorates (for torque actuator case) as the location of 
the actuator is changed toward the Mp where it reaches 
its maximum. 

For the multiple accuator case, it was shown that the 
single actuator analysis could be applied for each actuator 
separately, using the features obtained for tlie single 
actuator case, and making adjustments in case of conflicts 
for the same locations. Thus, the information obtained 
from the single actuator analysis is applicable to the mul- 
tiple actuator case, and can be used to set up preliminary 
distribution patterns for a chosen number and type of 
actuators. Then the degree of controllability can be 
evaluated, if necessary, in case some distributi'.as Involve 
locations for which information is not clear from the 
single actuator analysis (for thee, the weighting factors 
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and tima T say b« iafiortant) • Onea thia aaaauira of control- 
lability it avaluatad it can be deeidad aa to whicAi die- 
tributiun ia of the hi^heat quality (neat deairable) from 
the controllability point of view. 

In thia work the following can be conaidered aa the 
riin contributiona t 

1) The concept of the degree of controllability— 
a neaningfal defir ition for all linear time invariant 
ayataiBs (actually the definition is applicable to 
nonline 11 yeteiua also) • 

2} i n a goritha for a reaaonable approximation for 
the degree of controllability based on an approxi- 
mate bound for the recovery region. This approxi- 
mation satisfies the property that it is zero if fud 
only if the true degree of controllability is zero. 

lia algorithm is relative 'y single from the numerical 
point of view. Jaing this approximate measure various 
distributions of actuators can be ranked in descend- 
ing order of their desirability in a practical appli- 
cation. 

3} A logiccl and a very useful approach for pre- 
liminary distribution pettems for the actuators so 
that a designer need not ;lndulge in a blind search 
for some distributions to be studied. 

4) Relatively sii;q>le controllability testi for all 
linear time invariant systems. The minimum nundaer 
of actuators required for complete controllability of 
a system is a by-product of these t&sts. 



As a final ranark, ainca tha analysis for tha dagraa 
of controllability is basad on linaar tiaia invariant systans 
it can ba appliad to problans othar than tha spacacraft 
actuator distributions. This will involve identifying a 
paraaatar# analogous to actuator location, which will affect 
the degree of controllability of a given systasi. Of course, 
the equilibrium state must also be identified. 
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Table 5. 

Locations of Zeroes of Influence Curves for Force 
Actuators 


Mode of Motion of Bp 


Curve Symmetric 

Location Ir I 


0.0 

O.OA77 

0.0376 0.8166 



Ant isymmetr ic 


Locat ion 


0.255 

0.05A5 




0.815 


Table 6. 

Locations of Zeroes of Influence Curves for Torque 
Actuators 


Mode of Motion of Bp 

Symmetr i c 

Ant isymmetr Ic 

Location Ir ( 

Location Ir 1 
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8. AEraroicBS 

A. Matrix Algebra 

(i) By repeated postnultiplication of P'^APO in (2.4-1) by P*'a 
on the left side and JP"* on the right side it follows directly that 

•4c •••• 

(ii) Itie exponential of any matrix -Rt can be written as 



E - Rt + 


tl 


'•t 

7T 

e 


Suppose then premultiply both sides by P 


to obtain 


(A-2) 



2l ^ 

ll 31 



(A -3) 


v4iere (A-1) has been used for each term on the right side to obtain the 

series in J in the paranthosis (this Is obtained from (A-2) 
bv substituting R*J). 

(iil) From (2.4-1) we can inned lately write 
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where J >C , the unit matrix. Substitution of this on the right side 
of (A-2) for R>0 for all powers of e( , leads to 
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-r.t -T,t 
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where are the Jordan blocks of order . 

(iv) The matrix can be written as 

A 

where E|^ is the unit matrix of order , and is a nilpotent 

matrix of index as shown below. 


A 


O I O 

0 O I 
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o o 

o o 
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o 0 o O 


Because ^ £, and N, are commutative e"*^' can be written as 
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Expanding e 




by using (A>2) (R*\^) we obtain the result 
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Substituting this in (A-7 ) , we have 


-r„t -A^t -NkA: 

c « -e. c 


(A-1^ 


Again, by substituting R» N|^ in (A-2) and using (A-6) we get 


» J-l i^ir' 

w V- ^1- K X t_X ^ 3 i,i 
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K-»;! 

(A- »e) 


The powers of are easily evaluated by shifting the superdiagonal in 
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(A-5). Hence, by expansion of ^ 
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Nw, 


...,t)^^-l, in (A-IO), and adding up all the terms we obtain for e 

the following form: 
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ABSTRACT 


The unsolved problem of how to control the attitude and shape of future 
very large flexible satellite structures represents a challenging problem 
for modem control theory. One aspect of this problem is the question of 
how to choose the number and locations throughout the spacecraft of the 
control system actuators. Starting from basic physical considerations, 
this paper develops a concept of the degree of controllability of a con'> 
trol system, and then develops numerical methods to generate approximate 
values of the degree of controllability for any spacecraft. Tnese results 
offer the control system designer a tool which allows him to rank the 
effectiveness of alternative actuator distributions, and hence to choose 
the actuator locations on a rational basis. The degree of controllability 
is shown to take a particularly simple form when the satellite dynamics 
equations are in modal form. Examples are provided to Illustrate the use 
of the concept on a simple flexible spacecraft. 

INTRODUCTION 

In the last few years a large number of potential future satellite pro- 
jects have been Identified which require spacecraft of unprecedented size, 
and hence unprecedented flexibility. The attitude control problem for 
such a spacecraft is best characterized as simultaneous pointing control 
and shape control of the vehicle. In order to achieve shape control, or 
equivalently control of the various modes of oscillation of the flexible 
structure, it will be necessary to distribute actuators throughout the 

>This research was supported by NASA Contract NAS 8-32212 with the 
NASA Marshall Space Flight Center. 
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vehicle. How should Che number and locations of actuators be choseu in 
order to best control the flexible spacecraft? 

This problem has been recognized for some time, but to date little has 
appeared in the literature that would help guide Che control system de- 
signer in placing the actuators. Most of the known results identify the 
minimum number of actuators needed for a given set of modes to be con- 
trolled, and identify certain specific actuator locations which cannot be 
used because they result in an uncontrollable system. 

The concept of controllability in modern control theory is a binary 
concept, either a system is controllable or it is uncontrollable. Start- 
ing from a set of actuator locations which produce an uncontrollable sys- 
tem, but for which the number of actuators is sufficient to produce con- 
trollability, it will usually be the case that moving one of the actuators 
by a distance e > 0 can produce a controllable system, no matter how small 
the e . One expects that for a small e , even though technically the 
system is controllable, in some sense it will not be very controllable. 

A precise definition of this concept would prove useful for actuator 
placement. 

It is the purpose of this paper to generate, starting from basic physi- 
cal considerations, a definition of the degree of controllability. The 
definition obtained is certainly not the only possible definition, but it 
does have the advantage over a definition based on singular value decompo- 
sition that the physical reality of actuator saturation limitations is 
included in a fundamental way, and that time limitations on accomplishing 
one's control objective can be included. 

The definition is then applied to the actuator placement problem for 
flexible spacecraft. With this tool the control system designe;: can rank 
the desirability of various candidate actuator distributions, and thus he 
would have a rational way of picking which distribution to use. 

DEFINITION OF THE DEGREE OF CONTROLLABILITY 


Let us consider any general linear time invariant system in state var- 
iable form 


i*(t) - Ax*(t) + 3u*(t) (1) 

where x* « k and u* < r"*. It should be noted that although we focus our 
attention on this system, the degree of controllability definition which 
we adopt is also applicable to more general systems of the form 
i**(t)«f (x*,u*,t) having a solution z*(t)50 (f (0,0,t)*0) . 

It is instructive to discuss some of the candidate definitions of the 
degree of controllability which were considered and discarded— the process 
of starting with a blind attempt at a definition and progressing to a well 
formulated concept highlights the characteristics that a workable defini- 
tion must have. It is tempting to try to connect the degree of controlla- 
bility to properties of the standard controllability matrix 
Q ■ [B AB • • • and define degree of controllability as the square 

root of the minimum eigenvalue of QQ*^. Four apparent difficulties with 
this definition must somehow be handled before the definition becomes 
viable. These difficulties are as follows: 1) The degree of controlla- 

bility is affected by a transformation of coordinates (since the eigen- 
values of are not Invariant under changes in state variable representa- 
tion). 2) This candidate definition satisfies the basic requirement that 
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Che degree of concrollabillCy is zero when Che sysCen is unconcrollable, 
buc ic is noC innediacely clear vhac ocher physical meaning can be ac- 
Cached co che size of Che eigenvalues of QQ^. 3) The candidate definl- 
cion does noc involve a dependence on the amount of time T allotted to 
accomplish che concrol cask. It can be much easier to control the system 
state in some directions in che state space at one time than at another 
time, so the degree of controllability should depend on T. 4) It is not 
clear chat che amount of control effort needed to accomplish the control 
Cask is reflected in this definition. In the satellite described above 
where one actuator has been moved by a small amount c to produce control- 
lability, one expects Che "weak controllability" of Che system to be mani- 
fested in the need for very large control actions to accomplish certain 
small changes in the state, and hence Che control effort required is of 
fundamental Importance in making a definition. 

It is clear chat some type of limitation or standardization of che con- 
crol effort must be included in the definition. Consider a standardiza- 
tion which restricts the concrol to a unit Impulse, and consider systems 
with A in diagonal form and with u* a scalar. For distinct eigenvalues 
the system is controllable if none of the elements b{. of the column matrix 
B are zero. Furthermore, these components indicate how far a unit impulse 
control will move each state component instantaneously, so one might sug- 
gest che min|2?t) aa a degree of controllability. Here we are trying to 
generalize ^a second standard test for controllability Co obtain a degree 
of controllability definition. Among the apparent difficulties with this 
candidate definition is the fact that the concrol actions are so re- 
stricted that che components of the state cannot be affected indepen- 
dently. The control of all states by a single concrol u* relies on Che 
differences in che djmamlc behaviors of Che states. 

Both of these candidate definitions have difficulties; they do not 
appear to Include the effects of all pertinent variables. Hence, it will 
be necessary to build che definition from more fundaments] considerations. 
Ironically, when this is completed and interpreted properly, in certain 
special cases Che degree of controllability definition will essentially 
reduce to che second candidate definition above (and by employing a diffe- 
rent approach involving singular value decomposition of matrices something 
of Che general form of che first candidate definition can result). 

It is now evident that the definition of che degree of controllability, 
besides being in some sense a measure of how easy it is for Che controller 
to control che system, must in some way handle four things: 

1) It must have che property that che degree of controllability is 
zero when the system is uncontrollable. 

2) It must somehow consider dependence on t^tal time T. 

3) It must standardize Che concrol effort in some way. 

4) The concrol objective must be restricted. 

Concerning che last point, certainly different control objectives should 
influence che choice of the control system design, and hence the degree of 
controllability of a candidate design should be keyed to che objective 
involved. In a large class of problems (regulator problems), the equili- 
brium solution to equation (1) is of primary importance, and the 
control objective is to return x* to zero after a disturbance. Since this 
is che most common attitude and shape concrol problem for flexible space- 
craft, we will restrict ourselves to this objective. A companion 
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paper [l] develops Che concept of Che degree of concrollablllcy for var- 
ious saCelllce slew maneuver objectives. Concerning Che standardization 
of the control effort we will require that Che control components satisfy 
|u^|£ 1 for t*l,2,3,. . . ,m, which represents realistic physical limita- 
tions of the actuator capabilities. Note that the use of one as the 
bound for all control components Implies normalizing each component of u* 
to produce a new control vector u, and adjusting the B matrix to produce 
a new matrix B. 

Controllability requires the existence of a control function which can 
transfer any initial state to any final state in finite time. With our 
more limited control objective, the degree of controllability should be 
related to the volume of initial system states (or states resulting from 
disturbances) which can be returned to the desired state x**0 in time T 
using the bounded controls. Consider the nature of this volume in more 
detail. In an uncontrollable system there will be at least one direction 
in the state space for which initial conditions in this direction cannot 
be returned to the origin, and the volume will lose one or more dimen- 
sions. For a controllable system whose parameters are such that it is 
nearly uncontrollable, only Initial conditions very close to x**0 along 
the above mentioned direction could be returned to the origin in time T 
using the bounded controls. Hence, we will generate a definition of the 
degree of controllability based on the minimum distance from the origin to 
a normalized state that cannot be brought to the origin in time T. More 
loosely it is the minimum disturbance from which the system cannot recover 
in time T. 

The coordinates of a state space will very rarely all have the same 
physical units, and hence it is clear that comparing distances in the 
state space will require that each coordinate must be made unitless by 
normalization. How should one choose the normalization to use? Recognize 
that when comparing two controller designs for controlling the same dyna- 
mic system, the needed minimum distance for each design will usually cor- 
respond to a different direction in state space. Hence, ranking of the 
degree of controllability of the two systems will depend on comparison of 
distances in different directions, and this Implies that we must be 
equally interested in controlling deviations of the state from x*-0 in all 
directions in the state space. In order to accomplish this the control 
system designer must specify n-1 numbers which represent his degree of 
interest in controlling each component of the state. This could be done, 
for example, by determining the deviations of x^* , Xj* ,...,x*_, which 
would be considered of equal Importance to a deviation of Xq*> 1. The 
reciprocals of these numbers would then be used to produce normalization 
factors for each of the coordinates of the state space giving a new state 
vector X. The system equations expressed in terms of the normalized state 
X and normalized control u are then written as 

x(t) ■ Ar(t) + Bu(t) (2) 

I u^ I $ 1 r*l , 2 , 3 , . . . ,m 

Just as in optimal control theory where the control system designer 
must be specific about his goal by specifying a cost functional, in order 
to define the degree of controllability, it is necessary to be fully spe- 
cific not only about the objective of keeping x*0, but also about the rel- 
ative importance of keeping each component of x near zero. 

Relative to the normalized system (2) we are now ready to make the 
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following definitions: 

Definition 1 : The recovery region for time T for nomalized system 

(2) is Che sec 

«- |x(0)|3u(t). t€[0,T], |u.(t)|£ 1 for 
i-l,2,...,m » x(T)»0[ 

Definition 2 i The degree of controllability in time T of the x»0 solu- 
tion of normalized system (2) is defined as 
p - infl|x(0)|j V x(0) i R 
where || * || represents the Euclidian norm. 

Thus, the recovery region identifies all of the initial conditions (or 
disturbed states) which can be returned to the origin in time I using the 
bounded controls. And the degree of controllability is a scalar measure 
of the size of the region, where Che scalar is chosen as Che shortest dis- 
tance from Che origin to an initial state which cannot be returned to the 
origin in time T. 

The degree of controllability, as defined, is keyed to the state vector 
X employed. No transformations of coordinates can be allowed once the 
normalization has been specified (unless Che norm used in the definition 
is adjusted to compensate for the resulting distortion of the state space). 

Note the following property of Che recovery region: 

Remark : The recovery region R for system (2) is Che same as the sec of 

reachable states for time T for the system 
x(t) - -Ax(t)-Bu(t) t « [0,T] 

|u. |< 1 t*l,2,3, . . . ,m 
starting from xl0)«0. 

It should be pointed out Chat although Definition 2 incorporates all 
Che properties which were Identified as necessary in the definition of the 
degree of controllability, it is not necessarily unique in doing so. For 
example, a standardization of the control effort in terms of energy can 
also be employed, but the inequality saturation constraints on the con- 
trols used here represents Che more rea.lstic situation. 

An example section is provided in this paper in order to illustrate how 
the degree of controllability behaves as a function of actuator placement 
in a simple flexible spacecraft, and to demonstrate that the definition 
behaves according to our limited Intuitive notion of the degree of con- 
trollability. 


CONCEPTS FOR APPROXIMATING THE RECOVERY REGION 

In order to make the definition of Che degree of controllability use- 
ful, it is necessary to develop a simple algorithm to generate at least 
an approximation to the distance p. This necessitates approximating the 
recovery region R . 

Note chat 

x(T) - ♦(T,0)x(0) + ♦(T,0)4*(0,t)Bu(t)dt (3) 

where 0 is the state transition matrix for (2). The distance moved during 
time T is ^^(O)-x(T), and we are concerned with sending the system to the 
origin so that x(T)«0. Then the initial state C which reaches the origin 
in time T using control u(c) is elven by 

<•1 ♦(0,t)Bu(t)dt (4) 

By Che Caley-Hamilton theorem the state transition matrix can be written as 

♦(0,t) - e-^^ - <|»jj,(t)A“ (5) 
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V^Z. 


where the are scalar functions of time* Partition the B matrix Into 
column matrices and define the following matrices 

1.7 


B • [ 2?) ^2 2>s ••• ] 


■^sl 


« - I 

'!»” • [ '!'» 'I'j 'I'l ^ 

Q - [ B AB A^B ... A""^B ] 

Qg • [ 2?g Abg A*bg ... A 
Then 5 can be represented in the following alternative forms 

■“££</ 

m 

-V 

m 


t|>^jj(t)ug(t)dt I 


A% 


6 


. 


•*B**'"’8* 


n-1 


igjugdt 


< 6 ) 

r' 

't') 

(9) 


(10) 


( 11 ) 


( 12 ) 


For the purposes of illustrating certain concepts, let us restrict our> 
selves to the case of a scalar control so that the sunmatlons over 6 as 
well as the B subscripts in the above can be dropped, and B is a column 
matrix b. Also let n«2 for simplicity. Suppose the recovery region is as 
shown by Region I in Fig. 1. The maximum Xj component cf any state in the 
recovery region is obtained by using the control u equal to minus the 
signum function of the first component of the vector [Q<{>] in (12) , since 
this maximizes the x^ component of the Integrand at each time t. The 
right hand side (and left hand side) of the rectangle enclosing this re- 
covery region in Fig. 1 can thus be found by integrating the first compo- 
nent in (12) using this control. If desired the point at which the recov- 
ery region touches this side is obtained by integrating the second compo- 
nent of (12) using this control. The top and bottom of the rectangle are 
found similarly. 

The rectangle obtained in this manner might be considered an approxima- 
tion to the recovery region, and then the shortest distance from the origin 
to one of the sides might be considered an approximation, p, to the degree 
of controllability, Note that this necessarily produces a p which 

is an upper bound for the degree of controllability. In some cases This 
approximation is a tight one, but often it is not. Suppose the recoveiiry 
region was Region II of Fig. 1. This corresponds to a system which has a 
much poorer degree of controllability, t yet the approximation 23 

remains ttte same. In fact, suppose that P»*«^ in such a way that Region II 
degenerates to a line forming a diagonal ot the rectangle in FJ^. 1. Then 
the system is an uncontrollable system, but the approximation p still pre- 
dicts a good degree of controllability. Hence, this approximation must be 
rejected. 

For the case of a scalar control being considered, this shortcoming can 
be eliminated by using as expressed in (10) and maximizing components 
along A^b. The control 

u(t) - -sgn(i|>g(t)] (13) 

cxtremlzes the coefficient of the vector A 2> in (3).^ will simultaneous- 
ly produce some components a.'*'*ng the other vectors A b for This is 

a maximization of a component the vector 5* but it is a component as seen 
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in a nonorthogonal set of coordinates. Hence, the upper bounds obtained In 
the various di^rectlons define a parallelogram (more generally an n dimensional 
parallelepiped) which can be considered as an approximation to the recov- 
ery region, as shown In Fig. 2. As before there Is some point on each 
side of the parallelogram which Is In the recovery region, but no point 
outside the parallelogram is in the region. 

The minimum distance to a side of the parallelogram, l.e. the minimum 
perpendicular distance to a side, is an approximation p* to the degree of 
controllability, p. Vfhen the system becomes uncontrollable, the columns 
of Q become linearly dependent, and hence the perpendicular distance to 
one of the sides becomes zero. This means that this p* has the essential 
property that p*«0 whenever p«0. 

We conclude that for the scalar control case we have a viable method 
of approrlmatlng the degree of controllability. A simple method %rlll be 
presented in a later section to determine the needed minimum perpendicular 
distance. 

This approximation is still an upper bound, and it can be Improved, in 
fact made arbitrarily good, by considering more directions In the state 
space. Let e be any desired unit vector expressed as a column matrix of 
components. By examining (12) the state in the recovery region having 
a maximum component along the dlreetlon_e is obtained using the control 

u ■ -sgn(e^Qi|>] (14) 

and hence no points In the recovery region lie beyond the line perpendicu- 
lar to e and a distance ^ 

I |e^#|dt (15) 

from the origin (but at least one point in the recovery region lies on 
the line). Figure 3 Illustrates how use of three e’s (e^, e_ , and e ^ ) 

Identifies three tangents to the recovery region, and when t^en together 

they begin to approximate the region boundary. Let 0 Le the minimum 
value of 0^ for any set of directions e considered. Then an improved 
estimate of the degree of controllability is p**>mln(p*,0) , and p**>p can 
be made arbitrarily close to the true degree of controllability p by 
picking a sufficient number of directions e. This method of improving the 
approximation to the degree of controllability will also be generalized to 
the multiple control case. 

FUNDAMENTAL EQUATION FOR RECOVERY REGION APPROXIMATION 
IN THE MULTIPLE CONTROL CASE 

The previous section presented a procedure for generating an approxi- 
mation p** to the degree of controllability p in the case of a scalar 
control u. The procedure required the use of n carefully chosen direc- 
tions in the state space, 2 ?,a 2>, . . . ,A"~^i>, in the aoproxlmatlon to the 
recovery region in order to insure that p** had the property that p**«0 if 
and only if the system is uncontrollable. If the control vector is m dim- 
ensional with ir.>l it is no longer obvious how to obtain this property, 
since the columns of the Q matrix necessarily contain linearly dependent 
vectors. Instead, we will consider the eigenvectors and generalized 
eigenvectors of the A matrix of (2) as the n linearly independent direc- 
tions in the state space. Certainly some modifications must apply when 
these vectors are complex. In the single control case the value of p*** 
became zero when the system became uncontrollable because linear depen- 
dence of the vectors i,Ab, . . . ,A”"^i implies the collapse of at least one 


> 7 
I 
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dimension of the psrallelopiped. The vector:: chosen here Tor the 
multidimensional control case do not exhibit t:his reduction. Neverthe- 
less, It will be shown In the next section chat the desired property of 
the resulting p** can be demonstrated under fairly general assumptions. 

This section is devoted to generating the appropriate expression for ^ 
equivalent to equations (4,10-12), expressed in terms of components in 
these eigenvector directions. 

Let J be the Jordan canonical form of the matrix A, and let P be the 
matrix of eigenvectors and generalized eigenvectors so that 

P-1aP - J (16) 

J * dlag(J. , J2, . . . , Jy) 

where the are the square Jordan blocks of dimension Ux* Associated 
with each block is in eigenvalue f»l,2,3, . . . ,r, so tnat rSn and r is 
greater chan or equal to the number of distinct eigenvalues. Also, let 
pj be the n columns of P, and be the n rows of P”^ (the left eigen~ 
vectors and generalized left eigenvectors). 

The desired fundamental equation for ^ is given in the following theo- 
rem. The theorem is made significantly more complicated in order to 
handle repeated roots, but it is necessary to treat such roots since the 
rigid body modes for any spacecraft involve double roots. 

Theorem 1 : The displacement C “ x(0)-x(T) after time T of the svstem 

state for equation (2) resulting from control u(t), c < [O.X], 
can be written as « _ 

0 

where Ij is a possibly time dependent vector 
ij(t) ■ V'Ai 1'j2 ••• 


"J'"' V2 ••• > 


(fei-j)! 




jb - <7- Og uu; 

with given in (6). The values of are determined by 
the dimensions of the Jordan blocks as follows 

feo*0 ; ^ 

and the values of j associated with each t»l,2,...,r are 

J " ^'£—1^^* ^£—1^^’ *** ’ ^£ (22) 

The h. are given by " . 

^ h. • e“^£® 

for £ and j related by (22). 

For the special case of all distinct eigenvalues (r*n) equation (17) 
simplifies significantly, since there is no need for a distinction between 
£ and J, and the components of l> become the constants I 
Then (17) becomes ^ ^ j6 JB J B 

-C - / ^ (24) 

Due to space limitations the ^-roof of this theorem Is omitted, but 
can be found in reference [2] or the Journal version of this paper (to 
appear) . 
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FUNDAMENTAL EQUATION IN TERMS OF REAL VALUED FUNCTIONS 

Take one eigenvalue from e»ch complex conjugate pair and asaign the j 
or j*s associated with its eigenvector or generalized eigenvectors, p>, 
to the sut C, Let C* represent the set of j's associated with real eigen- 
values. 

For any complex eigenvalue the associated p., ft., and 1. in (17) vill 
be complex valued. Let their real and imaginary parts be given as follows 
together with the definitions of 6. and n* 

p . i rMsi 
V J J 

Zj ■ y ^ ... ,n (25) 

• c .Y - p .6 

'’j “ 

The vectors representing the semiaxes of the parallelepiped are given 

by 





3 *C 

(26) 



(27) 


3 t C* 

(28) 


' sit ■ JP' ' J 

(Note that these results can now be applied without regard for the reor- 
dering of the eigenvectors used in their derivation.) 

Let this set of n vectors be denoted by Vj, Vj, ... , v^. The minimum 
perpendicular distance to a side of this paralleloplped can be used as an 
approximation to the degree of controllability, ond as before the approxi- 
mation can be made arbitrarily tight by using distances in additional di- 
rections in the state space. 

Define the matrix F»[vi ... v ], Each vector goes from the ori- 
gin, or center of the parallelopipeS volume, to the center of one of its 
sides. Any surface of the parallelepiped consists of edges that are para- 
llel to n-1 of the v vectors. Let d., corresponding to v., represent the 
perpendicular distance from the origin to that surface for which no edge 
is parallel to ^•. Then d> is the component of v. normal to this surface. 
Theorem 2 ; Tne normal distances d., j«l,2, .A ,n, to the surfaces of 
the n dimensional parailelopiped prescribed by F are the 
reciprocals of the magnitudes of the column vectors given 
by (f’T)-I. 

Proof omitted due to space limitations (see [2] or journal version). 

The approximation p* of the degree of controllability p from the para- 
lleloplped boimd on the recovery region is 

p* ■ min d , (29) 

j ^ 

This approximate degree of controllability can be made arbitrarily tight 
by including the recovery region bounds in other directions as well as 
those of paralleloplped axes. This approximation p* will go to zero when 
the system becomes uncontrollable as indicated in the following theorem. 
Theorem 3; Suppose that all eigenvalues of the matrix A for system (2) 
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are diatlnct. Then p*«0 if and only if p*0^ That Is, 
the approximate degree of controllability based on the 
mlnimua perpendicular distance to the surface of the para- 
llelogram F obtained from (26-28) will be zero if and only 
if the system is uncontrollable. 

Proof omitted due to space limitations (see [2] or Journal version). 

SPECIALIZATION OF THE DEGREE OF CONTROLLABILITY CRITERION 
TO H(X)AL COORDINATES 

Consider a lightly damped flexible spacecraft with dynamic equations 
expressed in terms of spacecraft normal modes. Temporarily consider the 
shape control problem alone so that the rigid body mode is neglected. 

Then the equations are 

fi. + w.V. • r.^u t-l,2,...,(n/2) (30) 

The r. and u are defined so that each component of u is bounded by unity. 
Let tne numbers N. represent the relative importance we assign to the 
control of the He must also specify the importance we assign to 

Che control of in the state space trtilch will be generated. Since Che 
system Is lightlv damped, if is sln(b)^.C-t^), then will be u^.cos(a).c4^) 
approximately, so chat the relative importance of controlling /\.^ls repre- 
sented by Now gsnerate the state variables as follows ^ 


*2i-l " ’ *2i ■ *2i-l'“'t 

so Chat all unit deviations from the origin of Che state space will be 
considered equally serious. Then 

’H-f° ^ “1“ 

fli J hi ■^'i“i|L*2J 

Let the coefficient matrices be A. and B^. Then the normalized system 
equation (2) has coefficient matrices A and B, and eigenvector matrix P, 
which can be partitioned as follows 


^2i-l^“t 


A • dlflg [ f A 2 ) 

B - [B^"^ Bj^ . . . 
P - diag [P^, Pj, 

P"^ - diag [Pj^*^, 


^n/2^ 




where 


t • r ^ ^ 1 ; 1 r ] 

L^2i-l^“t ^2i“^2t-l L’^2i-1 “i J 




^2t • 


\ ' 
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T 

Using the rows of and the columns of B in ^ 


real part (Y2^ imaginary part 

,th 


6 * 


gives a zero 


^2i.S " 


where F^g is the 6 component of F^.. We wish to calculate (26) and (27) 
for j *Ct i.e. for one of the complex roots in each complex conjugate pair. 
From (25) g * ~^2i^2i B* ***** ^2t imaginary part of exp[-X2^t]. 

Tlien (26) becomes 


( f 

\ i ^ / 


(33) 


For a lightly damped system the exponential will not change significantly 
during one oscillation of the sine wave. Assuming that T is long compared 
with the period of the sine wave (so that the effect of partial completion 
of the final period of oscillation of the sine wave is negligible) • the 
absolute value of the sine can be replaced by 2/x, its average over a 
period. Then (33) becomes 


2 i). 


(34) 


II r^il 




(35) 


A ^ 

The analogous calculation for (27) gives the same result with r» . replaced 

hich are the 


V VI 
n 


by s,^. Together these vectors form the set v., 
columns of F. 1 Z 

fhe snd 82^ sre the real and imaginary parts of the appropriate 
column 01 P. LooKing only at the appropriate partition of F, call it F^. , 
we have 




1 




where a is the coefficient of 
from Theorem 2 result from tal 
given by 




— 1—1''^ 0 
oA:« 7 (_ -c -1 

in (34). The associated va 


F. 

z 


g the magnitude of the rows of F. 


ues.of d . 
an<T 


are 


°2i-l » 
the second of which is the smaller, 
controllability 0* is 





*2i 

From (29) the approximate degree of 


1+5^ 


P* 


■ .Inf^ II ■'illA - 1 \1 
^ Vi 

Ir X % 


(36) 


Note that when damping is present, the recovery region expands with T in 
such a way that p* grows exponentially. Furthermore, when the natural 
frequency of the minimizing mode is decreased the approximate degree of 
controllability increases. 

An important special case is that of a system model which has no 
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damping. The approximate degree of controllability in this case becomes 


2T 


P -T 


* ■ 


rain 

i 




(37) 


Note that p* Increases linearly with time, making T a scale factor which 
can be dropped. We conclude that the approximate degree of controllabi- 
lity is simply the minimum norm (in the sense of (35)) of the rows of the 
suitably normalized control coefficient matrix when the system is repre- 
sented in modal form. The normalization in (37) is that applying to the 
derivative of the modal coordinate hv* 

Now let us introduce the rigid body mode into the problem so that the 
control objective is simultaneous attitude and shape control of the 
flexible vehicle. For simplicity consider only one such mode, and let the 
variable Involved be 6. Let the normalization representing our degree of 
Interest in controlling 6, be N^^and that for 6 be, Ng. Then define state 

becomes 


variables s^^d/Ng 


and xs-6/Ns 


e 


0^ bliob for e be Nt . 
°and the equation O*!*^** 


ft] • [: V*] ft] . [,;„] 


and 


From (20) and 
1. Then (28) 


r “ ]• 

LO Ng/NgJ' 




(19), 

gives 


‘06 


• 0 , 


II rJI 


06 

2 


’^06^**0’ 


|-1 0 1 


'66 ‘ 06 ’ 


and h 




9"A 


N 


0' 




(38) 


Calculating d^ and dg using Theorem 2 results in l:he coefficients of these 
two vectors. 

The approximate degree of controllability p* when the rigid body mode 
is included is the minimum of the p* given by (37) (or (36) when damping 
is Included) and the two coefficients of the vectors in (3b): 


p* ■ min 1 ( 


II r,||* 

..ii'bIIa 1 

(39) 

L\ 

, IT / i N 

2 / Ng 

”0 



Note chat because of the modal representation used in the system equations, 
there is no coupling of the weighting factors N., N., and Ng from one mode 
to another. Thus, if one is not particularly interested in'^controlllng 
any specific mode, the corresponding N can be made very small. Then that 
particular mode is simply neglected in the minimization to determine p*. 

For example, in some applications one might feel that it is important to 
control the rigid body angle 0, but that if 0 is controlled well the value 
of 0 is unimportant. In such a case one simply Ignores the 
term in (.■’9). o a o 
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EXAMPLES 

^ The methods of this paper were applied to an example flexible space- 

craft. Due to space limitations only a few results will be cited. The 
spacecraft consists of a rigid symmetric central hub with two 60 foot 
radial stem typ« booms. The two lowest appendage modes are used to gene- 
rate two symmetric and two antisynnetrlc spacecraft modes. Both torque 
and force actuators were considered. Figure 4 applies to antisymmetric 

* spacecraft modes and plots the values of normalized by the 

torque actuittor strength. Equation (37) shows that the degree of control- 
lability is calculated by scaling each curve of Fig. 4 according to the 
normalization N^and the actuator strength, and caking Che minimum of these 
curves at the actuator location. The optimum location for the torque 
actuator is easily determined as that point which maximizes this minimum. 

* Any regions where one of the curves is near zero should be avoided. For 

a system which truly has only two modes, and with unit normalizations, Che 
best Corquer location is at Che end of the appendage. 
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ACTUATOR LOCATION 

Fig. 4. Decernlnatlon of the degree of controllability as a function 
of actuator location 



INSTRUMENT FLEXIBILITY LIMITATIONS FOR 


STABILITY OF MULTI-PURPOSE SPACECRAFT 
WITH FIXED CONTROL SYSTEM DESIGN 


9. AMVUSIS 


9.1 Introduction 

In part IWo, attention is focused on a component of a spacecraft 
while keeping the control system design fixed. Ihese ^cecraft are so 
designed that certain of their components are interchangeable with other 
physically compatible components. These interchangeable components may 
be designed to accomplish different tasks. These spacecraft can, 
therefore, be called multi-purpose spacecraft. 

A typical example is a Shuttle-based Instrunent Pointing System 
(IPS) . Such a ^cecraft will have a Space-Shuttle with a multi-purpose 
device called the Instrunent Pointing Houit (IPM) to which one of a 
family of scientific instruments is attached. The instruments are 
typically lightweight and their flexibility becomes significant due to 

pointing requirements. The shuttle and IPM ere modeled as rigid 

✓ 

components of the idealized system, while the instrunent is modeled as a 
flexible component. A control system can be specifically designed for 
such an IPS with a given instrunent. However, the central idea of the 
IPM is to make it a multi-purpose device, so that it works not simply 
for one instrunent but for a whole series of instrunents, many of which 
have yet to be designed. In fact, once the IPM characteristics arc 
defined, the scientific instrunents have to be designed so that they 
will function properly when attached to the 1PM. This is the inverse of 
the usual design problmn, in which the physical ^ant is given and the 
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control system is to be determined. Instead, here, the physical plait 
is variable because of the interchangeable instrunent. The approach 
taken, therefore, is to predesign a comprehensive control system which 
is shuttle-based, and hence ranains unchanged for the spacecraft for a 
whole series of instrunents that might be attached to the IFM. Severe 
accuracy requirements in the orientation of these jhuttle-based 
instrunents, sometimes to within fractions of an arcsecond, may be 
demanded of the control system design. And usually a dual control 
system will need to be adopted with one control system located in the 
shuttle, and the other in the IPii (for fine control) . In this case, the 
IPM can have motion relative to the shuttle. Now, the question is "what 
kind of instrunents are stably controllable with any given control 
system?" In re^nse to this, it is best to describe these instrunents 
in the most general terms possible so as to allow the designer ample 
freedom in their design. 







9.2 Parameter Plane Technique 






This method has its origins about a century ago in the 
idea of investigating the system response chracterist ics [jB^ 
by algebraic aoproach. Inawide variety of control problems 
the designer is interested not only in the stability of the 
system but also in the essential features of the system 
behavior over time. In its original form the approach began 
by treating two coefficients of a characteristic equation 
as variables, and studying how the roots of the equation are 
affected when these two coefficients are changed. By plotting 
the characterist ic curves in the plane of the variable coeffi- 
cients, the method enables adjustment of these coefficients 
so that the roots of the characteristic equation may be set 
at any desired locations. After the curves are plotted, the 
variable coefficients can be adjusted without any calculations. 
Several researchers have since. then extended the method and 
its applications to various problems such as sensitivity 
analysis of linear control systems, circuit synthesis, steady- 
state response analysis, sampled-data linear systems, and to 
other related problems of linear system design. Also, the 
method has been successfully applied to non-linear systems. 

The parameter plane technique for analysis and synthesis of 
linear and nonlinear control systems is amply described in 
Siljak's monograph P]. Once the system characteristic 
equation has been obtained, the parameter plane method enables 
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the designer to evaluate graphically the roots of the equation. 
Hence, he may design the control system in terms of the chosen 
performance criteria; e.g., absolute stability, damping ratio, 
and settling time. Thus, the technique is the mapping of the 
roots of the characteristic equation from the complex plane 
onto what Is known as a parameter plane. The design procedure 
has also been simplified by Siljak who introduced Chebyshev 
functions into the equations, thereby putting them in a 
suitable form for digital computer simulation. 

The parameter plane method requires that the control 
system be described by a character i st ic equation which is 
transformed into the complex domain (s-domain). Two adjustable 
parameters, o< , |J , (which are of interest to the designer) 
are selected, and the characteristic equation (essentially 
its coefficients) is recast in terms of V , p . Suppose the 
characteristic equation (CE) is 


C6 


.2 S • 0 




( 9 . 2 - 1 ) 


then 

"fj * P) ^ J * 0, I, t,.' (9.2-2) 

The points in the s-plane (complex domain) are described best 
suited to this method by the two system characteristic 
quantities 5* .where ^ damping ratio (Oil5l-^l) 

and is the undamped natural frequency. Let 


s»»x + i y*-Jw^+ iw 
Any power of s can be written as 




( 9 . 2 - 3 ) 





s « X . + I Y • 

J ■* 

where X- , Yj are the real and Imaginary parts 
there exist the following relationships 
Z . a - l T w, 2;-, - 

i 

j* X, ^,4. • • ' • 

Z.- X. or Y- 
a I ; X, * - 

Y. * 0 ; Y, - ‘-'a \ 

Substituting (9.2-4) in (9.2-1) 

CE - Z ( X • + i Y; ) - 0 
jsO 

from which we obtain 


(9.2-i-/ 
for which 


(9.2-5) 

( 9 . 2 - 6 ) 


n 

5 f. X- » 0 

• ” J J 

4 

fj Y; - 0 ^9-2-7) 

Now, the coefficients fj may be 1 inear or nonlinear in the 
parameters o< , p depending on the choice of the parameters. 
Consider a linear form for the fj , i.e., let 


f; . aj « + bj P + Cj 

where a. , b; , c- are constants. Then 

Z -fj (s- * + h r* 2; 


( 9 . 2 - 8 ) 


c D 

in wh ich Z . ■ X . or Y; 

i i i 


(9.2-9) 


we define 
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. B = XljXj C-X^-!’'.; (9.2-10) 

J«o <«o 4»0 

and similarly, A^ , B' , C^ for Y* , then from (9.2-9) 

J 


A*< + Bj*+C = 0 

A^«»< + B^p+ C = 0 (9.2-11 ) 

The coefficients A, B, C etc. are determined if "J and co^ are 
chosen. These two equations can then be sohed for o< and ^ 
if their discriminant is nonzero. 

In a similar fashion, two equations nonlinear in o< and ^ 
can be obtained if fj are nonlinear functions in o<, p , 

In this case, thert will be multiple images in the parameter 
plane corresponding to a root in the complex plane defined 


by -5 , u, . 

The characteristic curves in the parameter plane are plotted 
for various values of ^ • Theso curves are called 

“^-curves and w^-curves, because a set of curves can be 
plotted, each curve corresponding to a fixed value of ^ (orcj^) 
and the other quantity (or J ) varying along the curve. 
Similar curves are plotted for real roots r of the Ct. For 
real roots one of the equations (connected with the imaginary 
part Yj ) in (9.2-11) vanishes, and for various values of <T 
we obtain the (T-lines in the parameter plane. Having plotted 
these curves one can read the roots of the CE corresponding 
to any point in the o< - p plane. Generally, these curves 
are shaded to help relate the crossing of these curves in the 



parameter plane to crossing of the axes In the complex plane. 

The details on the shading convention adopted for the parameter 
plane curves can be found in Siljak's monograph [ 93 . 

One of the most important regions in the parameter plane 
is the unstable region which the designer must avoid in 
designing the system. But, though stability is necessary In 
a wide variety of control problems it Is not sufficient, and 
the designer might be interested in other features of the system 
behavior. The parameter plane can be used in several ways as 
one wishes to study the system response. Thus, for instance, 
a designer might want tc keep $>0.2 for all the roots of 
the equation, and a similar restriction on . Depending on 
his requirements he can mark out the regions of interest and 
adjust his parameters t to obtain best performance for 
the system. 

The following are the main steps in the application of 
the parameter plane technique: 

1) Obtain the characteristic equation (CE) for 

n 

the system in the s-domain, i.e., ^ f* s * 0. 

2) Identify two parameters of interest o< , p and 

recast the coefficients of the CE in terms 

of o( , p . 

3) Use the substitution s « X* i Y* and obtain 

two algebraic equations, f- 
— 5*0 



k) Using the recursive relationships from (9.2-5) 
for Xj , Yj in terms of ^ , solve the 

two algebraic equations for o< , p for various 
values of “5 * • 

5) Plot the '5 -curves or t*>|j-curves in the parameter 
plane ( p plane) . 

6) For real roots use the CE with s»CT, and for 
various values of plot the d"-lines in the 

- p plane. 

7) Having plotted these curves mark the unstable 
regions in the parameter plane, and use the 
shading convention. 

In the next chapter we will study the flexibility 
characteristics of the shuttle-based Instrument of a multi- 
purpose spacecraft applying the parameter plane analysis. 



10. ANALYSIS OF A SHUTTLE-BASED INSTRUMENT POINTING 


SYSTEM OF A MULTI-PURPOSE SPACECRAFT 
10.1 Introduction 

The parameter plane analysis has been applied to study 
stability of attitude control of spinning skylab £l 03 , £l 0 » 
and recently Seltzer and Shelton C'O applied the analysis 
to study the rigidity of spacecraft flexible appendages for 
a two-body model, a rigid body and a flexible Instrument, 
with a spacecraft attitude controller of a standard positlon- 
integral-der ivat ive (PID) state feed back type. 

In the present work we w i 1 analyze a three-body system, 
a rigid shuttle, a rigid instrument pointing mount ( IPM) and 
a flexible instrument. There will be two PID controllers, one 
shuttle-based and the other IPM-based. Once the equations are 
derived two cases will be examined: 1) the IPM is locked to 

f 

the shuttle and only the shuttle-based PID controller is in 
operation, thus reducing this system essentially to a two-body 
system, and 2) the IPM is allowed to have motion relative to 
the shuttle with its PID controller also in operation, which 
is a three-body system. The model of the two-body system here 
differs from the model in 023 in Its feed back control law. 



10,2 Planar Model 
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Figure 12 shows a typical planar model of a three-body 
system in which is the rigid shuttle, B,^^. rigid 

IPM and Bp is the flexible instrument whose characteristics 
are of interest to us. The system is inertially at rest and 
the rotational dynamics is assumed to be confined to a single 
plane of motion with control torques applied in this plane. 

Two attitude controllers of the proport ional- Integra i -di ffere- 
ntial (PID) type are assumed, one located in the shuttle (B,^| ) 
and the other in the 1PM The shuttle-bvised controller 

applies torque 7^ to 3^^ and is meant for large attitude 
control, and the IPM-based controller applies an interaction 
torque, T to B. and -7" to B,, , and is meant for finer control 

of the instrument. The inertial attitude of B,^j is denoted by 
(h)j and the relative rotation of Bj^j^ with respect to B^^^ is 
denoted by The flexible appendage is assumed to be an 

elastic beam and is characterized by distributed coordinates, 
and the hybrid coordinate approach [ 5 ] is used. In the 
stability analysis the appendage is characterized by a single 
assumed mode. 



'«»w< » 


10.3 Dynamics 

The kinetic energy of the system is derived in e similar 
manner as In part I, but is more involved due to the relative 
motion of It is given by 

T, - (a), + -fiC®.)* fiCi)] 

+ C,(n)] + ©I 

• 7m 

00.3-1) 

where 

■*■ pr'’* 

= i ^-= ) -^ C t-'+3[ r) h; 







[ S', (Xy -'♦) + SL( <4- xjJ 


f,cn) * i 

■fa®J = (■x,-='() 


(10.3-2) 

in which 

^ J * ^ <©j_ ” “‘'l ’ ®a- 

- c^» C«®^ H. \ ^fn@^ 

Ky c Y® <=xrt<3^-^v J S**n(5)^ 

c^(S>^ — r f 

(. 3-3) 

"’*, . "p .Jiare masses of Bp, . B, and total system, 
respectively,^^ ^ areas given in (3.3-J») in parti, 




and other quantities are related to the geometry and mass 
distribution in the system. 

The potential energy of the system is only due to the 
elastic strain energy and is, as in part 1, given by 



( 10 , 3 - 4 ) 


From the Lagrangian L« T^- , the equations of motion are 

derived, and after linearizing these equations about a nominal 

p 1 « * • * 

state 0| , , and assuming 0^ » » 0^ • ®, *■ ^ 


we obtain 
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( 10 . 3 - 5 ) 

where only a single mode for the appendage has been used, and 
•j , ^ are small perturbations about , respect ivel y, 

is a damping coefficient for the appendage, 





(10.3-6) 
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Ik 


z ^ — 2l£. 
tH 

and p * Pj » f| system constants based on mass distribution 

and geometry. Now, let the torques . 7J of the PID 

controllers be 

7 ^ = - [ + *^ 11 . ®| ‘^13 } 

• t n 

r r - [ Kil j 

Taking the transform of (10.’-'') and (10.;i-7) into the 
s-domain (complex domain), we obtain 
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The characteristic equation for the system is obtained by 
equating the determinant of the coefficient matrix in this 
equation to zero. The result is 







where 


( 10 . ^- 9 ) 


f, • 

■f, . J -)- e, ^ 

= Cj *’ + ^ a' 

fj = +<=3 

fy * «tJ + =y :» +S/r'- 

• S ^ + S/*'- 

S ^ + Cj/*'- k,^ sy. ),,*• 

f, ' - ifi <f ‘ h K 'f/''" ^ 


f / 0 - 3 - » 0 ^ 

\ 
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in which 


Cj « k,j 

c * k . K *+ » 

s a k„ Ka., + *<|^ *'>•*> 

S * 2Lf>^ +i)^ k,3 + •^H J^II.*** 

Cj * ^hj •'»! ■'■ ’^'2-'^'^^ 

S “ •’« " 

S = A» - * 
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The above* character i st ic equation is for a three-body system 
The CE ^or a two-body S''Stem in which the IPM girrbal is 
locked to the shuttle so that 9^ » 0, and the PID controller 
in the IPM is absent ( contro gains Ky 's for 9^ are zero} ^ 
;s obtained by striking out the second row and column in the 
coefficient matrix in (10. ^-8) as follows: 

i|*3 kn |)| » 

|>j j** 



(10.3-12) 


from which the CE for a two-body system is 


I 



c^e^O 


f • S 
4 


where 
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*fi = *h ^ 


(10.3-13) 
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At this point since we have obtained the CE in the s-domain 

we can apply the parameter plane analysis. We will consider 

a two-body system and a three-bod” system seoarately. 





cxes< / 


! f> 


i J 


10.4 Two- Body System 

Once the CE is obtained we have to identify the parameters 
and p . if we write the equation fc*- the appendage modal 
coordinate from (10.T-5), we obtain 


(10.4-1) 


Def in in,g 


1, 




i.jTy' 

S. - 

A' 

we can write (10.4-1) as 


(10.4-2) 




(10. 4-?) 

which is a standard form. Hence, dividing (10.3-13) by , 
we can write the CE as 


I fi 5 » O 

j-0 


(10.4-4) 


where 




■f j. m I— ^ • X s- 


(10.4-5) 
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Now, from the three instrument characteristics ^ , Tfr -id 
we can identify two as the , p parameters. Usually, 
is assumed, henc»>, and ^ are the candidate? for the «< , 
I® parameters. The rest of the quantities in fj are 
to be chosen including . The control gains a^ , , a^ 

can be chosen ^n any appropriate manner bearing In mind 
that ’■he choice of these should give a stable system if the 
appendage were rigid, i.e., the system 


® p X (10.4-6) 

should be a stable system for any choice of the gains. 
Following the procedure for the parameter pla.ie analysis we 
recast f- in terms of <n and as follows: 

j f j 


i — 


-f . - -K>i o-p Ys 1 f 

Table 7 shows the coefficients «<; , A.- , Y- 

obtain the two algebraic equations 

u 


X. — 


(10.4-7) 
J- . Then 




* A, L 


' ' T (10.4-8) 

Solution of these equations for \arious values of "J , cj^ gives 
and ^ Figures 14 through 17 show these parameter 
plane curves. Figure 14 shows the unstable region; Fig. 15 
shows the ^-curves for 5* 0, 0.1,...., 1; Fig. 16 shows 
a few curves of Fig. 15 for clarity, and Fig. 17 shows the 
effect of varying the . As decreases the unstable 
region expands. 
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To Illustrate the use of these curves., refer to Fig. 16. 
If a designer wants to keep ^^0.6 and damping ratio S 
between J = 0.1 and 0.2 he could, for instance, work in 
the elbow region shaded in the diagram. By adjusting the 
system parameters ( essentially the instrument characteristics) 
he could obtain point A or B if he so chooses. 

For the two bodv system considered we see the parameter 
olane lending itself completely to the designer's choice of 
the system parameters. We assumed only the control gains 
and the instrument in the whole system in plotting 

these curves. All other system character i st i cs especially 
rhe flexible Instrument can be varied over a wide range of 
cho i ce , 
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10.5 Three- Body System 
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The characteristic equation for the three-body system is 
given in (10.3-9) through (10.3-11). As in the case of l 
two-body system we will try to identify two parameters , 

The equation for the appendage modal coordinate from (10.3-5) 
i s 




the equation 





r' 


(10,5-1) becomes 


( 10 . 5 - 2 ) 




+ 






(10.5-3) 


This equation is similar to the case of a two-body system 
(see (10.4-3)) except that there is an additional term on 
the right side due to the added motion in the system. This 
indicates that the appendage motion is controlled by both 
motions and 0^^ . 

Examination of the coefficients f* in (10.3-10) tells 
us that fj cannot be expressed in terms of only two flexi- 
bility^ parameters o( , and of the instrument, and 

the control gains. There are additional terms which would 






not disappear. Thus, for the three-body system the coeff- 

1 

Icients do not simplify as in the case of a two-body system. j 

Nevertheless, a parameter plane study can be carried out for | 

any chosen pair of parameters. All quantities except the 1 

'I 

two parameters must then be specified. The usefulness of ] 

the curves depends on what and how many 'quantities are 

1 

specified in the system. ideally we would like to leave i 

the entire flexible instrument unspecified (except Sf ) • 

In the following three different forms are shown for writing 
the fj in terms of some possible parameters. In all cases 
the CE is the same, i.e., (10.3-9). 

1 ) Form One 

fi • 

•fr + 

Yy 


4 





in which <Sc is one parameter and the normalized mass of 

/ (M 

the instrument is the other p^irameter. The are 

functions (polynomials) of the normalized mass. All other 
physical attributes except :hese parameters have to be specified. 
If the f. are written in the form 

j 


where 
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(10.5-5) 


(10.5-6) 


then the (3- , etc. are given in Tables 8 and 9. 

'j * J J 

Table 8 is for the second form and Table 9 is for the third 


form lor these f. of the CE. 
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In exactly a similar manner as in the case of a two-body 
system the control gains are chosen to stabilize a full three 
rigid-body system. Since 0 5 N S 1 the rigid system can be 
stabilized such that it is stable for any value of N. The value 
N depends on only mass distribution. In the forms two and three 
-a or is an additional parameter which is dependent on the 
parameter h’ . Hence, it is not very useful to plot planar 
curves of 0^ and S* . 

In the present work form one is used to plot the parameter 
plane curves. All quantities specified are normalized with 
respect to some system quantities. Figures 18 and 19 show 
these curves. The unstable region is shown shaded. In Fig. 18 
the effect cf varying is also seen. Due to the existence 
of multiple roots there are multiple “^-curves. Figure 19 
shows some of these ^-curves. 





10.6 Summary and Conclusions 

In the case of a two-body system the parameter plane 
proved to be a very elegant and powerful method to study 
the flexibility characteristics of the class of instruments 
that may be fitted to the shuttle of a mul t i -purpose 
spacecraft. Once the model is specified and the control gains 
are obtained in terms of ratios to some system quantity, the 
parameter plane curves can be plotted for any given value 
of the modal damping . No other data need be specified. 
This offers plentv of freedom in designing the flexible 
instrument. Also, the physical characteristics of the 
shuttle and the IPM can also be adjusted. Indeed the 
parameter plane curves are very useful to the designer 
if the system is represented by a two-body model. 

In the case of a three-body system complications arise 
and separating the flexibility characteristics of the instr- 
uments into two parameters is not possible. There must be 
three parameters (excuding ) to portray the flexibility 
of the instruments. But two of them are dependent through 
the physical characteristics of the instruments. Hence, 
even a three dimensional plot is not very useful. One 
could plot sets of two parameter plane curves for a wide 
range of the third parameter. But this may not be convenient 
or reasonably straightforward. The alternate approach is to 





1 


specify as much data as is essential regarding the shutt.ie 
and the IPM and try to express the flexibility characteristics 
that are most variable from the designer's point of view 
as parameters. This will require specification of some 
physical attributes of the instrument itself. This undoubt- 
edly lessens the freedom of design of the instruments. Some 
ingenuity is called for In choosing the proper parameters 
so that least number of data about the instrument need only 
be specified and the curves obtained are most useful in 
direct application «n designing the rest of the charactef-.- 
istics of the instruments. In short, three- body system 
leads to limitations in the use of the parameter plane curves. 
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